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APRESENTAÇÃO 


A “Engenharia Agronômica: Ambientes Agrícolas e seus Campos de Atuação” é 
uma obra que apresenta dentro de seu contexto amplas visões que reflete em ambientes 
agrícolas e seus campos de atuação trazendo inovações tecnológicas e sustentáveis que 
proporciona em melhorias sociais, ambientais e econômicas para toda comunidade agrária. 

A coleção é baseada na discussão científica através de diversos trabalhos que 
constitui seus capítulos. Os volumes abordam de modo agrupado e multidisciplinar 
pesquisas, trabalhos, revisões e relatos de que trilham nos vários caminhos da Engenharia 
Agronômica. . 

O objetivo principal foi apresentar de modo agrupado e conciso a diversidade e 
amplitude de estudos desenvolvidos em inúmeras instituições de ensino e pesquisa do 
país. Inicialmente são apresentados trabalhos relacionados a sustentabilidade, envolvendo 
questões agroecológicas, produção orgânica e natural, e suas relações sociais. Em 
seguida são contemplados estudos acerca de inovações tecnológicas do meio rural, que 
abrange qualidade de sementes, nutrição mineral, mecanização, genética, dentre outros. 
Na sequência são expostos trabalhos voltados à irrigação e manejo do solo, envolvendo 
processos hídricos, sistemas agroflorestais e adubação. 

A obra apresenta-se como atual, com pesquisas modernas e de grande relevância 
para o país. Apresenta distintos temas interessantes, discutidos aqui com a proposta de 
basear o conhecimento de acadêmicos, mestres, doutores e todos que de algum modo se 
dedicam pela Engenharia Agronômica. Abrange todas regiões do país, valorizando seus 
diferentes climas e hábitos. 

Inicialmente são apresentados trabalhos relacionados a sustentabilidade, 
envolvendo questões agroecológicas, produção orgânica e natural, e suas relações sociais. 
Em seguida são contemplados estudos acerca de inovações tecnológicas do meio rural, que 
abrange qualidade de sementes, nutrição mineral, mecanização, genética, dentre outros. 
Na sequência são expostos trabalhos voltados à irrigação e manejo do solo, envolvendo 
processos hídricos, sistemas agroflorestais e adubação. 

Assim a obra Engenharia Agronômica: Ambientes Agrícolas e seus Campos de 
Atuação expõe um conceito bem fundamentado nos resultados práticos atingidos pelos 
diversos educadores e acadêmicos que desenvolveram arduamente seus trabalhos aqui 
apresentados de modo claro e didático. Sabe-se da importância da divulgação científica, 
portanto ressalta-se também a organização da Atena Editora habilitada a oferecer uma 
plataforma segura e transparente para os pesquisadores exibirem e disseminarem seus 


resultados. 


Tamara Rocha dos Santos 
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RESUMO: A bacia do Salgado insere-se em 
uma região semiárida do estado do Ceará, 
Brasil. Essa constitui importante região 
produtora de alimentos e enfrenta restrição 
na disponibilidade de recursos hídricos. 
Neste estudo, o desempenho da agricultura 
irrigada foi analisado utilizando-se critérios 


Engenharia Agronômica: Ambientes Agrícolas e seus Campos de Atuação 3 


de segurança hídrica, produtiva, econômica e 
social. Baixos índices de desempenho sempre 
estiveram associados a inadequação de mais 
de um indicador de segurança. De uma forma 
geral os resultados indicaram a necessidade 
de redução do volume de água aplicado pela 
irrigação, aumento da produtividade, da renda e 
de empregos. Culturas ocupando a maior parte 
da área cultivada (e.g., arroz, banana, cana- 
de-açúcar, coco verde, goiaba, maracujá, milho 
verde) requerem volumes elevados de água para 
irrigação. Para essas culturas, são necessárias 
intervenções técnicas relativas ao manejo do 
sistema solo-água-planta visando a melhoria da 
produtividade da água. Em condições de fortes 
restrições hídricas e grande competição pelo uso 
da água, como é o caso da bacia do Salgado, 
sugere-se que não sejam dadas prioridades aos 
cultivos irrigados que apresentem baixo índice 
de desempenho, por exemplo, inferior a 0,70. O 
estudo aqui apresentado poderá dar suporte às 
políticas públicas no âmbito da irrigação e das 
técnicas agronômicas necessárias à melhoria do 
desempenho da agricultura irrigada. 
PALAVRAS-CHAVE: Coeficiente de uso 
consuntivo da irrigação, produtividade física da 
água, produtividade econômica da água. 


IRRIGATION WATER PRODUCTIVITY IN A 
SEMI-ARID REGION — A CASE STUDY IN 
THE SALGADO BASIN — CE, BRAZIL 
ABSTRACT: The Salgado Basin belongs to 
the semi-arid region of Ceará State, Brazil. It is 
an important food producing region that faces 
constraints of water resources. In this study, the 
performance ofirrigated agriculture was analyzed 
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based on water, productivity, economic and social security. Low levels of performance always 
were related to more than one index of security. Overall, results indicated that is necessary 
to reduce the volume of water for irrigation, to increase yield, incomes, and employment. 
Crops occupying most of cultivated land (e.g., rice, banana, sugarcane, coconut, guava, 
passion-fruit, corn), require high volumes of water for irrigation. For these crops, technical 
interventions related to management of the soil-water-plant system are necessary to improve 
water productivity. Under severe shortage of water resources and disputes among water 
users, which applies for the Salgado River Basin, irrigated crops of low performance index 
(e.g., smaller than 0.70) should not be prioritized. This study may support public policies 
related to irrigation and agronomical techniques necessary to improve the performance of 
irrigated agriculture. 

KEYWORDS: lrrigation consumptive use coeficiente, physical water productivity, economic 
water productivity. 


11 INTRODUÇÃO 


A bacia do rio Salgado é uma das cinco sub-bacias do rio Jaguaribe e se posiciona 
na porção meridional do Estado do Ceará. Possui uma extensão de 308 km e drena uma 
área de 12.623.89 km?, o equivalente a 9% do território cearense. Depois da bacia do Baixo 
Jaguaribe, é uma das mais importantes regiões produtora de alimentos, sob irrigação, 
do estado do Ceará. Possui baixa capacidade de acumulação de águas superficiais 
(447.210.000 mº). A disponibilidade de águas subterrâneas é 127.600.758 mº. O clima é do 
tipo BSW'h' pela classificação de Kôppen, quente e semiárido, com grande variabilidade na 
distribuição espacial das chuvas (INESP, 2009). 

Historicamente a escassez hídrica no Estado do Ceará foi justificada apenas por 
condições eminentemente climáticas. Atualmente, porém, deixou de ser exclusivamente um 
fator natural para ser, também, consequência da crescente demanda pelos múltiplos usos, 
dentre os quais a agricultura irrigada, turismo, crescimento populacional e indústria, são 
dependentes de grandes quantidades de água. Melhorar a gestão dos recursos hídricos na 
agricultura é objetivo prioritário em todo o mundo (HUANG e ZHANG, 2020). 

Os irrigantes estão sob pressão crescente para garantir a segurança alimentar 
bem como a sustentabilidade ambiental e econômica de longo prazo (KOECH; LANGAT, 
2018). A disponibilidade de água afeta diretamente a produção de alimentos e o suprimento 
suficiente de água é vital para garantir o crescimento das culturas e a sobrevivência do 
homem (FERNÁNDEZ et al, 2020). Esse desafio é maior nas regiões áridas e semiáridas, 
onde a escassez é o fator limitante mais importante para a produtividade das plantas. 
Essas condições podem se tornar mais severas pelo aquecimento global (SUN et al., 
2017). Pela falta de recursos hídricos nesses agro-ecossistemas, grande parte das áreas 
cultivadas não é irrigada e a agricultura é manejada sob condições de sequeiro (HADEBE 
et al., 2017), com baixa produtividade. 
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Aumentar a eficiência do uso da água na fazenda requer decisões sábias sobre 
o sistema de irrigação, a estratégia de irrigação e o método para programar a irrigação, 
entre outros fatores relacionados ao gerenciamento da água. Desde o início dos anos 
2000, a abordagem da produtividade da água tem sido amplamente usada para resolver 
esse problema. Ele fornece indicadores úteis para a produtividade física da água e o 
desempenho econômico da irrigação. A análise da literatura, entretanto, mostra confusão 
no uso dos termos e falta de acordo nas equações (HEYDARI, 2014). Este artigo é focado 
no uso racional da abordagem de produtividade da água para o irrigante para melhorar a 
produtividade física e econômica da água em escala de campo, na bacia do rio Salgado, 
no estado do Ceará. 


21 MATERIAL E MÉTODOS 


O coeficiente de uso consuntivo da irrigação (ICUC) é sugerido para definir a fração 
da água de irrigação aplicada a um campo, fazenda ou projeto que é convertida em vapor 
ou consumida (transpiração mais evaporação da superfície do solo ou das plantas) (BURT 
et al. 1997). O são medidas físicas de uma dada tecnologia de irrigação presumindo um 
nível de manejo e, portanto, não é comparáveis ao termo eficiência de uso de água ou 
produtividade da água. A fração não consumida é 1- ICUC, representa a porção recuperável. 

Em qualquer estudo de balanço hídrico de um projeto ou bacia hidrográfica ou ao 
estimar o impacto de alguma intervenção, deve-se considerar ambos os indicadores, a 
fração consumida e a não consumida (JENSEN, 2007). O indicador ICUC é apropriado ao 
considerar a água consumida (evapotranspiração da cultura — ET) na produção do efeito 
desejado (produção da cultura), mas é um termo inapropriado se considerar a água não 
consumida como desperdiçada, uma vez que frequentemente essa água é recuperada e 
reutilizada em escala de bacia hidrográfica (GRAFTON, et al. 2018). 

Nos sistemas de produção agrícola a WFP é usada para definir a relação entre 
produção comercial das culturas e a quantidade de água consumida nesta produção (ET), 
e tem servido como um indicador para quantificar o impacto dos calendários de irrigação 
com relação ao manejo da água. Dessa forma, a produção total de biomassa (matéria seca) 
é transformada em produção comercial da cultura (VAZIFEDOUST et. al. 2008; PERRY et 
al. 2009) e define-se WFP com respeito à ET (WFP..), conforme Eq. (1). AWFP,., constitui 
a chave para avaliação das estratégias de irrigação com déficit (GEERTS e RAES, 2009). 


Rendimento comercial da cultura (kg) 


WFPpr = 
Br Evapotanspiração da cultura (m?) 


(1) 


Na escala de campo, o uso de água representado no denominador das Eg. (1) é 
difícil de ser determinado com precisão. Assim, em algumas situações, outros substitutos 
para WFP são usados por muitos profissionais de irrigação e, como consequência, resultam 
diferentes valores. Se a quantidade total de água aplicada [irrigação (IR) + precipitação 
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efetiva (PE)] é considerada como água consumida pela cultura, então a Eg. (1) pode ser 
Eq. (2)]. O denominador 
da Eq. (2) é um substituto do uso de água para se obter o rendimento comercial 


usada para determinar a produtividade física da água (WFP p..e) [ 


correspondente. Sob condições de muito baixa precipitação, como em regiões áridas e 


semi-áridas, pode-se converter WFP em WFP,, (Eq. 3). Nesses casos negligencia-se 


IR+PE 
a variação do armazenamento de água no solo durante o ciclo da cultura, a percolação, a 
ascensão capilar e o escoamento superficial. Muitos profissionais usam as Egs. (2) e (3) 
com a finalidade de identificar diferenças entre os métodos de irrigação e/ou manejo de 


irrigação. 


Rendimento comercial da cultura (kg) 
WFPr+pe = (2) 


Volume de irrigação + precipitação (mº) 


Rendimento comercial da cultura (kg) 
WFPjp =——— >>> > >—>—>——>——— 


(3) 


Observa-se relativo consenso sobre o numerador de WFP ser o rendimento 


Volume de irrigação (m?) 


comercial. A biomassa total seca ou fresca ou produto colhido pode ser usado no 
numerador, expresso em termos físicos ou econômicos. No entanto, como os valores 
econômicos de diferentes produtos agrícolas não são os mesmos, a produtividade da água 
deve ser definida economicamente (HEYDARI, 2014). Por outro lado, algumas culturas 
são de base úmida, outras de base seca, algumas culturas como alfafa e milho podem ser 
reportadas em ambas as bases e para algumas outras o produto final é mais importante (por 
exemplo, cana-de-açúcar, beterraba açucareira, arroz em casca, etc.). Já a definição e os 
componentes do denominador diferem entre os autores. Alguns consideram apenas a água 
transpirada pela cultura (BOUMAN, 2007), mas a maioria inclui tanto a água evaporada 
do solo quanto a transpirada pela cultura. Pereira et al. (2012) sugeriram o uso do volume 
total de água envolvido na produção, isto é, além da ET somam-se a fração de lixiviação 
quando houver risco de sal na zona radicular e o uso não benéfico de água, ou seja, a água 
que é percolada, escoada do local de cultivo e derivada pelo vento quando a irrigação é por 
aspersão (isto é, volume total de agua aplicado). 

O conceito de produtividade da água pode se referir à produção física ou ao seu valor 
econômico (RODRIGUES e PEREIRA, 2009; HEYDARI, 2014). Um termo adequado para o 
último é a produtividade econômica da água (WEP), que se refere à razão entre produtos 
e insumos em termos monetários. Para culturas com baixos custos de investimento, por 
exemplo cereais, a produtividade econômica bruta da água de irrigação (WEPB,), que 
considera a margem bruta (receita bruta menos os custos variáveis) é aceitável [Eq. (4)] 
(FERNÁNDEZ et al., 2020): 


Margem bruta ($ ha!) (4) 


WEPBp =———————————————— 
IR Volume de irrigação (m3ha-?) 
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em que o denominado representa apenas o uso da água aplicada por irrigação. 

No caso de culturas lenhosas e outras culturas que requerem substancial 
investimento desde o início, a produtividade econômica líquida da água (WEPL,) é um 
indicador mais adequado (FERNÁNDEZ et al., 2020), pois considera a margem líquida em 
vez da margem bruta, ou seja, inclui os custos variáveis e os fixos [Eqg. (5)]: 


Margem líquida ($ ha?) 


WEPLjp =» —>>»>»>>—>—>—>—»>———> 
IR Volume de irrigação (mºha-!) 


(3) 


Ainda assim, nem WEPB, nem WEPL, consideram os custos de oportunidade, 
definidos como os benefícios perdidos ao longo da vida útil da cultura a uma determinada 
taxa de juros (FERNÁNDEZ etal., 2020). Uma análise econômica adequada deve considerar 
os custos de oportunidade. Portanto, a produtividade econômica da água total (WEP q.pe); 
definida pela Eq. (6) é recomendada: 


lucro ($ ha"5) 


WEB, >|] DWDWw > —>—>w>w>———»>»>——— 
IR+PE " Volume de irrigação + PE (m3ha-!) 


(6) 


Se apenas o volume de irrigação for considerado no denominado, tem-se a 
produtividade econômica da água de irrigação (WEP ,). 
O lucro no numerador da Eq. (6) é definido como a receita bruta menos a soma dos 


custos variáveis, fixos e de oportunidade. Então WEP e WEP,, são adequados para 


IR+PE 
tomar decisões sobre o manejo de irrigação de culturas lenhosas (FERNÁNDEZ et al., 
2020). Ressalta-se que o tempo de vida da cultura afeta o valor WEP resultante, devido 
ao seu impacto nos custos fixos e de oportunidade. Isso é particularmente relevante para 
pomares de árvores frutíferas. Tanto WEP ..,, quanto WEP,, são particularmente úteis 
para irrigantes que precisam tomar decisões sobre como gerenciar a irrigação da maneira 
mais lucrativa, ou seja, quando a meta de produção é aumentar a lucratividade e não a 


produtividade física da água. Entretanto, um cálculo preciso do WEP , pe € WEP , não pode 


IR+PE 
ser feito antes do final da colheita, quando a receita e os custos ainda não são conhecidos. A 
receita é dada pelo rendimento e pelo valor de mercado e, para o cálculo dos custos totais, 
devem ser conhecidos os custos fixos, variáveis e de oportunidade (ALCON et al., 2013). 


Isso limita a utilização de WEP e WEP,, para a tomada de decisão na irrigação, uma 


IR+PE 
vez que a avaliação econômica deve ser feita antes do início da temporada de irrigação. O 
desafio é maior quando o valor do rendimento depende da qualidade do produto e quando o 
preço de alguns insumos, como energia, fertilizantes e defensivos, varia de uma temporada 
para outra e mesmo durante o período de crescimento. 

Indicadores de desempenho da agricultura irrigada são parâmetros importantes para 
se analisar o impacto da irrigação sobre a produção agrícola, receita líquida e geração de 
empregos em uma região. Para analisar o desempenho das culturas irrigadas são propostos 
neste estudo indicadores de desempenho agrupados em quatro classes (Tabela 4), cada 


uma composta por dois indicadores: (a) segurança produtiva — indicadores: produtividade 
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da terra (kg ha!) e produtividade da água (kg m?), (b) segurança econômica — indicadores: 
produtividade econômica da terra (R$ ha”) e produtividade econômica da água aplicada 
(R$ mº), (c) segurança social — indicadores: quantidade de empregos gerados por unidade 
de área (empregos ha!) e por unidade de volume de água aplicada (empregos mº) e (d) 
segurança hídrica — quantidade de água utilizada na irrigação por unidade de área (mº ha!) 
e ciclo da cultura. 

Para formular um índice geral de desempenho da agricultura irrigada (|) foram 
definidas as seguintes relações adimensionais: 

(a) Segurança hídrica: os dois indicadores são: 

(a) Coeficiente de uso consuntivo da irrigação (ICUCij): 

ICUC;; = A (7) 

sendo VR, o volume líquido de irrigação requerido (mº ha!) pela cultura i na sub-bacia 
je VA, o volume bruto de irrigação aplicado (mº ha”). Para culturas irrigadas com déficit 
hídrico o valor de será maior que 1. No entanto, na análise que segue para desempenho 
das culturas irrigadas, o será assumido no máximo igual a 1. 

(a,) Ciclo da cultura: 

Culturas permanentes: C, = 0,5 

Culturas temporárias de ciclo longo (ciclo 180 dias): C, = 0,75 

Culturas temporárias de ciclo curto (ciclo <180 dias): C,= 1,0 

sendo C, o peso atribuído à cultura i, referente ao comprimento do ciclo. 

(b) Segurança produtiva: os dois indicadores são: 

(b,) Razão de produtividade física da terra (LFPR;): 

LFP,; 

LFP, 


imáx 





LFPR;; = 


(8) 


sendo LFP, a produtividade física da terra (kg ha") da cultura i na sub-bacia j, e LFP, 
máx à Máxima produtividade física da terra (kg ha!) observada na bacia para a cultura i. 
(b,) Razão de produtividade física da água (WFPR): 


baga) (9) 
WEFP; máx 





sendo WFP, a produtividade física da água (kg m?º) da cultura i na sub-bacia j, e 
WFP «à Máxima produtividade física da água (kg m*) observada na bacia para a cultura i. 

(c) Segurança econômica: os dois indicadores são: 

(c,) Razão de produtividade econômica da terra (LEPR,): 





LEB,; 
LEPRj = rp (10) 


sendo LEP,, a produtividade econômica da terra (R$ ha!) da cultura i na sub-bacia 


Engenharia Agronômica: Ambientes Agrícolas e seus Campos de Atuação 3 Capítulo 1 6 


j, e LEP, ., a máxima produtividade econômica da terra observada na bacia (R$ ha!) para 
a cultura i. 
(c,) Razão de produtividade econômica da água (WEPR): 
WEB,; 


WEPR;; É WEP; max (11) 


sendo WEP,, a produtividade econômica da água (R$ mº) da cultura i na sub-bacia 
j e WEP, 


imáx 


a máxima produtividade econômica da água (R$ m*º) observada na bacia para 
a cultura i. 

(d) Segurança Social: os dois indicadores são: 

(d,) Razão de número de empregos gerados por unidade de área (LLABR): 


LLAB;j (12) 


LLABR; =———+— 
) LLAB; máx 


sendo LLAB,, o número de empregos gerados por unidade de área (empregos ha” 
?) pela cultura i na sub-bacia j, e LLAB, 


i máx 


o número máximo de empregos gerados por 
unidade de área (empregos ha!) observado na bacia para a cultura i. 
(d,) Razão de número de empregos gerados por unidade de volume de água 
aplicado (WLABR): 
WLAB,; 
WLABR; = WLAB, ns; (13) 
sendo WLAB,, o número de empregos gerados por unidade de volume de água 


aplicada (empregos m*º) pela cultura i na sub-bacia j, e WLAB.. o número máximo de 


empregos gerados por unidade de volume de água aplicada (empregos m*) observado na 
bacia para a cultura i. 
O índice de desempenho (1) para uma cultura irrigada i na sub-bacia j foi definido da 


seguinte forma: 


lj = €,ICUC; + «> ICUC;j + a3LFPRij + qq WFPRi; + as LEPRi; + oWEPRi; + «,LLABR;; (14) 
+ as WLABRi;; 

sendo a, os pesos atribuídos pelo gestor a cada razão de segurança de acordo com 
suas prioridades; a + a, +...+ 0q,=1. Sea, =q,=..= q, =a, então q=1/8 e , corresponde à 
média aritmética das razões de segurança. O valor de |, varia entre zero e 1. Quanto mais 
próximo de 1 melhor é o desempenho da cultura irrigada e, conforme |, Vai reduzindo, 
significa que uma ou mais razões de segurança estão baixas e a atividade agrícola 
necessita de alguma intervenção de melhoria. O índice , também serve para priorizar as 
atividades agrícolas irrigadas. Para culturas irrigadas com déficit hídrico obtém-se ICUC 

maior que 1. Para normalizar |, entre zero e 1, considerou-se Os ICUC sl. 


Para o cálculo do indicador de produtividade econômica considerou-se a margem 
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líquida, isto é, da receita bruta subtrafram-se os custos fixos e as variáveis. As informações 
sobre áreas cultivadas, produtividade física da terra, volume bruto de irrigação aplicado 
aos cultivos, empregos gerados por unidade de área pelas atividades agrícolas, custos 
de produção e valor da produção foram obtidas na base de dados da Agência de 
Desenvolvimento do Estado do Ceará (ADECE, 2018) e Instituto Brasileiro de Geografia e 
Estatística (IBGE, 2017). O estudo de caso foi feito com as informações do ano de 2017. 

Utilizou-se o método de Penman-Monteith divulgado no Boletim FAO 56 (ALLEN 
et al., 1988) para cálculo da evapotranspiração de referência (ETo). Os coeficientes 
de cultura (Kc), utilizados para o cálculo da evapotranspiração de cultura (ETc) foram 
aqueles divulgados pela FAO 56 adaptados para as condições climáticas da região. Os 
dados utilizados para o cálculo da ETo na bacia do Salgado foram obtidos de estações 
meteorológicas pertencentes ao Instituto Nacional de Meteorologia (INMET), localizadas na 
região de estudo. O sistema SQI (LIMA et al., 2015) foi o SSD utilizado para a determinação 
da ETc nos municípios da bacia. O volume de irrigação requerido pelas culturas (VR) foi 
estimado com base na evapotranspiração potencial de cultura (E Tc), já que a precipitação 
pluvial na bacia no ano de 2017 foi muito baixa, isto é, em média 32,7 mm. 


31 RESULTADOS E DICUSSÕES 


Na Tabela 1 mostra-se que, em 2017, as 19 culturas irrigadas na bacia do salgado 
ocuparam uma área irrigada 4.742 ha a cultura da banana ocupou a maior área (34,4%). 
Observa-se que a cultura do feijão também se destaca com 958 ha irrigados (20,2%). A 
área total irrigada recebeu um volume bruto de irrigação de 64.977.730 mº, enquanto o 
volume requerido foi de 44.269.610 mº, resultando num ICUC global de 0,681. As atividades 
agrícolas geraram uma produção total de 120.321.056 kg, produziram uma receita líquida 
de R$ 95.191.188 e geraram 3.429,63 empregos diretos. Três culturas (banana, feijão- 
caupi e manga) ocupam 63% da área, recebem volume bruto de irrigação de 64,8% do 
volume total e produzem 57,5% da receita líquida total da bacia. A banana é o produto que 
rendeu maior receita líquida total para a bacia, entretanto a uva é a cultura que apresenta 
a maior rentabilidade por unidade de área. De acordo com informações da ADECE (2018) 
a área irrigada na bacia está distribuída em 21 municípios, tendo o município de Barbalha 
a maior área irrigada (1.069 ha). 
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Culturas Área VA VR LFP LEP Empregos 
(ha) (mê ha?) (mºha?!) (kg ha?) (R$ ha?) (por ha) 





Abacate 16 19.000 7.850 10.000 19.017 0,70 
Acerola 50 18.000 9.960 11.500 45.460 1,70 
Arroz 57 28.000 10.920 9.528 4.754 0,73 
Banana 1.630 18.000 11.950 28.034 24.704 0,54 
Batata doce 159 10.490 8.950 25.400 21.856 1,30 
Cana-de-açúcar 140 19.000 6.220 59.533 8.807 0,16 
Coco 114 15.000 8.470 12.300 11.043 0,26 
Feijão-caupi 958 7.500 7.090 3.840 5.410 0,90 
Goiaba 106 15.000 7.460 33.000 33.004 0,32 
Laranja 84 16.500 9.960 11.250 15.047 0,34 
Limão 27 14.000 9.960 9.400 20.578 0,35 
Macaxeira 353 7.200 7.970 31.000 17.879 0,36 
Mamão 100 15.000 8.460 77.325 33.569 0,55 
Manga 397 14.000 9.960 20.475 23.251 0,40 
Maracujá 209 14.000 9.950 17.300 23.063 0,38 
Melancia 52 12.000 5.810 60.200 12.599 0,70 
Milho verde 73 12.000 6.710 3.600 10.566 0,82 
Tomate 208 10.000 4.980 82.000 46.540 3,22 
Uva 9 17.600 8.470 34.625 86.333 2,19 
Total 4.742 


Tabela 1 — Culturas e respectivas áreas cultivadas, volumes brutos de irrigação aplicados (VA), 
volumes de irrigação requeridos (VR), produtividade da terra (LFP) e receitas líquidas (LEP), na 
bacia do rio Salgado. 


As estimativas de (Tabela 2) feitas aqui são aproximadas porque os volumes 
brutos de irrigação aplicados às culturas representam informações gerais apresentadas 
no relatório da ADECE (2018), revelando a necessidade de contabilidade abrangente da 
água em escala do campo e de bacia hidrográfica. Conforme salienta Burton (1999) a maior 
restrição para o entendimento sobre o uso da água na agricultura é a dificuldade associada 
com sua medida e quantificação. Medições e registros de vazão do sistema, de volume 


aplicado e variáveis de clima são raros e providos de erros potenciais. 
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Culturas ICUC WFP WEP Empregos I 





(kg mº) (R$ mº) (por 1000 mº) 
Abacate 0,413 0,579 1,001 0,037 0,793 
Acerola 0,553 0,639 2,526 0,094 0,724 
Arroz 0,390 0,340 0,170 0,026 0,776 
Banana 0,664 1,557 1,372 0,030 0,717 
Batata doce 0,855 2,426 2,084 0,132 0,982 
Cana-de-açúcar 0,327 3,133 0,464 0,008 0,797 
Coco verde 0,565 0,820 0,736 0,017 0,674 
Feijão-caupi 0,945 0,512 0,721 0,120 0,925 
Goiaba 0,497 2,200 2,200 0,021 0,695 
Laranja 0,604 0,682 0,912 0,021 0,766 
Limão 0,711 0,671 1,470 0,025 0,755 
Macaxeira 1,107 4,306 2,483 0,050 0,718 
Mamão 0,564 5,155 2,238 0,037 0,802 
Manga 0,711 1,463 1,661 0,029 0,845 
Maracujá 0,711 1,236 1,647 0,027 0,658 
Melancia 0,484 5,017 1,050 0,058 0,785 
Milho verde 0,559 0,300 0,881 0,068 0,601 
Tomate 0,498 8,200 4,654 0,322 0,845 
Uva 0,481 1,967 4,905 0,124 0,860 





Tabela 2 — Coeficiente de uso consuntivo da água da irrigação (ICUC), produtividade física da 
água de irrigação (WFP), produtividade econômica da água de irrigação (WEP), geração de 
empregos por 1000 mº de água aplicada e índice de desempenho das culturas irrigadas (|) na 
bacia do Salgado para q=1/8. Ano de 2017. 


Entre as 19 culturas irrigadas, apenas 6 (31,6%) apresentaram ICUC superior a 
0,70, indicando a necessidade de adoção de técnicas para um manejo racional da irrigação. 
A cultura da macaxeira, com ICUC superior a 1 indica que a cultura está sendo manejada 
com déficit de irrigação. Valores muito baixos de ICUC são encontrados para as culturas do 
arroz e cana de açúcar. O ICUC contabiliza a fração da água total aplicada que é utilizada 
para usos consuntivos, no caso, especialmente evapotranspiração das culturas. Não são 
contabilizados aqui outros usos consuntivos como evapotranspiração de plantas daninhas, 
evaporação da água aspergida nos sistemas de aspersão, evaporação de reservatórios 
e água exportada com o produto colhido. A fração 1 — ICUC representa a água que, de 
alguma forma, transita para fora da área de cultivo, mas não se perde para a atmosfera 
na forma de vapor. A água de uso não consuntivo pode ser reutilizada na própria área, ou 
em outra área de jusante na bacia hidrográfica ou retornar à fonte. Está relacionada com a 
água que drena abaixo da zona radicular, a água que escoa superficialmente e a que infiltra 
nos canais (BURT et a., 1997). 
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Os valores tão baixos de ilustram problemas típicos de aplicação de água. Esses 
problemas enfatizam a necessidade de considerar rigorosamente técnicas de manejo da 
irrigação, seleção de solos aptos ao cultivo sob inundação, projetos adequados dos sistemas, 
técnicas adequadas de operação dos projetos de irrigação, programa de manutenção dos 
sistemas de irrigação e treinamento dos irrigantes. Oportunidades para aumentar ICUC 
são adoção de sistemas de programação ad irrigação baseados na demanda das culturas, 
uso de programação ad irrigação com déficit, uso de técnicas de manejo e conservação 
da água do solo que permitam a redução das perdas por evaporação, uso de sistemas de 
irrigação que minimizam a fração de solo molhado, melhor manutenção dos equipamentos 
de irrigação e estruturas de condução de água (WALLACE e BATCHELOR, 1997). 

Na Tabela 2 são apresentadas as produtividades física e econômica da água para 
as culturas irrigadas na bacia do salgado, no ano de 2017. Para a produtividade física 
da água foi considerado volume de água aplicado por irrigação (VA) aos cultivos. Para a 
produtividade econômica da água considerou-se a receita líquida advinda da produção das 
culturas. Não foram computados os custos de oportunidade. 

A cultura do arroz apresenta LFP (9.528 kg ha”) relativamente boa. Mas a WFP é 
baixa (0,340 kg m*) por receber um volume bruto de irrigação alto (28.000 mº ha”). Este 
valor é inferior ao valor médio 0,43 kg mº verificado por Costa et al. (2005) com aplicação 
de um volume de água 13.760 mº ha! em solos de textura argilo-siltosa no perímetro 
irrigado Morada Nova, pertencente à sub-bacia do Baixo Jaguaribe. Observa-se ainda que 
WEP para o arroz é baixa, em média R$ 0,17 mº, porque a receita líquida por hectare (LEP) 
advinda do cultivo não é elevada e o volume bruto de irrigação aplicado é alto. Destaca- 
se ainda uma contribuição social do cultivo do arroz na região apenas regular, gerando 
durante o ano na bacia, em média, 0,73 empregos ha"! ou 0,026 empregos por 1000 mº de 
água aplicada. Estes valores levaram a um índice global de desempenho (1) regular (0,776). 
O aumento da WFP pode ser conseguido tanto pela redução do volume de água aplicado 
como pelo aumento da produtividade para a mesma quantidade de água evapotranspirada 
pela cultura (Coelho et. al. 2015). 

O problema se repete com a cultura da cana-de-açúcar, que embora apresente ligual 
0,797 necessita de redução no volume de irrigação aplicado, pois, apresenta o menor ICUC 
da bacia (0,327). Teodoro et al. (2013), conduzindo experimentação de campo no município 
de Rio Largo — AL, região Nordeste do Brasil, apontam que a produtividade de colmos de 
cana-planta de 160.000 kg ha”, utilizando irrigação por gotejamento superficial, pode ser 
obtida com aplicação de 17.370 mº ha”! de água e 30 kg ha” de nitrogênio, resultando WFP 
de 9,22 kg m. De acordo com Oliveira et al. (2011), em pesquisa realizada no município de 
Carpina — PE utilizando sistema de aspersão convencional, a irrigação plena proporcionou 
maior produtividade física da água, com acréscimo na produção média de colmos de 7,0 kg 
m% em comparação ao regime de sequeiro. 

O problema é agravado com as culturas do coco (| = 0,679), goiaba (| = 0,696), 
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maracujá (| = 0,658) e milho verde (| = 0,601), indicando um desempenho marginal da 
dessas culturas. Os resultados mostram que são necessárias ações para melhoria 
expressiva do manejo da irrigação e adoção de práticas agronômicas que promovam 
aumento da produtividade da terra e uso de água de irrigação. Em condições de fortes 
restrições hídricas e grande competição pelo uso da água, como é o caso da bacia do 
Salgado, sugere-se que não sejam dadas prioridades aos cultivos que apresentem baixos 
índices de desempenho, por exemplo, inferiores a 0,70. Ou, de outra forma, que seja dado 
suporte técnico ao produtor rural para melhoria do desempenho da irrigação. O sucesso 
da aplicação de técnicas agrícolas modernas depende também de ações institucionais 
que visem estimular a assistência técnica e a difusão do conhecimento, a educação e a 
formação do irrigante e a promoção de incentivos para o uso eficiente da água e penalidades 


para o uso ineficiente. 


41 CONCLUSÕES E RECOMENDAÇÕES 


Neste artigo analisou-se desempenho da agricultura irrigada na bacia do rio 
Salgado, no estado do Ceará. Foram analisados critérios de segurança hídrica, produtiva, 
econômica e social. Baixos índices de desempenho sempre estiveram associados a 
inadequação de mais de um indicador de segurança. De uma forma geral os resultados 
indicaram a necessidade de redução do volume de água aplicado pela irrigação, aumento 
da produtividade, da renda e de empregos. Várias culturas com importantes áreas cultivadas 
e, portanto, demandando elevados volume de irrigação, como arroz, banana, cana-de- 
açúcar, coco verde, goiaba, maracujá e milho verde necessitam de intervenções técnicas 
relativas ao manejo do sistema solo-água-planta visando a melhoria da produtividade com 
menos água aplicada. Em condições de fortes restrições hídricas e grande competição 
pelo uso da água, como é o caso da bacia do Salgado, sugere-se que não sejam dadas 
prioridades aos cultivos irrigados que apresentem baixos índices de desempenho, por 
exemplo, inferiores a 0,70. O estudo aqui apresentado poderá dar suporte às políticas 
públicas no âmbito da irrigação e das técnicas agronômicas necessárias à melhoria do 
desempenho da agricultura irrigada. 
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CIDADE DE GOIÁS (GO) 


RESUMO: A determinação de atributos de chuvas 
para determinada localidade permite planejar 
práticas de prevenção e controle da erosão. O 
objetivo do presente trabalho foi determinar a 
distribuição da erosividade, padrões hidrológicos 
e período de retorno das chuvas máximas diárias 
anuais na cidade de Goiás, GO. Foram utilizadas 
séries de dados pluviográficos referente ao 
período de 2002 a 2011, fornecidos pelo INMET, 
e série de dados pluviométricos do período de 
1949 a 2011, disponível na base HIDROWEB/ 
ANA. Os pluviogramas foram digitalizados e 
seus dados processados em planilha eletrônica 
e, posteriormente, submetidos ao software 
Chuveros a fim de determinar os índices de 
erosividade El, e KE>25, altura precipitada, 
quantidade de chuvas e respectivos padrões 
de precipitação. Os resultados encontrados 
demonstraram que a precipitação média anual 
foi de 1.648,6 mm e o período chuvoso vai de 
outubro a abril, com 90,8% da precipitação anual. 
Os valores médios para os índices de erosividade 
El, e KE>25 foram de 7.112,22 MJ mm ha” h' 
ano! e 110,3 MJ ha ano”, respectivamente. O 
padrão de chuvas avançado (49,8%) apresentou 
maior ocorrência, seguido dos padrões atrasado 
(35,1%) e intermediário (15,1%). Os valores de 
chuva máxima diária anual para a cidade de 
Goiás, esperados para períodos de retorno de 2, 
5, 10, 20, 50 e 100 anos são, respectivamente, 
de 93,9; 112,6; 131,0; 154,2; 171,5 e 188,6 mm. 
PALAVRAS-CHAVE: Erosão hídrica, fator R, 
enchentes, chuva de projeto. 


Capítulo 2 


HYDROLOGICAL CHARACTERISTICS OF RAINFALL TO SOIL AND WATER 
CONSERVATION AT GOIÁS CITY (GO) BRAZIL 


ABSTRACT: The determination of rainfall attributes for a given locality allows planning 
erosion prevention and control practices. The objective of the present work was to determine 
the distribution of erosivity, hydrological patterns and period of return of the annual maximum 
daily rainfall in the city of Goiás, GO. A series of pluviographic data for the period from 2002 to 
2011, provided by INMET, and a series of rainfall data from the period 1949 to 2011, available 
at the HIDROWEB/ANA database, were used. The pluviograms were digitized and their 
data processed in a spreadsheet and later submitted to the Chuveros software in order to 
determine the erosivity indexes El, and KE>25, precipitation height, amount of rainfall and 
respective precipitation patterns. The results showed that the average annual precipitation 
was 1,648.6 mm and the rainy season goes from October to April, with 90.8% of the annual 
precipitation. The averages for the El, and KE>25 erosivity indices were 7,112.2 MJ mm ha” 
h' yr! and 110.3 MJ ha” yr, respectively. The advanced rainfall pattern (49.8%) presented 
the highest occurrence followed by delayed (35.1%) and intermediate (15.1%) patterns. The 
annual maximum daily rainfall values for the city of Goiás, expected for return periods of 2, 5, 
10, 20, 50 and 100 years are, respectively, 93.9; 112.6; 131.0; 154.2; 171.5 and 188.6 mm. 
KEYWORDS: Water erosion, R factor, floods, rainfall of project. 


11 INTRODUÇÃO 


As perdas de solo e água por erosão hídrica proporcionam a diminuição da 
produtividade das culturas, além de causar problemas ambientais, como o assoreamento 
de rios e poluição dos recursos hídricos (CASSOL et al., 2008). 

A quantificação de fatores relacionados à erosividade das chuvas são informações 
cruciais para adequar sistemas de manejo de solos, adoção de práticas conservacionistas, 
auxiliar planos de recuperação de áreas degradadas, além de poder subsidiar políticas 
públicas e assim, permitindo o uso de recursos financeiros para fins conservacionistas, 
priorizando as regiões mais críticas (OLIVEIRA et al., 2011; SILVA et al., 2008). 

Nesse contexto, a equação universal de perda do solo (USLE) é um dos modelos mais 
utilizados para estimar a quantidade média de perda de solo a longo prazo, mensalmente 
e/ou anualmente (LEE e HEO, 2011), e ajustar o sistema de produção agrícola para que as 
perdas de solo fiquem dentro de limites toleráveis. A USLE foi desenvolvida por Wischmeier 
e Smith (1965, 1978) nos Estados Unidos, a qual é composta por fatores erosividade das 
chuvas — R; erodibilidade do solo — K; comprimento da inclinação — L; e inclinação — S; uso 
e cobertura — C; e práticas conservacionistas — P (MEUSBURGER et al., 2011). 

O El, é o índice preconizado pela USLE e um dos mais utilizados para a 
quantificação do fator erosividade, sendo a multiplicação da energia cinética (E) pela 
intensidade máxima da chuva durante um período de 30 min (1,,) para cada tempestade 
de chuva (WISCHMEIER e SMITH, 1958, 1978). Para determinação do fator R da USLE 
necessita-se de séries históricas de dados de pluviógrafos de pelo menos 22 anos. Séries 
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temporais com mais de 20 anos são escassas no Brasil, logo, modelos que relacionam 
o El, com dados pluviométricos diários vêm sendo utilizados para se obter o fator R da 
USLE (ANGULO-MARTÍNEZ, 2009; MACHADO et al., 2013). Outro índice também utilizado 
e recomendado para regiões tropicais é o KE>25, o qual considera a energia cinética de 
chuvas com intensidade de precipitação superior a 25 mm h'! (HUDSON, 1973). 

O conhecimento do índice de erosividade tem importância tanto prática quanto 
científica, já que do ponto de vista prático, esse índice permite a identificação do período do 
ano com maior risco ou potencial de erosão hídrica, ou seja, permitindo o planejamento de 
medidas de controle. Já no ponto de vista científico, esse índice permite saber o motivo das 
diferenças nas quantidades de erosão em determinado ano e local (SCHICK et al., 2014). 

No Brasil, o parâmetro El, vem sendo amplamente estudado e utilizado 
(ALBUQUERQUE et al., 2005; MACHADO et al., 2017; PINHEIRO et al., 2018; ALMEIDA 
et al., 2019). Segundo Oliveira et al. (2012), o índice de erosividade El, no Brasil varia 
entre 1.672 e 22.452 MJ mm ha”! h! ano”, enquanto em países como a Europa, Grécia e 
Suíça, foram encontrados índices entre 51,4 a 6.228,7 MJ mm ha” h! ano! (PANAGOS 
et al., 2015a), 84 a 2.825 MJ mm ha” h” ano! (PANAGOS et al., 2015b) e 124 a 5.611 
MJ mm ha! h! ano” (MEUSBURGER et al., 2011) de erosividade, respectivamente. Essa 
grande diferença entre os valores pode ser atribuída à extensão dos territórios e variações 
climáticas dos locais (OLIVEIRA et al., 2012). 

De acordo com sua intensidade, as chuvas podem ser classificadas em padrão 
avançado, intermediário, e atrasado, quando a localização dos picos de maior intensidade 
está no início, meio e fim do período de duração das chuvas, respectivamente (HORNER 
e JENS, 1941). Conforme Panachuki et al. (2006), maiores perdas de solo são previstas 
para os padrões intermediário e atrasado. Isso ocorre devido à maior umidade antecedente 
do solo que favorece o processo de desagregação pelo impacto das gotas de chuvas, 
ocasionando selamento superficial e, posteriormente, redução da taxa de infiltração de 
água e aumento do escoamento superficial. 

O conhecimento das chuvas máximas diárias anuais associadas a períodos de 
retorno é uma das formas de obter a chuva de projeto, sendo também ferramenta bastante 
útil no dimensionamento de obras hidráulicas e conservacionistas como terraços agrícolas 
e bacias de infiltração (barraginhas), a fim de permitir o dimensionamento adequado dessas 
e outras práticas mecânicas (ELTZ et al., 2013). 

Acidade de Goiás (GO) localiza-se em região suscetível a processos erosivos devido 
a problemas com infiltração de água das chuvas/escoamento superficial. Esse problema é 
resultante do desmatamento, uso e manejo inadequado do solo, relevo predominantemente 
acidentado à montante, grande área de afloramentos rochosos e solos rasos, decorrentes 
de geomorfologia complexa. Isso vem acarretando problemas de cheias históricas nos 
afluentes e no rio Vermelho, com diversos prejuízos na área urbana incluindo o patrimônio 
histórico tombado pela UNESCO (CAVALCANTI et al., 2008; SANTOS, 2014). 
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Diante do exposto, caracterizar as chuvas da cidade de Goiás permite obter 
parâmetros que contribuam para o incentivo ao planejamento e uso de práticas 
conservacionistas visando reduzir a erosão em áreas rurais nos seus diferentes usos 
da terra, de forma a contribuir para redução de enchentes na bacia do rio Vermelho, e 
consequentemente, nos impactos na área urbana da cidade de Goiás (GO). O presente 
trabalho teve por objetivo determinar a distribuição da erosividade anual e do fator R da 
USLE, dos padrões de precipitação pluvial, e da chuva máxima diária anual associada a 
períodos de retorno para a cidade de Goiás — GO. 


21 MATERIAL E MÉTODOS 


A pesquisa foi realizada com base nos dados pluviográficos da estação meteorológica 
denominada de Cidade de Goiás, pertencentes ao Instituto Nacional de Meteorologia 
(INMET/10º DISME) sob o código 83374, onde coletou-se dados de 10 anos de registros 
referente ao período de 2002 a 2011. Os dados obtidos, em formato de pluviogramas, 
com intervalo de registro de no máximo 10 mm de precipitação e de 24 horas de duração, 
com unidades de 0,2 mm para o volume e 10 minutos para o tempo, foram digitalizados 
com uso de scanner portátil e, posteriormente, convertidos para formato digital, utilizando 
planilha eletrônica. Para tanto, foram anotados os horários (horas e minutos) de início e fim, 
seguidos da lâmina correspondente de cada segmento de chuva de mesma intensidade e 
a precipitação total do pluviograma, possibilitando assim, possíveis correções posteriores. 

Em sequência, os dados foram organizados por meio da individualização das chuvas 
erosivas, arquivados em DAT e submetidos ao software CHUVEROS (software elaborado 
pelo prof. Elemar Antonino Cassol). Dos resultados de saída do software, utilizou-se os 
índices de erosividade El, e KE>25, altura precipitada; e a classificação das mesmas de 
acordo com o pico de intensidade, em padrões avançado, intermediário ou atrasado. No 
software, o índice de erosividade El,, é obtido por meio do produto da energia cinética (EC) 
e a intensidade máxima da precipitação em trinta minutos (l,,). O índice de erosividade 
KE>25 é calculado pelo produto da energia cinética pela intensidade de segmentos de 
precipitação maior que 25 min h'. A determinação da energia cinética (EC) utilizada no 


Chuveros é a mesma do El.., ou seja, Wischmeier e Smith (1958) e modificada por Foster 


ai 
et al. (1981) descrita por (Equação 1): 

EC = 0,119 + 0,0873 log | (1) 

Onde: 

| = intensidade da chuva em mm h'! 

Com o objetivo de obter o fator R da USLE, o qual requer série de dados pluviográficos 
de 22 anos no mínimo (WISCHMEIER e SMITH, 1978), foi feita a ampliação dos dados 
utilizando série pluviométrica da mesma estação pluviográfica, obtidas da base de dados 
Hidroweb/ANA. Com isso, foi feito o ajuste de equações utilizando índices de erosividade 
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das chuvas (El, e KE>25) obtidos de dados pluviográficos e características pluviométricas 
com base no Rc (coeficiente de chuva) para igual período de avaliação. O Rc foi calculado 


pela Equação 2, proposta por Fournier (1960) e modificada por Lombardi Neto (1977), 


conforme a seguir: 
Pm? 


Rc =— (2) 


Rc: coeficiente da chuva (mm); 

Pm: precipitação média mensal (mm); 

Pa: precipitação média anual (mm). 

Os ajustes mais adequados foram encontrados com base na precipitação média 
mensal (p), com modelos lineares para os dois índices de erosividade, conforme Figura 1. 








A 
1400 º » B 
= 1200 | Elto= 20,377 Ro + 103,28 KE>25 = 0,3109 Re + 1,7551 
t R2=0,876 o N ad 
é 16 0867. + 
q 1000 pa 
E 800 4 
É 

: 2 
E co A 
3 400 N 
E N 
m 200 x 4 

0 0 

0 20 40 60 0 20 40 60 


Coeficiente de chuva (Rc) Coeficiente de chuva (Rc) 


Figura 1. Correlação entre os valores mensais médios dos índices de erosividade El, (A) e 
KE>25 (B) e o coeficiente de chuva. 


Visando garantir a validade das estimativas dos índices, foi aplicada a análise de 
intervalos de confiança buscando analisar o maior número de anos contínuos da série 
pluviométrica de 1949 a 2011 que fosse semelhante à série pluviográfica de 2002 a 2011, 
para emprego nos modelos ajustados para estimativa dos índices de erosividade. Houve 
sobreposição das médias e/ou dos IC entre dados pluviográficos e pluviométricos para 
cada mês do ano, ao testar inicialmente toda a série pluviométrica, indicando assim, a 


adequação do uso do período de 1949 a 2011 para as estimativas de El,, e KE>25 (Figura 
2). 
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Figura 2. Médias e intervalo de confiança entre séries pluviométricas e pluviográficas da 
estação Cidade de Goiás — GO. 


Os índices de erosividade, mês a mês, dentro de cada ano, foram estimados 
conforme Moretti et al. (20093), utilizando planilha com dados de Rc mensais (calculados 
para a precipitação mensal e anual de cada ano com a equação 2) para o período de anos 
definidos na etapa anterior (1949 a 1999), e aplicando às equações dos ajustes obtidos 
para cada índice (Figura 1) seguido de regra de três simples usando a equação 3: 

Ei= (RCij.IEm) (3) 

RCm 

em que, 

Eij = índice de erosividade de um mês i contido num ano j da série de dados (MJ 
mm ha” h' ano"?); 

RCij = Rc total em um mês i contido num ano j da série de dados (mm); 

IEm = índice de erosividade (El, ou KE>25) médio mensal, no qual o mês i está 
contido na série de dados; e 

RCm = Rc total mensal médio no qual o mês i está contido na série de dados (mm). 

Os índices de erosividade mensais, calculados por meio da soma do El,, de cada 
chuva erosiva ocorrida no mês, foram utilizados para determinar a erosividade anual (El, e 
KE>25) para cada ano da série de 1949 a 2011 da localidade em estudo. 

Com base no padrão hidrológico de cada precipitação, classificado pelo software 
Chuveros, foi calculada frequência das chuvas em cada padrão hidrológico, bem como, a 
distribuição relativa de altura precipitada e erosividade medida pelos índices El, e KE>25 
em cada padrão hidrológico. 

O cálculo da chuva máxima diária anual foi realizado a partir dos dados da 
série histórica, obtidos da mesma estação meteorológica (HIDROWEB/ANA), a qual foi 
organizada em planilhas eletrônicas. Ajustou-se os dados ao modelo de distribuição de 
probabilidades de Gumbel (Equações 4 a 7), utilizando a chuva máxima diária anual 
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média e seu desvio padrão (PINTO et al., 1996). Assim, obteve-se a precipitação da chuva 
relacionada a diferentes períodos de retorno (2, 5, 10, 25, 50 e 100 anos). 


r=[-n(- 5) (4 





— y=m 
K = sn 6) 
P=X+KS (6) 
o 1 
T= Ts (7) 
Onde: 


Y = variável reduzida da distribuição de Gumbel; 

In = logaritmo neperiano; 

T = período de retorno (anos); 

yn e Sn = valores tabelados em função dos números de anos (n) da série de dados; 
P = valor da variável (precipitação) associado à frequência T; 

X = média aritmética da amostra; 

K= coeficiente de frequência; 

S= desvio padrão da amostra. 


31 RESULTADOS E DISCUSSÃO 


A precipitação média anual com base na série pluviométrica (1949 a 2011) foi de 
1.648,6 mm. As chuvas e, consequentemente, o potencial erosivo concentrou-se entre 
os meses de outubro e março no período estudado, sendo considerado o “período das 
águas” (Figura 3). O período de maio a setembro com baixos totais mensais precipitados, 
consequentemente apresenta potencial erosivo pouco relevante. Desta forma, o período 
de outubro a março reúne 88,3% das precipitações e 87,5% e 86,87% da erosividade pelos 
índices EI30 e KE>25, respectivamente. Observa-se que os índices El,, e KE>25 possuem 
comportamento semelhante ao longo do ano (Figura 3). 

Resultado semelhante foi encontrado por Galdino (2015) para o estado de Goiás 
e Distrito Federal, onde o período de maior incidência de chuvas ocorreu de outubro a 
março, totalizando 87% do índice total anual pluviométrico na região e 94% do total anual 
da erosividade da chuva. De acordo com Cassol et al. (2008), a tendência é que os maiores 
valores do índice El,, coincidam com as maiores precipitações, sendo um período de alto 
risco à perdas de solo por erosão, se o solo estiver menos protegido. Isso demonstra a 
importância do estudo da distribuição do potencial erosivo da chuva (MORETTI et al., 2003) 
com relação ao manejo do solo e das culturas, visando garantir a cobertura do solo para 


diminuir o escoamento superficial. 
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Figura 3. Valores médios mensais percentuais dos fatores de erosividade El, e KE>25 e do 
volume total das chuvas obtido entre 1949 e 2011 (exceto os anos 1988, 2000 e 2001), para a 
cidade de Goiás (GO). 


A erosividade anual para a cidade de Goiás (GO) variou entre 3.872,4 MJ mm ha” 
h' em 1963 a 11.535,9 MJ mm ha” h'! em 1985 (Tabela 1). Mesmo no período chuvoso com 
que apresentou menor variabilidade (menores desvios padrão e coeficiente de variação) 
que no da seca, observa-se que essa foi muito alta, sendo, no entanto considerado normal 
já que cada evento de chuva na natura é considerado único, especialmente nas regiões 
tropicais. Umas das causas de alta variabilidade espacial e temporal das chuvas no Brasil 
está relacionado aos diferentes tipos de chuva (FISCH et al., 2007), como convectivas, 
frontais e orográficas. 

Meses de maior erosividade observados podem ter relação direta com cheias 
históricas reportadas na cidade de Goiás (GO). De acordo com trabalho de Cavalcanti et al. 
(2008), há documentadas total de nove cheias responsáveis por grandes danos à histórica 
cidade de Goiás (GO). De acordo com esses autores as causas das enchentes incluem o 
fator uso da terra nas cabeceiras. Das poucas cheias que tiveram as quantidades de chuvas 
registradas oficialmente nessa cidade, verifica-se que três delas datam 19 de fevereiro de 
1980 (121,8 mm), 29 de março de 1981 (150,6 mm) e 31 de dezembro de 2001 (140,9 
mm), destacando que já havia elevada precipitação acumulada na semana antecedente às 
cheias. Essas cheias citadas coincidem com os máximos índices de erosividade mensais 
da série estudada para o mês de fevereiro de 1980 (4673,8) e março de 1981 (3815,2 MJ 
mm ha” h'). Para o mês de dezembro, provavelmente, o máximo valor de erosividade da 
série pode ter coincidido com a cheia de 2001, mas devido à falha na série, isso não pôde 
ser confirmado. 


A erosividade anual média (fator R) obtida de dados pluviométricos para a localidade 
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estudada foi de 7.1 MJ mm ha” h?. Esse índice tende em ser inferior ao obtido em outras 
localidades como Brasília (DF) (8.319) (DEDECEK, 1988), Goiânia (GO) (8.353) (SILVA et 
al. (1997)), Flechas (MT) (7.830) (MORAIS et al., 1991) e Cuiabá (MT) (8.810 MJ mm ha” 
h' ano”) por Almeida et al. (2011), todas no mesmo tipo de clima (Aw). 








o Jan Fev Mar Abr Mai Jun Jul Ago Set Out No Dez a 
1949 28619 23987 1225 47,1 74,2 0,0 0,0 0,0 0,0 2584 2667 6322  6661,7 
1950 3123 7,8 2834 1649 30,6 0,0 0,0 0,0 0,0 965,2 20428 1931,7 5738, 
1951 16348 4943 27716 4,2 444,8 0,0 0,0 0,0 0,0 1020 4125 747,0 661,1 
1952 481,0 980,6 11647 3840 0,0 0,0 0,0 0,0 0,7 45,2 4888 24180 59630 
1953 143,1 517,3 6892 4131 0,0 0,0 136,4 0,0 237,7 8825 435,7 23548 58098 
1954 459,5 830,1 407,9 1933 354 0,0 0,0 0,0 30,9 38 13727 17809 51145 
1955 36139 247,6 3145 2731 25 41,5 0,0 0,0 0,0 469,4 2288 19762 71674 
1956 1369 10675 442,2 12791 3322 14209 7,9 225 40,0 130,0 10998 14682 7447,1 
1957 21327 15162 8563 2240 68,7 0,0 0,0 38,2 440,3 1454 8084 327,6  6557,8 
1958 35566 973,7 317,0 631,9 4226 0,0 8730 0,0 35,4 173,2 66,7 2873,7 99238 
1959 1801,7 351,5 21064 67,1 20,2 139,5 0,0 211,6 0,2 523,8 787,2 797,7 68067 
1960 16676 28468 10555 1065 21,0 0,0 0,0 0,0 0,0 373,2 13122 6646 80473 
1961 39686 443,5 2346 101,1 15,2 0,0 0,0 7,0 0,0 281,2 8920 20712 80142 
1962 35536 260,0 7274 14,2 0,0 0,0 0,0 9,9 194,1 0,0 453,1 37590 8971,2 
1963 7164 14379 495,1 3764 0,0 0,0 0,0 0,0 5,6 0,0 400,9 440,1 38724 
1964 69952 920,5 121,1 59,9 7757 0,0 0,0 0,0 49,2 6722 11937 4437  11231,2 
1965 17106 834,2 11395 8036 0,0 0,0 0,0 0,0 17,2 1494,3 19951 197,1 81915 
1966 24916 29552 13136 1,7 78,2 0,0 0,0 0,0 0,1 174,5 131,2 19886 9134,7 
1967 831,0 7527 18431 1914 00 0,0 0,0 0,0 57,7 155,4 7566 10885 56764 
1968 450,0 23641 7233 1951 1,6 0,0 0,0 37,1 125,9 80,9 14376 960,4  6976,1 
1969 1530,7 0,0 0,0 0,0 1587,6 0,0 0,0 0,0 0,0 340,5 30268 14873 7922,9 
1970 12832 451,6 7020 339,1 15,1 0,0 0,0 2,6 2840 1995 14922 9253 56947 
1971 298,1 494,9 3158 1495 350 0,0 0,0 0,0 5,8 1141,7 1921,2 11658 5528,0 
1972 229,7 11608 2900 70,7 68,7 0,0 177,6 4087 232 311,6 7204 27830 62444 
1973 480,6  1481,7 15285 2104 146 0,0 0,0 0,0 16,1 350,1 7464 461,2 52896 
1974 18179 314,6 17568 9095 0,8 0,0 0,0 3,7 4889 0,0 0,0 0,0 5292,1 
1975 457,3 977,0 7485 10486 748 12,7 0,0 0,0 30,3 195,5 961,1 3637 48695 
1976 299,5 1370,7 7659 14,2 654,1 0,0 0,0 0,1 6788 9763 14992 7792  7087,9 
1977 43838 335,7 210,1 141,0 7541 139 0,0 0,0 576,0 175,1 5831 4845 76578 
1978 12239 634,8 6898 3992 24,1 0,0 12093 0,0 20,0 330,2 15362 18657 79334 
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1979  3791,9 22598 366,6 1,5 1289 0,0 62,3 911,5 142,9 748 310,3 1294,1 9344,6 
1980 16828 46738 179,0 1398 2,1 20,3 0,0 0,1 91,2 32,0 693,5 961,1 8475,5 
1981 10730 78,4 3815,2 167,4 453,6 260,7 0,0 0,0 24,1 1159,8 11342 433,4 8599,8 
1982 34734 660,7 12866 79,7 254,6 0,0 1,7 412,5 462,1 442,5 343,0 777, 8194,4 
1983 3577,0 10714 389,0 350,6 125,8 1,5 562,8 0,0 45,1 359,8 630,6 177,4 8291,0 
1984 4684 289,0 14769 234,2 86,3 0,0 0,0 287,5 611,5 1627,9 582,5 605,9 6270,0 
1985 79681 5546 521,8 176,8 27,9 0,0 8,1 11,9 330,0 199,7 216,6 1520,3 11535,9 
1986  1142,6 445,2 642,3 55,0 43,9 0,0 40,1 1473,1 149,4 397,3 415,2 18580 6662,1 
1987 4829 6985 0,0 575,1 35,6 0,0 0,0 0,6 86,2 318,1 14892 3250,6 6936,7 
1989 902,2 5738 724,8 45,7 0,0 196,2 2303,1 1385 6,6 835,0 505,1 5247,6  11478,6 
1990 1237,7 6847 3127 10026 2120 0,0 18,6 103,8 245,4 4180 388,6 1016,8  5640,7 
1991 929,2 681,0 27703 7330 0,0 0,0 0,0 0,0 96,2 695,0 0,0 1503,1 7407,8 
1992  2248,7 1967,8 0,0 10069 167,3 0,0 0,0 0,3 82,8 1030,0 1131,2 12326  8867,6 
1993 0,0 2229,8 0,0 1141,0 189,6 0,0 0,0 0,3 93,8 167,2 1281,8 1396,8  7500,3 
1994 14859 26238 152,7 127,8 0,0 13,6 0,0 2204,9 0,0 0,0 0,0 0,0 6608,7 
1995 0,0 0,0 0,0 230,3 19,8 1092,2 554,0 0,0 130,8 2294  3841,3 15269 7624,7 
1996 13014 1730,3 459,3 2258 44,3 0,1 0,0 0,0 146,2 132,3 332,1 3734,7 8106,5 
1997 348,1 1075,6 4934 481,2 312 0,7 2,5 18,4 56,7 136,7 22624 30269 7905,9 
1998 1531,8 947,1 10863 276,9 549,7  2931,)0 0,0 0,0 309,9 97,5 152,0 997,6 8879,6 
1999 802,5  1042,3 180,3 6060 33,4 0,4 6,0 10,6 75,4 183,1 1089,3 11936  5223,0 
2002 20232 8539 340 43,2 8,3 0,0 0,3 0,0 92,5 46,8 718,4 734,5 4555,1 
2003 27170 431,6 980,3 3393 0,7 0,0 0,0 4,9 30,0 139,5 349,4 800,2 5792,8 
2004 526,1 2447,8 611,2 29,6 41,3 0,0 167,2 0,0 0,0 423,5 10428 539,7 5829,1 
2005 15835 5334 6840 33,8 29,4 3,2 0,0 8,0 166,1 119,7 948,0 1989,3  6098,4 
2006 3826 659,2  1075,6 266,8 0,0 0,0 0,0 0,3 200,7 18132 639,4 741,1 5778,9 
2007 107,7 38908 5769 331,6 0,1 0,5 A 0,0 0,0 120,6 130,5 806,2 6036,0 
2008 1224,8 11887 1101,7 427,6 80 0,0 0,0 0,1 87,4 77,8 647,1 1925,6 6688,8 
2009 676,1 796,2 17811 476,6 17,3 34,5 0,0 145,8 429,7 6835 521,3 1904,9 74671 
2010 371,3 0,0 245,4 00 0,0 6,3 28,0 0,0 0,0 520,3 17487 15533 44734 
2011 2921,3 137,0 18830 77,3 0,0 308,5 0,0 0,0 2,9 711,5 664,5 967,5 7673,6 
Média 1641,8 1077,5 7994 308,3 133,1 108,3 103,8 107,9 125,9 419,1 888,0 13992 7112,2 
D.P. 1601,6 962,98 762,8 3188 265,24 438,1 358,8 360,88 171,37 4297 734,7 1008,3 1666,1 
CM. 97,6 89,4 95,4 103,4 199,3 404,5 345,6 334,4 136,1 102,5 82,7 72,1 23,4 





Tabela 1. Valores mensais, anuais e médias de erosividade El, (MJ mm ha” h! ano”) com 
estatísticas descritivas, obtidos de séries pluviométrica e pluviográfica para o período de 1949 a 


2011 (exceto 1988, 2000 e 2001). 
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Os meses de dezembro e janeiro destacaram-se na série histórica com os maiores 
valores de erosividade (El,.). Segundo Mazurana et al. (2009), meses de pico excepcional 
na erosividade das chuvas são responsáveis por grandes danos às lavouras cultivadas, 
em decorrência do processo de erosão hídrica. Um atenuante nesse período é que 
normalmente nos meses de dezembro e janeiro a maioria das culturas agrícolas em Goiás 
já está bem desenvolvida, cobrindo em parte ou totalmente a superfície do solo pelo seu 
dossel, e nesse caso, o fator cobertura do solo passa a ser fundamental para diminuir os 
danos do maior potencial erosivo das chuvas nesse período. 

Conforme a classificação de Carvalho (1994) a erosividade para a cidade de Goiás 
é considerada de média a forte (4.905 < R < 7.357), o que indica a necessidade do uso de 
práticas conservacionistas nas áreas agrícolas. 

Quanto ao índice KE>25 (Tabela 2), seus valores variaram entre 60 MJ em 1962 a 
186,3 MJ ha” ano”, em 1990. A média anual para esse índice foi de 110,3 MJ ha”! ano”, 
sendo os maiores valores coincidentes com os maiores valores do índice El,, de outubro a 
março, indicando o período de maior potencial erosivo.$ 

Valores de KE>25 menores que o desse estudo foram encontrados por Carvalho 
et al. (2010) de 99,2 MJ h'! ano! para o município de Seropédica (RJ). Valle Junior et al. 
(2012), encontraram 102,05 MJ h! ano”! para a região de Uberaba (MG). Já em estudo de 
Silva et al. (1997) para região de Goiânia (GO) foi encontrado 129,8 MJ h ano”!, ou seja, 
superior ao presente estudo. 

Em relação à distribuição das chuvas quanto aos padrões hidrológicos (Tabela 
3), constata-se que o padrão de maior ocorrência na cidade de Goiás - GO é o 
avançado (49,8%), seguido dos padrões atrasado (35,1%) e intermediário (15,1%). 
Consequentemente, as perdas de solo esperadas podem ser menores que poderiam ser, 
pois as chuvas não predominam no padrão atrasado. Nesse padrão apresentam pico de 
intensidade no terço inicial de sua duração em momento que solo pode possuir menos 
umidade quando comparado aos demais padrões, ocasionando menor desagregação, 
selamento e transporte de solo (MEHL et al., 2001). Em São Borja, RS, o padrão hidrológico 
predominante foi o avançado em trabalhos realizados por Mehl et al. (2001) e Cassol et al. 
(2008), assim como na região de Ribeirão das Lajes, RJ (MACHADO et al., 2008). 
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Es Jan Fev Mar Abr Mai Jun Jul Ago Set Out Nov Dez Bia 
1949 44 37 2 1 1 0 0 0 0 4 4 10 102,8 
1950 5 0 4 3 1 0 0 0 0 15 32 30 89,1 
1951 25 8 43 0 7 0 0 0 0 2 6 11 102,7 
1952 7 15 18 6 0 0 0 0 0 1 8 37 92,2 
1953 2 8 11 7 0 0 2 0 4 14 7 36 90,7 
1954 7 13 6 3 1 0 0 0 1 0 21 27 79,1 
1955 56 4 5 4 0 1 0 0 0 7 4 30 110,7 
1956 2 16 7 21 6 28 0 0 1 2 17 23 117,0 
1957 33 23 13 4 1 0 0 1 7 2 13 5 102,1 
1958 55 15 5 10 7 0 15 0 1 3 1 44 155,3 
1959 28 5 33 1 0 2 0 4 0 8 12 12 105,7 
1960 26 44 16 2 0 0 0 0 0 6 20 10 124,3 
1961 61 7 4 2 0 0 0 0 0 4 14 32 123,6 
1962 55 4 11 0 0 0 0 0 3 0 7 58 138,4 
1963 11 22 8 6 0 0 0 0 0 0 6 7 60,0 
1964 107 14 2 1 13 0 0 0 1 11 18 7 174,3 
1965 26 13 18 13 0 0 0 0 0 24 31 3 127,5 
1966 38 46 20 0 1 0 0 0 0 3 2 31 140,9 
1967 13 12 29 3 0 0 0 0 1 2 12 17 87,9 
1968 7 36 11 3 0 0 0 1 2 1 22 15 98,8 
1969 24 0 0 0 26 0 0 0 0 5 47 23 124,6 
1970 20 7 11 5 0 0 0 0 5 3 23 14 88,6 
1971. 5 8 5 2 1 0 0 0 0 18 30 18 85,9 
19722 4 18 4 1 1 0 3 7 0 5 11 43 97,4 
1973 7 23 24 3 0 0 0 0 0 6 12 7 82,0 
1974 28 5 27 15 0 0 0 0 8 0 0 0 82,9 
1975 7 15 12 17 1 0 0 0 1 3 15 6 76,1 
1976 5 21 12 0 11 0 0 0 11 15 23 12 110,9 
1977 67 5 3 2 13 0 0 0 10 3 9 Fá 120,0 
1978 19 10 11 6 0 0 21 0 0 5 24 29 124,7 
1979 58 35 6 0 2 0 1 15 2 1 5 20 145,8 
1980 26 72 3 2 0 0 0 0 2 1 11 15 130,8 
1981 16 il 59 3 8 4 0 0 0 18 18 7 134,2 
1982 53 10 20 1 4 0 0 7 8 7 5 12 128,1 
19e3 55 17 6 6 2 0 10 0 1 6 10 18 129,2 
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1984 7 4 23 4 1 0 0 5 10 26 9 9 99,0 
1985 122 9 8 3 0 0 0 0 6 3 3 23 178,2 
1986 18 7 10 1 1 0 1 25 3 6 6 29 105,4 
1987 7 q 0 9 1 0 0 0 1 5 23 50 107,6 
1989 14 40 5 7 0 7 0 0 0 8 10 15 106,7 
1990 19 9 11 1 0 3 39 2 0 13 8 81 186,3 
1991 14 E 5 16 4 0 0 2 4 E 6 16 83,8 
1992 35 9 37 10 0 0 0 0 1 9 0 20 120,9 
1993 O 34 0 18 3 0 0 0 2 18 20 22 117,3 
1994 23 40 2 2 0 0 0 37 0 0 0 0 105,2 
1995 O 0 0 4 0 17 9 0 2 4 59 24 119,6 
1996 20 27 7 4 1 0 0 0 2 2 5 57 125,3 
1997 5 Ed 8 8 0 0 0 0 1 2 35 47 122,4 
1998 24 15 Fá 4 9 47 0 0 5 2 2 15 139,7 
1999 12 16 3 10 1 0 0 0 1 3 17 18 81,2 
2002 31 13 1 1 0 0 0 0 2 1 11 11 70,3 
2003 42 7 15 5 0 0 0 0 1 2 5 12 89,6 
2004 8 38 9 0 1 0 3 0 0 7 16 8 90,4 
2005 24 8 1 1 0 0 0 0 3 2 15 31 94,4 
2006 6 10 irá 4 0 0 0 0 3 29 10 11 90,3 
2007 2 60 9 5 0 0 1 0 0 2 2 12 93,4 
2008 19 18 TF 7 0 0 0 0 1 1 10 30 103,6 
2009 10 12 28 8 0 1 0 2 x 11 8 29 116,6 
2010 6 0 4 0 0 0 0 0 0 8 27 24 69,3 
2011 45 2 29 1 0 5 0 0 0 11 10 15 118,8 
Média 25,2 16,7 12,4 4,8 22 1,8 1,8 1,8 2,1 6,5 13,6 21,4 110,3 
D.P. 24,6 15,1 11,5 4,9 4,5 7,0 6,1 6,1 2,9 6,7 11,4 15,5 26,1 
C.v. 97,6 89,9 93,0 102,3 2033 380,3 345,6 334,4 1385 103,1 83,5 72,5 23,7 





Tabela 2. Valores mensais, anuais e médias de erosividade KE>25 (MJ ha”! ano!) com 
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a ” Número de chuvas pr disendo dá El, KE>25 
Absoluto % mm % Absoluto* % Absoluto* % 
Avançado 726 49,8 653,8 48,9 2551,7 49,6 33,9 43,52 
Intermediário 220 15,1 304,4 22,8 1141,4 22,2 19,6 25,16 
Atrasado 512 35,1 378,1 28,3 1453,7 28,2 24,4 31,32 
TOTAL 1458 100,0 1336,30 100,0 5146,8 100,0 TT, 100,0 





* MJ mm ha” h'! ano"!; ** MJ h'! ano”. 


Tabela 3. Número de chuvas, altura precipitada e erosividade acumuladas nos padrões de 
chuva avançado, intermediário e atrasado na cidade de Goiás (GO), obtidos de pluviógrafo no 
período de 2002 a 2011. 


Pode-se observar ainda na Tabela 3, que o maior volume precipitado ocorreu no 
padrão avançado, além dos índices de erosividade (El, e KE>25) apresentarem o mesmo 
comportamento. Eltz et al. (2011), ao trabalhar com um simulador de chuva em Argissolo 
Vermelho-Amarelo, observaram maiores perdas de solo em chuvas de padrão atrasado 
com relação às perdas de solo em chuvas de padrão intermediário e avançado, o que foi 
atribuído a maior umidade antecedente. 

A respeito da chuva máxima diária anual (CMDA), valores médios e associados a 
diferentes períodos de retorno estão apresentados na Tabela 4. Observa-se que os valores 
encontrados ficaram abaixo dos obtidos por Mesquita et al. (2009) para algumas cidades 
vizinhas à Goiás (GO), para os períodos de retorno coincidentes. Já em trabalho realizado 
por Pereira et al. (2017) no sul do estado (Ipameri - GO) para os mesmos períodos de 
retorno, os valores de CMDA da cidade de Goiás (GO) tenderam em ser maiores. As 
diferenças entre as localidades podem estar relacionadas ao efeito do relevo sobre as 
chuvas, tamanho e período das séries de dados. 

O presente trabalho permite definir o valor da chuva de projeto para diferentes 
exigências de período de retorno de obras hidráulicas no meio rural (Tabela 4). No 
dimensionamento de terraços agrícolas em nível, por exemplo, é recomendado usar CMDA 
com período de retorno de 5 a 10 anos (CRUCIANI, 1988). Pode-se calcular também, com 
base na série histórica, que período de retorno determinada chuva intensa está associado 
(Equação 7). No caso da cheia ocorrida em 2001 com valor precipitado de 140,9 mm 
(CAVALCANTI et al., 2008) por exemplo, tal evento pode ser igualado ou ultrapassado pelo 
menos uma vez em média a cada 14,7 anos. Deve se ressaltar que no referido evento a 
cheia ocorreu quando nos seis últimos dias já havia chovido 124,1 mm segundo os autores 
supracitados, ou seja, o solo já se encontrava com muita umidade, e com isso, baixa 
capacidade de armazenamento. 

Acredita-se que a redução de cheias na bacia do rio Vermelho no tocante ao uso 
da terra, está no aumento da capacidade de infiltração de água no solo, bem como, no 
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armazenamento máximo possível do escoamento superficial (enxurrada) nas propriedades 
rurais. Existem diferentes práticas conservacionistas que podem ser empregadas para isso 
podendo ser citados alguns exemplos como, melhoria da cobertura do solo (pastagens), 
plantio direto, reflorestamento, descompactação do solo, plantio em nível, terraceamento 
e construção de bacias de infiltração dos quais esses últimos podem utilizar informações 


para dimensionamentos com base neste trabalho. 








Localidade Média 2 5 10 25 50 100 
Goiás 89,3 93,9 112,6 131,0 1542 171,5 188,6 
Itaberaí" - - 135,6 1640 202,5 - - 
Itapuranga! - - 122,9 145,1 174,5 - - 
Ipameri? = 80,96 97,07 107,74 121,23 131,23 141,15 





*Mesquita et al. (2009). 2Pereira et al. (2017) 


Tabela 4. Chuva máxima diária anual (mm) média e associada a diferentes períodos de retorno 
(anos) e localidades. 


41 CONCLUSÃO 


O fator R da USLE para a cidade de Goiás (GO) nas condições do estudo é de 
7.112,2 MJ mm ha” h' ano”, servindo de base para estudos de previsão de perdas de solo. 
O valor anual médio pelo índice KE>25 é de 110,3 MJ h'! ano”. Dezembro e janeiro são os 
meses de maior erosividade. 

O padrão hidrológico avançado concentra os maiores valores de precipitação e 
índices de erosividade. 

Chuvas máxima diária anual para a cidade de Goiás, com recorrência de 2, 5, 10, 
20, 50 e 100 anos, são respectivamente, de 93,9; 112,6; 131,0; 154,2; 171,5 e 188,6 mm. 
A chuva de projeto com duração de 24 h para determinação de terraços em nível na região 
sob influência da estação Cidade de Goiás é de 131,0 mm nas condições do estudo. 

Chuvas máxima diária anual para a cidade de Goiás, com recorrência de 2, 5, 10, 
20, 50 e 100 anos, respectivamente, de 93,9; 112,6; 131,0; 154,2; 171,5 e 188,6 mm. A 
chuva de projeto com duração de 24 h para determinação de terraços em nível na região 
sob influência da estação Cidade de Goiás é de 131,0 mm. 
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RESUMO: Neste trabalho foram calculadas 
doses de fertirrigação, com água residuária da 
bovinocultura (ARB) anaerobiamente digerida 
em reator UASB, para 12 culturas agrícolas 
quando cultivadas em dois diferentes solos: 
Argissolo e Planossolo. Como resultado temos 
os valores de recomendação de lâminas de 
fertirrigação com ARB, equivalentes à aplicação 
de fontes convencionais de nitrogênio, como 
a ureia. Concluiu-se, que a água residuária 
de ovinocultura possui potencial para a 
complementação ou substituição da adubação 
nitrogenada convencional, sendo uma alternativa 
econômica, circular e sustentável. 
PALAVRAS-CHAVE: Economia circular, reuso, 
compostos nitrogenados. 
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UASB 


RECOMENDATION FOR FERTIGATION 
WITH CATTLE WASTEWATER 


ABSTRACT: In this work, fertigation rates were 
calculated, with anaerobically digested cattle 
wastewater (CWW) from UASB reactor, for 12 
agricultural crops when grown in two different 
soils: Argisolo and Planosol. As a result, the 
recommendation values for fertigation blades 
with CWW are suggested, equivalent to the 
application of conventional nitrogen sources, like 
urea. The conclusions were that the anaerobically 
digested CWW has the potential to complement 
or replace conventional nitrogen fertilization, 
being an economical, circular and sustainable 
alternative. 

KEYWORDS: Circular economy, reuse, nitrogen 
compounds. 


INTRODUÇÃO 


O tratamento e descarte de águas 
residuárias da ovinocultura (ARB) requer 
atenção devido a seu potencial poluidor. 
Uma criação intensiva de gado com 1.000 
cabeças possui potencial poluidor semelhante 
ao | aproximadamente 14.000 | pessoas 
(MENDONÇA et al., 2016; MENDONÇA et 
al., 2018). A disposição final incorreta destes 
resíduos pode causar contaminação de solo, 
das águas superficiais e subterrâneas (ERTHAL 
et al., 2010). 

Uma alternativa viável e econômica 
para a disposição final destes efluentes é a sua 
utilização para fins agrícola. Nesse sentido, a 
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fertirrigação de culturas com a ARB torna-se uma opção que contribui para a redução da 
poluição, tendo efeito positivo na melhoria dos atributos do solo, reduzindo os custos com 
a adubação convencional (MENDONÇA et al., 2016; MENDONÇA et al., 2018). 

No entanto, a aplicação de maneira inapropriada dos biofertilizantes também 
pode ocasionar efeitos deletérios aos solos, cursos d'água e à saúde pública (MACIEL 
et al., 2019). Dessa forma, o procedimento deve ser bem manejado, aplicando-se doses 
de maneira eficiente de forma a fornecer os nutrientes necessários para a planta e seu 
desenvolvimento equilibrado (BARBOSA, 2019). 

Diante disso, o presente trabalho teve por objetivo delinear a montagem de uma 
tabela com recomendação de adubação orgânica para 12 culturas agrícola, por meio do 
cálculo das doses de aplicação de ARB previamente tratadas em Reator UASB - Upflow 
Anaerobic Sludge Blanket Reacior, quando cultivadas em dois tipos de solos: Planossolo e 


Argissolo, coletados no município de Seropédica — RJ. 


MATERIAL E MÉTODOS 

Ambos os solos foram coletados na área experimental da UFRRJ. Foram realizadas 
2 campanhas de amostragem composta, oriunda de 8 pontos. As amostras de solos 
foram analisadas conforme Manual de Métodos de Análise de Solo (EMBRAPA, 1997) no 
Laboratório de Fertilidade do Solo, Instituto de Agronomia da UFRRJ (Tabela 1). 





Tipo de Profundidade Na Ca Mg K HH Al S T 











Solo (em) Cmol./dm?---------- 
iisissol 0:40 0,49 315 225 036 3,38 0,20 6,25 9,63 
(0,27) (0,15) (005) (0,14) (1,32) (0,00) (0,50) (1,82) 
Pissóssoio 0:40 0,04 1,00 0,55 013 264 0,25 1,72 436 
(0,01) (0.00) (005) (001) (096) (005) (003) (0,93) 
Tipo de Profundidade V m n  PHágu MO Pp K p 

Solo (em Yg-===- 1:25 % --mg/L-- g/em' 

Asbisdio 0:40 66,30 312 475 44 116 3131 140,44 133 

(31) (024) (188) (010) (011) (169) (5343) (0,02) 

4141 12,67 0,84 51 L03 57,57 51,23 1,58 

Planossolo 0-40 
(947) (241) (007) (000) (036) (1143) GM (0,08) 


Valores entre parênteses indicam o desvio padrão. 


TABELA 1. Resultados das análises químicas — Rotina de fertilidade do solo. 


Água residuária (biofertilizante) e tratamento prévio 


A ARB utilizada para o cálculo das doses de recomendação foi coletada após reator 
UASB com tempo de retenção hidráulica de 7,5 (+0,50) dias, na “Fazendinha Agroecológica” 
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(coordenadas: 22º45'22” S: 43º40'28” W), campus da UFRRJ. Anterior ao reator anaeróbio 
havia separação de sólidos grosseiros em esterqueira, que serve como um sedimentador 
primário. Os parâmetros físico-químicos foram determinados em triplicatas em conformidade 
com as metodologias do Standard Methods (APHA, 2012) e são apresentados na Tabela 2. 





pH CE  NH$ No NOy NOy P Mg” Ca” Na ST SD 


(UpH) (uS/cm) nana UA 
6,2 1,496 150 126 2 0,75 368 32 53 22 897 7155 
(03) 354) 69) 4) EC) O DD a O) a4)ia» 9 





Valores entre parênteses indicam o desvio padrão. 


TABELA 2. Caracterização físico-química da ARB após reator UASB 


Cálculo da taxa de aplicação (dose de referência) 


Foram calculadas conforme metodologia proposta por Matos (2017), representada 
pela Equação 1, com base na produtividade máxima de cada cultura. 


+(N + 


amônio 


DA,, = 1000 * [N,.- (TM1* MO * p,* p*107* 0,05* n/12)] /[TM2* (n/12) * N 
N 


org 
Hirata) E TR] (1) 

Em que: DA,,: taxa de aplicação ou dose, mº ha; N.,.: absorção de nitrogênio 
pela cultura, kg ha"; TM1: taxa anual de mineralização da MO no solo, kg.kg.ano!; MO: 
conteúdo de matéria orgânica do solo, kg kg"; p,: massa específica do solo, t m?*; p: 
profundidade do solo considerada, m; n/12: fração anual relativa ao período de cultivo, ano; 
TM2: taxa de mineralização do nitrogênio orgânico, kg kg ano”; Nor nitrogênio orgânico, 
mg Lº; Nano: Nitrogênio amoniacal, mg L*; No: Nitrogênio nítrico, mg L”, e TR: taxa de 
recuperação do nitrogênio mineral pela cultura, kg kg'ano”. 

As taxas anuais de mineralização da matéria orgânica (TM1), nitrogênio orgânico 
(TM2), utilizados na Equação (1), foram 0,01 e 0,5 respectivamente (MATOS, 2017), 
enquanto a taxa de recuperação de nitrogênio mineral extraído pela cultura (TR) variou 
conforme a capacidade de cada vegetação. 


RESULTADOS E DISCUSSÃO 


Na Tabela 3 estão apresentadas as lâminas de fertirrigação recomendadas para 
cada cultura quando cultivadas nos solos aqui analisados. 
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Remoção Ciclo da Dose 








Cultura ne o cultura (kg no Ei (mm) 
é her) tmesas) Argissolo - Planossolo 

Milho 14 200 4 0,5 185,99 185,44 
Arroz 16 141 3,3 0,3 201,01 200,30 
Trigo 10 160 4 0,5 146,78 146,23 
Feijão 1 31 3 0,5 24,07 23,63 
Soja 2 250 2 0,5 267,60 267,30 
Batata 62 147 3 0,5 143,96 143,53 
Cana-de-açúcar 300 254 12 0,5 154,96 153,79 
Algodão 2 110 4 0,5 97,76 97,21 
Capim elefante 46 800 3 0,7 613,46 613,13 
Capim-guiné 35 560 4 0,7 409,01 408,59 
Capim-pangola 31 400 4 0,7 289,97 289,55 
Sorgo 16 200 4 0,7 141,16 140,74 





TR = taxa de recuperação do nitrogênio mineral pela cultura, adotados como 0,5 para culturas 
anuais, 0,3 para arroz inundado e 0,7 para culturas de sistema radicular volumoso (MATOS, 
2017). Os rendimentos dos cereais estão expressos em base de massa seca. 


TABELA 8. Doses de ARB para o fornecimento exigido de N para cada cultura. 


Tendo em vista as diferentes necessidades de nitrogênio de cada cultura, verificou- 
se as doses de água residuária necessárias para atender tal demanda com vistas à 
produtividade máxima. Dessa forma, foram obtidos os valores da lâmina (mm) de ARB a 
serem aplicadas no Planossolo e Argissolo quando cultivados com as culturas descritas na 
tabela acima. 

As menores lâminas de ARB obtidas para ambos os solos, foram para o feijão. 
Como trata-se de uma cultura fixadora de nitrogênio, podemos observar que, dentre as 
culturas estudadas, esta possui uma menor taxa de remoção de nitrogênio e consequente 
menor valor de dose. 

O capim elefante por sua vez, obteve os maiores valores de dose, pois sua taxa de 
absorção de nitrogênio é maior do que a das outras culturas, recebendo assim uma maior 
lâmina de irrigação, podendo ser utilizado como uma cultura para a produção de biomassa 
e extração do nitrogênio do solo proveniente da água residuária. 

O uso do biofertilizante proveniente de ARB pode apresentar resultados satisfatórios, 
como demonstrado no estudo realizado por MENDONÇA et al (2016) e MENDONÇA et al 
(2018) no cultivo de algumas cultivares de cana-de-açúcar. O biofertilizante produzido por 
meio da ARB foi utilizado como fonte de nitrogênio, substituindo totalmente a adubação 
convencional com ureia, sem prejudicar a produtividade e qualidade das cultivares. 

Outras vantagens do uso da ARB como biofertilizante na agricultura é a economia 


Engenharia Agronômica: Ambientes Agrícolas e seus Campos de Atuação 3 Capítulo 3 


com fontes convencionais de nitrogênio (ureia, sulfato de amônio etc) além de evitar que 


estes efluentes cheguem aos cursos de água. 


CONCLUSÕES 


A ARB anaerobiamente digerida por reator UASB pode ser considerada como 
alternativa econômica para substituição de fertilizantes convencionais. A menor lâmina 
calculada foi para o feijão e a maior para o capim elefante, demostrando que este último 
pode ser utilizado para extração de N da solução do solo de forma mais eficiente. 
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ABSTRACT: The objective of this study was 
to evaluate the effects of gypsum doses on 
physical-hydric attributes, root growth and 
soybean productivity. The experiment was 
conducted at the Agronomic Institute of Paraná 
(IAPAR) in Santa Tereza do Oeste-PR. The soil 
was classified as Typic Haplortox. Five doses of 
agricultural gypsum were evaluated: 0; 3; 6; 9 
and 12 t ha”, in outline randomized block design 
with six repetitions. Soil density, total porosity, 
macroporosity, microporosity and saturated 
hydraulic conductivity were evaluated at layers 
of 0.0 - 0.1; 0.1 - 0.2 and 0.2 - 0.3 m. Soybean 
productivity and root growth were also evaluated. 
Data were submitted to regression analysis. The 
physical attributes soil density, macroporosity 
and saturated hydraulic conductivity did not differ 
significantly with the application of the gypsum 
doses in the 0.0-0.1 and 0.2-0.3 m layers. 
However, in the 0.1-0.2 m layer, due to pressures 
imposed by the machines and agricultural 
implements deforming the soil, the were 
significant differences in the physical attributes 
of the density, macroporosity and saturated 
hydraulic conductivity. There was no significant 
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SOYBEAN PRODUCTIVITY 


difference in grain productivity and root growth of 
soybean. 
KEYWORDS: Aluminum. Glycine max. Soil 
structure. 


IMPLICAÇÕES DE DOSES DE GESSO 
AGRÍCOLA EM ATRIBUTOS FÍSICO- 
HÍDRICOS DE UM LATOSSOLO 
ARGILOSO E NO CRESCIMENTO 
RADICULAR E PRODUTIVIDADE DA 
SOJA 
RESUMO: O objetivo deste estudo foi avaliar 
o efeito de doses de gesso nos atributos 
físico-hídricos do solo, crescimento radicular 
e produtividade da soja. O experimento foi 
realizado no Instituto Agronômico do Paraná 
(IAPAR) município de Santa Tereza do Oeste- 
PR. O solo foi classificado como Latossolo 
Vermelho distroférrico típico. Foram avaliadas 
cinco doses de gesso agrícola: 0; 3; 6; 9 e 12 
t ha!, em delineamento de blocos ao acaso 
com seis repetições. Avaliou-se a densidade 
do solo, porosidade total, macroporosidade, 
microporosidade e condutividade hidráulica 
saturada nas profundidades de 0,0 - 0,1; 0,1 
- 0,2 e 0,2 - 0,3 m. Também foram avaliados 
a produtividade e crescimento radicular da 
soja. Os dados foram submetidos à análise de 
regressão. Os atributos físicos densidade do 
solo, macroporosidade e condutividade hidráulica 
saturada não diferiram de forma significativa com 
a aplicação das doses de gesso nas camadas 
0,0 -0,1 e 0,2 - 0,3 m. No entanto, na camada 
de 0,1 - 0,2 m, devido à pressão imposta pelas 
máquinas e implementos agrícolas deformando 
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o solo, houve diferenças significativas nos atributos físicos da densidade, macroporosidade 
e condutividade do solo saturado. Não ocorreu diferença significativa para a produtividade e 
crescimento radicular da soja. 

PALAVRAS-CHAVE: Alumínio. Glycine max. Estrutura do solo. 


INTRODUCTION 


Agricultural gypsum is composed basically of calcium and sulfur and acts as a soil 
conditioner. High solubility, when applied to the soil, reduces aluminum saturation in depth, 
and translocate nutrients ffom the superficial layers to the subsurface, mainly calcium, thus 
allowing greater efficiency and area explored by the roots (LEITE et al., 2007; ZAMBROSI 
et al., 2007b; BROCH et al., 2008; RAIJ, 2008; SORATTO, CRUSCIOL, 2008b; NEIS et 
al., 2010; BROCH et al. (2011); ZANDONÁ et al., 2015). Gypsum could be applied to acid 
soils to complement limestone to favor the root system growth and development in depth, 
because of its soil conditioning effects (MEURER et al., 2004). 

According to Soratto, Crusciol (20084), it is large the amount of information on the 
agricultural gypsum effect to improve the radicle environment of plants, due to the calcium 
movement to subsurface layers of the soil or decrease in the toxic effects of high aluminum 
contents. Still for Soratto, Crusciol (2008b), subsurface layers with low calcium levels and/ 
or high exchangeable aluminum contents may cause decrease of harvests, especially in 
regions prone to low rainfall (dry spell), as they lead to less root system deepening, resulting 
in less soil volume explored by the roots, and in turn, less nutrients and water available to 
the plant. 

Within this optics to Raij (2008), the activity of free AI in the solution is a more 
consistent indicator of aluminum toxicity in soil solutions. Thus, the Ca displacement in the 
soil profile is much greater when gypsum is the source (CaSO,.2H,0). SO? anion that is 
part of the agricultural gypsum formula is important in reducing AI** activity, which helps 
root development (Zambrosi et al., 2007a; Nora et al., 2014). However, the incorrect use of 
plaster can cause chemical imbalance to the soil, causing long-term damage (PAULETTI 
et al., 2014). According to Costa et al. (2007), the reduction of soil density (Sd) can be 
attributed to the fact that the gypsum aggregation action is due to the supply of cations that 
possibly neutralize part of the negative charges occurring in the medium. 

For Sanchez et al. (2014), the maintaining soil with good fertility, physical, chemical, 
biological and water properties can provide good crop productivity. 

In western Paraná, soybean stands out in the economy, with much of the region's 
income coming from the culture (CONAB, 2017). The technology of the use of agricultural 
gypsum has been an important handling strategy to aid in grain production, mainly as soil 
conditioner for the soybean crop. 

Therefore, this study aimed to evaluate the effect of agricultural gypsum doses on 
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soil physical attributes, root growth and soybean productivity in a no-till system in west 
Paraná. 


MATERIAL AND METHODS 


The experimental area is in the municipality of Santa Tereza do Oeste-PR, located 
on the third Paraná plateau at latitude 25º08” (South) and longitude 53º58” (West), with an 
average elevation of 750 m. The experiment occurred in a Typic Haplortox, clayey texture, 
basalt substrate, smooth-wavy relief (EMBRAPA, 2013), cultivated under no-tillage system 
(NTS) for more than eighteen years. The soil is clayey textural class, with 29% sand, 11% silt 
and 60% clay. Table 1 shows the soil chemical determinations, collected in 2013, prior to the 
experiment, that is, before the application of the agricultural gypsum. Fifteen subsamples 
were collected in the experiment area. 

According to Kôppen's classification, the climatic characteristic of the region is 
mesothermal humid subtropical, Cfa, with average in the hottest month above 22ºC and 
in the coldest month below 18ºC, with no defined dry season, hot summer and frost 
less frequently. Rainfall is abundant and well distributed throughout the year; the annual 
averages are between 1,800 and 2,000 mm (IAPAR, 2000). During the whole cycle, the 
rainfall volume was around 1,118 mm and the average temperature was 22.61ºC. Figure 1 
shows the rainfall data, temperature and evapotranspiration evaluated every 15 days during 
soybean cycle. 

Initially, the area was cultivated with no-tillage system with winter and spring coverage; 
later, gypsum was applied to the soil surface (December 2013), to the following late growing 
season of corn (February 2014), beans in sequence (October 2014), late growing season 
of maize (February 2015), and finally soybean (October 2015) when the experiment started. 

The in outline was randomized block design with six repetitions. The plots (5 x 6 m) 
received the treatments that consisted of five doses of agricultural gypsum (0, 3, 6, 9, 12 
tha) in a randomized block design with six replicates. Agricultural gypsum dose was 3 t 
ha! (60% clay, 50 kg = 3.0 tha”). From the 3.0 t ha! value, the other gypsum doses were 
defined: zero, once, twice, three and four times the recommended dose, applied manually 
depending on the soil clay content Embrapa's recommendation EMBRAPA (2005). 

The cultivar BMX Apollo was utilized for the experiment and treated with insecticide 
TS Cruieser 350 FS 200 mL 100 kg seeds", and fungicides Vitavax-Thiram 300 mL100 kg 
of seeds”. Soybean sowing occurred on October 5, 2015, with 0.45 m spaced, 16 seeds 
per meter. The control of weeds, pests and diseases occurred according to the crop need. 
During August 2015, the area remained in fallow to collect the soil physical samples. To 
analyze soil density, soil microporosity, soil macroporosity, total porosity and saturated 
hydraulic conductivity, trenches were opened in each experimental unit and samples were 
taken from undisturbed samples on stainless-steel volumetric rings with a 98 cm? volume 
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(5 cm diameter and 5 cm height), in three soil layers (0.0 - 0.1; 0.1 - 0.2 and 0.2 - 0.3 m), 
with the aid of pedological hammer and soil extractor and there was no interaction between 
the factors. 

Soil volume in the samples was adjusted properly and saturated in water for 24 
hours in a tray with a water level at 2/3 of their height. Saturated samples were weighed 
and placed in a 0.6 m.c.a sand column and remained there for two days, draining the water 
in the macropores. After that we determined the saturated hydraulic conductivity (K,ç) in a 
constant-load permeameter. The samples were again saturated for 24 hours. Afterwards, this 
sample was allocated to the constant-load permeameter so that it could read the collected 
volume as a function of time, after which the soil hydraulic conductivity was calculated, 
according to the methodology recommended by EMBRAPA (2011). 

After the K,s determination, the samples were placed in an oven at 105ºC for 48 
hours to determine the dry soil mass. To calculate the soil density, the dry samples mass 
was divided by the volume of volumetric ring at 105ºC, according to the equation: Sd = Mss/ 
Tv, in which Sd is soil density and Mss, mass of the sample of dry soil at 105ºC, and Tv is 
total volume of the ring, EMBRAPA (2011). 

Calculating total soil porosity occurred with the relationship between soil density and 
particle density, as equation: Tp = (1 - Sd)/Pd *100 in which Tp is total porosity, measured 
in (mº mº) and Sd is soil density, measured in (Mg mº) and Pd is particle density measured 
in (Mg mº), EMBRAPA (2011). 

Microporosity was determined using water content retained in the samples in 
equilibrium with the tension of 0.6 m.c.a; this tension is enough to remove the water in the 
macropores, and the remaining water represents the volume of micropores obtained by 
equation: Micro = (Msu-Mss)/Mss *100, Micro = Soil microporosity (mº? m); Msu or Pa 0.6 
m.c.a = soil mass after 0.6 m.c.a tension (Mg); Mss = mass of dry soil (Mg). 

Thus, with total porosity and microporosity, calculating the macroporosity was 
possible (Camargo et al., 2009) according to equation: Macro = Pt-Micro [Macro = soil 
macroporosity (mº m3); Tp = total porosity (m3 mº); Micro = soil microporosity (m* m?3)], 
EMBRAPA (2011). 

To analyze soybean root growth, a 0,50 x 0,50 m, wooden grid was made with 
several nylon lines forming small 5 cm square strips on each side, which was used to 
estimate the root size and area. For the analysis, 1.0 x 0.6 m depth trenches were opened 
for each treatment, where soybean roots were exposed in parallel rows by trench. Thus, the 
roots were exposed for evaluation, and, with the squared grid placed in front of the root, the 
measurements of the soybean roots were carried out, profile method (Bôhm, 1979). 

The soybean harvest occurred on February 10, 2016 using a combine harvester 
(Winter steiger Classic?) with three 0.45 m rows. Productivity was determined by the useful 
area, 30 m?, transformed in kg ha”! and corrected to 13% moisture. 

Data were submitted to analysis of variance and the effect of the agricultural gypsum 
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doses evaluated by regression analysis. The model based on the coefficient significance 
of the adjusted regression equation as well as on the coefficient of determination (R?) 
associated with each model, using the Assistat software version 7.7 beta (SILVA, 2016). 


RESULTS AND DISCUSSION 


The significant results with the application of the gypsum doses on the physical 
properties of the soil occurred only in the layers of 0.1 to 0.2 m for all attributes worked. 
These results are usually expected in this layer, due to the pressures of the machines, 
agricultural implements and the influence of the chemical reactions of the soil in the 
presence of the gypsum. 

The area with the system of direct manipulation with the gypsum has contributed to 
the results of the processes of comparison with the root systems of the soybean crop, with 
the use of gypsum improved physical soil properties. When applying 3 t ha! gypsum, the 
best result of productivity corresponds to the recommendation of Embrapa (2015), which 
had a better effect. Thus, as a result, we can see Figures 2, 3 and 4 attributes. 

Soil density is the property representing the soil compaction state of the soil where, 
generally, in areas handled with NTS and agricultural gypsum, the highest values occur 
in the 0.1 - 0.2 m layer. In addition, the pressure imposed by agricultural machines and 
implements on the surface soils may negatively impact the root growth of soybeans. 

Soil density in the 0.0 - 0.1 m and 0.2 - 0.3 m layers did not differ significantly, 
whereas for the 0.1 to 0.2 m layer a significant difference occurred (Figure 2A, Be C). Soil 
density had an average value of 1.02 Mg m%. Treatments on the 0.0 - 0.1 m and 0.2-0.3m 
layers showed results close to the general average. For the 0.0 - 0.1 m layer, due to higher 
organic matter, root and biopor content, and because of greater soil rotation by seeder 
furrows, compacted layer is not formed. For the 0.1 m - 0.2 m layer, there was a significant 
increase in relation to the other layers, because agricultural gypsum improved the roots 
of the previous crops (Figure 2B). Pressures exerted by the machines and agricultural 
implements that promote deformation, compaction and root debris from previous crops 
promote increased soil density by occupying and pressing the soil porous spaces. The 0.2 - 
0.3 m layer has not changed, because the tire pressure of the agricultural machinery would 
not reach this layer. 

According to Reichert et al. (2003) the restrictive value of density for Haplortox with 
clay contents greater than 70% is close to 1.40 Mg m*. For Reichert et al. (2003) soil density 
values clayey soils of 1.35 Mg m* and for sandy soils of 1.55 Mg m* indicate compaction, 
which may hinder the root system development. 

In the present study soil density values were below those reported as limiting or 
potentially causing root growth difficulties. These values are similar to those reported by 
Araújo et al. (2007) found value 0.84 in between 1, 3 Mg mº. 
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For soil macroporosity in the 0.0 - 0.1 m and 0.2 - 0.3 m layers there was no 
significant difference (Figure 3A and 3C) depending on the gypsum doses. For the 0.1 to 0.2 
m layer, significant difference occurred (Figure 3B). 

In this study, the average macroporosity values were 16.46% for the 0.0 - 0.1 m 
layer and 15.19% for the 0.2 - 0.3 m layer. These values were higher than those considered 
critical to the crop development according to Jong Van Lier (2010). In the 0.1 - 0.2 m layer, 
the value was 17.79 %, which is above the critical value. Therefore, root growth of the 
previous cultures improved, which were benefited with the agricultural gypsum allowing an 
increase of macropores in this layer. Critical are 10% Reichert et al (2009). 

According to Jong Van Lier (2010), the critical value for the agricultural cultivation 
development is around 10% for soil macroporosity. Thus, all treatments present adequate 
amount of macropores (Figures 2, 3, 4). 

The saturated hydraulic conductivity in the 0.0 - 0.1 m and 0.2 - 0.3 m layers did not 
differ significantly. For the 0.1 to 0.2 m layer, the difference was significant (Figure 4B). As 
soil macroporosity was favored with the organic matter of the previous crops roots in this 
layer, the saturated hydraulic conductivity also had a significant effect. This higher Ks,s value 
associates with the higher volume of macropores in this layer (Figure 4C). 

There is a possible tendency for saturated hydraulic conductivity (K,ç) to increase in 
depth (Figure 4B). The results did not show significant differences between the 0.0 - 0.1 m 
and 0.2 - 0.3 m layers, whereas for the 0.1 - 0.2 m layer a significant difference occurred. 

Regarding the root growth, a root evaluation was performed at stage R8 (full 
maturation) for gypsum dose responses, in which there was no positive response to vertical 
and lateral root growth in the application of agricultural gypsum doses in the soybean crop 
(Figure 5). 

Caires et al. (2008) explains that the absence of soybean response to gypsum 
application may occur because of the growth of its soybean root system, in the absence 
of water deficit (Figure 2), not being influenced by the reduction of Al saturation in the soil 
subsurface. However, according to Raij (2008), applying agricultural gypsum allows better 
adequate conditions to the subsoil, and thus it can reduce the soil compaction, which is 
usually favorable to the roots. 

In the periods between 2015 and 2016, climatic conditions allied to management 
practices in the no-tillage system with surface layer rich with organic matter and a year with 
a good rain distribution during the entire vegetative cycle, without mechanical restriction, 
good soil structure in satisfactory physical conditions, for having a larger number of porous 
space allowed the root system to explore the studied soil layers. Good crop development is 
possible with no effect on the application of the agricultural gypsum. 

No significant differences between treatments for grain productivity were observed 
(Figure 6). 

There were probably no significant differences between the gypsum doses for the 
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grain productivity of soybean because of climatic conditions favorable to crop development, 
such as good fertility and structural conditions of the soil, good rainfall distribution, and 
presence of organic matter under the no-tillage system. 

Other authors also found no effect on grain productivity compared to gypsum doses, 
as in the works of Neis et al. (2010); Souza et al. (2010); Caires et al. (2011a); Caires et al. 
(2011b). There were probably no significant differences between the gypsum doses for the 
grain yield of soybean because of climatic conditions favorable to crop development, such 
as good fertility and structural conditions of the soil, good rainfall distribution (Figure 2), and 
presence of organic matter under the no-tillage system, as demonstrated by Rampim et al. 
(2011). 

In an application experiment of gypsum in the soybean crop, conducted in a 
dystrophic Red Latosol, Cardoso et al. (2014) observed the absence of soybean response 
to gypsum application, pointing to the non-occurrence of water deficit as a limiting factor to 
the gypsum action for this culture. 

Likewise, Caires et al. (2011a) did not obtain increase of grain productivity in the 
soybean crop when applying 9 t ha! gypsum in a dystrophic Red Latosol with clay texture, 
both in no-tillage system and in conventional tillage. 


CONCLUSIONS 


Were often evaluated to respond to the application of gypsum in this layer only the 
0.1 - 0.2 m layer presented a statistical difference for all the soil physical attributes. 


ACKNOWLEDGEMENTS 


To Capes for granting scholarship to the first author. 


REFERENCES 


ARAÚJO, R.; GOEDERT, W. J.; PINTO, M. P. C. Qualidade de um solo sob diferentes usos e sob 
cerrado nativo. Revista Brasileira de Ciência do solo, v. 31, n. 5, p. 1099-1108, 2007. 


BÔHM, W. Methods of studying root systems. Berlin: Springer-Verlag, 1979. 188p. 


BROCH, D. L.; NOLLA, A.; QUIQUI, E. M. D.; POSSENTI, J. C. Influência no rendimento de plantas de 
soja pela aplicação de fósforo, calcário e gesso em um latossolo sob plantio direto. RECEN-Revista 
Ciências Exatas e Naturais, v. 10, n. 2, p. 211-220, 2008. 


BROCH, D.L.; PAVINATO, P.S.; POSSENTTI, J.C.; MARTIN, T.N.; QUIQUI, E.M.D. Produtividade da 
soja no cerrado influenciada pelas fontes de enxofre. Revista Ciência Agronômica, v.42, n.3, p.791-796, 
2011. 


Engenharia Agronômica: Ambientes Agrícolas e seus Campos de Atuação 3 Capítulo 4 5 


CAIRES, E. F.; GARBUIO, F. J.; CHURKA, S.; BARTH, G.; CORRÊA, J. C. L. Effects of soil acidity 
amelioration by surface liming on no-till corn, soybean, and wheat root growth and yield. European 
Journal of Agronomy, v. 28, n. 1, p. 57-64, 2008. 


CAIRES, E. F.; JORIS, H. A. W.; CHURKA, S. Long-term effects of lime and gypsum additions on no-till 
corn and soybean yield and soil chemical properties in southern Brazil. Soil Use and Management, v. 
27,n.1, p. 45-53, 2011a. 


CAIRES, E. F.; MASCHIETTO, E. H. G.; GARBUIO, F. J.; CHURKA, S., JORIS, H. A. W. Surface 
application of gypsum in low acidic Oxisol under no-till cropping system. Scientia Agricola, v. 68, n. 2, p. 
209-216, 2011b. 


CAMARGO, O. A.; MUNIZ, A. C.; JORGE, J. A.; VALADARES, J. M. A. S. Métodos de análise química, 
mineralógica e física de solos do Instituto Agronômico de Campinas. Campinas: Instituto Agronômico, 
2009. 77 p. (Boletim técnico, 106, Edição revista e atualizada). Disponível em: <http:/Avww.iac.sp.gov. 
br/publicacoes/porassunto/pdf/Boletim106.pdf> Acesso em: 8 set. 2017. 


CARDOSO, J. A. E.; PERES, G. C. M.; LAMBERT, R.A. Influência da aplicação de calcário e gesso 
na cultura da soja (Glycine max (L.) Merril). Revista Enciclopédia Biosfera, v. 10, n. 18, p. 1980-1987, 
2014. 


COMPANHIA NACIONAL DE ABASTECIMENTO. Acompanhamento da safra brasileira de grãos. 
Brasília: v 10 Disponível em: < https:/Ayww.conab.gov.br/index.php/institucional/publicacoes/compendio- 
de-estudos-da-conab/item/2898-compendio-de-estudos-da-conab-v-10-pordutividade-da-soja > Acesso 
em: 23 abril 2018. 


COSTA, M. J. ROSA JUNIOR, E. J.; ROSA, Y. B. C. J.; SOUZA, L. C. F.:; ROSA, C. B. J. Atributos 
químicos e físicos de um Latossolo sendo influenciados pelo manejo do solo e efeito da gessagem. 
Acta Scientiarum. Agronomy, v. 29, n. 5, p. 701-708, 2007. 


EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA - EMBRAPA. Manual de métodos de análise 
de solo. Rio de Janeiro: Centro Nacional de Pesquisa de Solos, 2011. 230 p. 


. Sistema brasileiro de classificação de solos. Rio de Janeiro: Centro Nacional de Pesquisa de 
Solos, 2013. 353 p. 


INSTITUTO AGRONÔMICO DO PARANÁ - IAPAR. Cartas Climáticas do Estado do Paraná. Londrina: 
IAPAR, 2011. Disponível em: < http:/Avww.iapar.br/pagina-677.html > Acesso em: 8 set. 2017. 


JONG VAN LIER, Q. Física do Solo. Viçosa: Sociedade Brasileira de Ciência do Solo, 2010. 298 p. 


LEITE, E. M.; CAVALCANTE, L. F.: DINIZ, A. A.; SANTOS, R. V.; ALVES, G. S.; CAVALCANTE, |. H.L. 
Correção da sodicidade de dois solos irrigados em resposta à aplicação de gesso agrícola. Irriga, v. 12, 
n. 2, p. 168-176, 2007. 


MEURER, E. J.; RHENHEIMER, D.; BISSANI, C. A. Fenômeno de sorção em solos. In: MEURER, J. E. 
(Ed.). Fundamentos de química do solo. 2. ed. Porto Alegre: Gênesis, 2004. p. 131-179 


NEIS, L.; PAULINO, H. B.; SOUZA, E. D.; REIS, E. F.; PINTO, F.A. Gesso agrícola e rendimento de 
grãos de soja na região do sudoeste de Goiás. Revista Brasileira de Ciência do Solo, v. 34, n. 2, p. 409- 
416, 2010. 


Engenharia Agronômica: Ambientes Agrícolas e seus Campos de Atuação 3 Capítulo 4 EO 


NORA, D. D. AMADO, T. J. C.; BORTOLOTTO, R. P.; FERREIRA, A. O.; KELLER, C.; KUNZ, J. 
Alterações químicas do solo e produtividade do milho com aplicação de gesso combinado com calcário. 
Magistra, v. 26, n. 1, p. 1-10, 2014. 


PAULETTI, V.; PIERRI, L de.; RANZAN, T.; BARTH, G.; MOTTA, A. C. V. Efeitos em longo prazo da 
aplicação de gesso e calcário no sistema de plantio direto. Revista Brasileira de Ciência do Solo, v. 38, 
n. 2, p. 495-505, 2014. 


RAIJ, B. V. Gesso na agricultura. Campinas, Brasil: Instituto Agronômico/Fundação IAC, 2008. 233 p. 


RAMPIM, L. Atributos químicos de um latossolo vermelho eutroférrico submetido à gessagem e 
cultivado com trigo e soja em semeadura direta. Revista Brasileira Ciência do Solo, v. 35, n. 5, p. 1687- 
1698, 2011. 


REICHERT, J. M.; REINERT, D. J.; BRAIDA, J. A. Qualidade dos solos e sustentabilidade de sistemas 
agrícolas. Ciência Ambiental, v. 27, p. 29-48, 2008. 


SANCHEZ, E.: GENÚ, A. M.; MAGGI, M. F; MULLER, M. L. Propriedades físicas do solo e 
produtividade de soja em sucessão a plantas de cobertura de inverno. Revista Magistra, v. 26, n. 3, p. 
262-271, 2014. 


SILVA, F.A. S.; AZEVEDO, C. A. V. The Assistat Software Version 7.7 and its use in the analysis of 
experimental data. African Journal of Agricultural Research, v. 11, n. 39, p. 3733-3740, 2016. 


SORATTO, R. P.; CRUSCIOL, C. A. C. Atributos químicos do solo decorrentes da aplicação em 
superfície de calcário e gesso em sistema plantio direto recém-implantado. Revista Brasileira de 
Ciência do Solo, v. 32 n. 2, p. 675-688, 2008a. 


SORATTO, R. P.; CRUSCIOL, C. A. C. Produção de fitomassa e acúmulo de nutrientes pela aveia-preta 
em função da aplicação de calcário e gesso em superfície na implantação do sistema plantio direto. 
Ciência Rural, v. 38, p. 9228-935, 2008b. 


SOUZA, F. R.; ROSA JUNIOR, E. J.; FIETZ, C. R.; BERGAMIN, A. C.; VENTUROSO, L. R.; ROSA, 
Y.B. C. J. Atributos físicos e desempenho agronômico da cultura da soja em um Latossolo vermelho 
distroférrico submetido a dois sistemas de manejos. Ciência e Agrotecnologia, v. 34, n. 6, p. 1357-1364, 
2010. 


ZAMBROSI, F. C. B.; ALLEONI, L. R. F.; CAIRES, E. F. Aplicação de gesso agrícola e especiação 
iônica da solução de um Latossolo sob sistema plantio direto. Ciência Rural, v. 37, n. 1, p. 110117, 
.2007a. 


ZAMBROSI, F. C. B.; ALLEONI, L. R. F.; CAIRES, E. F. Nutrient concentration in soil water extracts and 
soybean nutrition in response to lime and gypsum applications to an acid Oxisol under no-till system. 
Nutrient Cycling in Agroecosystems, v. 79, n. 2, p. 169-179, 2007b. 


ZANDONÁ, R. R.; BEUTLER, A. N.; BURG, G. M.; BARRETO, C. F.: SCHMIDT, M. R. Gesso e 
calcário aumentam a produtividade e amenizam o efeito do déficit hídrico em milho e soja. Pesquisa 
Agropecuária Tropical, v. 45, n. 2, p. 128-137, 2015. 


Engenharia Agronômica: Ambientes Agrícolas e seus Campos de Atuação 3 Capítulo 4 





Layers AP Ca? Mg? k* V m 





em O cesacacecennanaama cmol, dm? ===... ooo nnnnnnn manaaaaaao Yo ====000.— 
0-10 0.68 4.20 1.083 0.25 36.09 11.21 
10-20 1.20 3.09 0.92 0.22 28.31 23.36 
20-40 1.33 2.14 0.91 0.19 23.28 30.62 
40-60 1.08 1.97 0.99 0.18 25.48 25.53 
60-80 0.36 1.82 1.30 0.19 33.83 9.37 





Table 1. Content of Al, Ca, Mg, K, base saturation (V) and aluminum (m) at different soil layers 


mms Rainfall(mm) — —s-Eto (mm day-1)  —m--Temperature (ºC) 
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Vegetative and reproductive stages of soybean plants 


Figure 1. Rainfall, temperature and evapotranspiration (ETc) of the crop in the soybean 
vegetative cycle ffom October 2015 to February 2016 in the western region - PR. 
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Figure 2. Density of the soil in the layers 0.0 - 0.1 (A), 0.1 - 0.2 (B) and 0.2-0.3m(C)asa 
function of agricultural gypsum doses. 
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Figure 3. Macroporosity of the soil in the layers 0.0 - 0.1 (A), 0.1 - 0.2 (B) and 0.2-0.3m (C) as 
a function of agricultural gypsum doses. 
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Figure 4. Soil hydraulic conductivity (Ksat) in the layers 0.0 - 0.1 (A), 0.1 - 0.2 (B) and 0.2 - 0.3 


m (C) as a function of agricultural gypsum doses. 
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Figure 5. Lateral and vertical growth of soybean roots as a function of agricultural gypsum 
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Figure 6. Grain productivity in relation to agricultural gypsum dosages applied in the soybean 
crop. 
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RESUMO: As substâncias húmicas (SHs) 
são compostos da matéria orgânica (MOS), 
que constituem uma parte significativa do C 
e N orgânico total do solo. Pode haver uma 
diminuição progressiva nas concentrações de 
SHs em solos que foram convertidos de mata 
nativa para sistemas de produção convencionais. 
Portanto, é fundamental que se adote sistemas 
conservacionistas de manejo do solo, como 
os sistemas agroflorestais (SAFs), a fim de 
promover a manutenção da qualidade do solo. 
O presente trabalho teve como objetivo avaliar 
as características da MOS influenciadas pelo 
uso do solo: sistema agroflorestal (SAF), pousio 
e mata nativa, no município de Bom Jesus, 
estado do Piauí. O sistema agroflorestal (SAF) 
regenerativo, do tipo quintal agroflorestal, foi 
implantado em outubro de 2016 em uma área 
experimental na Universidade Federal do 
Piauí, Campus de Bom dJesus-Pl. Adjacente 
à área de implantação do SAF, foi mantida um 
área com vegetação nativa e outra deixada em 
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PERMANENTE 


pousio. A vegetação nativa predominante são 
espécies florestais do bioma Cerrado. O solo foi 
classificado como Latossolo Amarelo distrófico 
típico, textura média. A amostragem do solo foi 
feita em quatro profundidades: 0-10, 10-20, 20- 
30 e 30-40 cm, e quatro repetições. Foi realizado 
o fracionamento químico da MOS pela técnica 
da solubilidade diferencial para obter as frações 
humina (HU), ácido húmico (AH), e ácido fúlvico 
(AF). Não houve interação entre o tipo de uso e 
profundidade do solo para nenhuma das SHs da 
matéria orgânica e para o carbono orgânico total 
do solo (P>0,05). Para a fração HU foi observado 
efeito simples (P<0,05) tanto para o tipo de uso 
quanto para a profundidade avaliada, enquanto 
para o carbono orgânico total a variação dos 
dados foi significativa somente em função 
da profundidade. Considerando somente a 
profundidade do solo não foi observada variação 
das frações AF e AH até a profundidade de 40 
cm. Já a camada superficial (0-10 em) do solo 
apresentou maior percentual da fração HU, 
seguida das profundidades de 10-20 e 20-30 
cm, enquanto a 30-40 cm a proporção desta 
fração foi significativamente menor. Conclui-se 
que a distribuição de carbono entre as diferentes 
frações pouco variou entre os diferentes usos 
do solo e profundidade. As frações, AF e AH 
foram estáveis e semelhantes para os três 
usos e profundidades. Por outro lado, a fração 
HU discriminou as mudanças no uso do solo. 
Em comparação com o solo da mata nativa, os 
usos SAF e pousio mantiveram praticamente os 
mesmos percentuais de SHs e carbono orgânico 
total, indicando que ambas as práticas são 
eficientes em manter a estabilidade da matéria 
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orgânica do solo. 
PALAVRAS-CHAVE: Agricultura conservacionista, Sistemas Agroflorestais (SAFs), 
Substâncias húmicas. 


ABSTRACT: Humic substances (SHs) are composed of organic matter (MOS), which make 
up a significant part of the total organic C and N of the soil. There may be a progressive 
decrease in the concentrations of SHs in soils that have been converted from native forest 
to conventional production systems. Therefore, it is essential to adopt conservationist 
soil management systems, such as agroforestry systems (SAFs), in order to promote 
the maintenance of soil quality. This study aimed to evaluate the characteristics of MOS 
influenced by land use: agroforestry system (SAF), fallow and native forest, in the municipality 
of Bom Jesus, state of Piauí. The regenerative agroforestry system (SAF), of the agroforestry 
type, was implemented in October 2016 in an experimental area at the Federal University of 
Piauí, Campus de Bom Jesus-PlI. Adjacent to the SAF implantation area, an area with native 
vegetation and another one that was set aside were maintained. The predominant native 
vegetation are forest species from the Cerrado biome. The soil was classified as a typical 
dystrophic Yellow Latosol, medium texture. The soil sampling was done at four depths: 0-10, 
10-20, 20-30 and 30-40 cm, and four replications. The chemical fractionation of the MOS was 
performed by the differential solubility technique to obtain the fractions humine (HU), humic 
acid (AH), and fulvic acid (AF). There was no interaction between the type of use and depth 
of the soil for any of the SHs of organic matter and for the total organic carbon of the soil (P> 
0.05). For the HU fraction, a simple effect (P <0.05) was observed for both the type of use and 
the depth assessed, while for the total organic carbon the variation in the data was significant 
only as a function of depth. Considering only the depth of the soil, no variation of the AF and 
AH fractions was observed up to a depth of 40 cm. The superficial layer (0-10 cm) of the soil 
presented a higher percentage of the HU fraction, followed by depths of 10-20 and 20-30 cm, 
while at 30-40 cm the proportion of this fraction was significantly lower. It is concluded that 
the carbon distribution among the different fractions varied little between the different uses 
of the soil and depth. The fractions, AF and AH were stable and similar for the three uses 
and depths. On the other hand, the HU fraction discriminated against changes in land use. 
Compared to the native forest soil, the SAF and fallow uses maintained practically the same 
percentages of SHs and total organic carbon, indicating that both practices are efficient in 
maintaining the stability of the soil's organic matter. 

KEYWORDS: Conservation Agriculture; Agroforestry Systems (SAFs); Humic substances. 


11 INTRODUÇÃO 


A matéria orgânica do solo (MOS) desempenha papel fundamental na manutenção 
da produtividade dos solos tropicais. Fazem parte da MOS substâncias húmicas e não- 
húmicas. Substâncias não húmicas são compostos lábeis com ciclagem relativamente 
rápida no solo (SCHMIDT et al., 2011). As substâncias húmicas (SHs) são compostos da 
matéria orgânica mais estáveis, que constituem uma parte significativa do C e N orgânico 
total do solo (MILORI et al., 2002). As SHs podem melhorar a capacidade de tamponamento 


do solo, aumentar a retenção de umidade e fornecer micronutrientes disponíveis às 
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plantas. Além disso, esses compostos também podem se ligar a metais, atenuando tanto 
a toxicidade dos metais pesados quanto a deficiência de metais nos solos (MCCARTHY, 
2001). 

As SHs incluem três frações distintas: fração humina é a fração insolúvel das SHs; 
o ácido húmico (AH) é a fração solúvel em condições alcalinas; e o ácido fúlvico (AF) é a 
fração solúvel em condições alcalinas e ácidas (SUTTON; SPOSITO, 2005). 

Diferentes usos do solo podem alterar as propriedades químicas das SHs do 
solo (SANTOS et al., 2013). Em geral, estudos relatam uma diminuição progressiva nas 
concentrações de SHs em solos que foram convertidos de mata nativa para sistemas de 
produção convencionais. Portanto, é fundamental que se adote sistemas de manejo do 
solo conservacionistas, como os sistemas agroflorestais (SAFs), que se assemelham as 
condições de mata nativa, a fim de promover a manutenção da qualidade do solo. 

No Nordeste do Brasil, principalmente na região do semi-árido, devido as altas 
temperaturas, o teor de MOS nos solos cultivados é baixo. Nesta região, existem poucos 
estudos que abordam o impacto da mudança do uso do solo nos níveis de SHs nos solos. O 
presente trabalho teve como objetivo avaliar as características da MOS como influenciadas 
pelo uso do solo: sistema agroflorestal (SAF), pousio e mata nativa, em um Latossolo 
Amarelo no município de Bom Jesus, estado do Piauí. 


21 METODOLOGIA 


A área avaliada consiste de um experimento com sistema agroflorestal (SAF) 
regenerativo implantado em outubro de 2016 em uma área experimental na Universidade 
Federal do Piauí, Campus de Bom Jesus, estado do Piauí (09º 04'26” Se 44º 21'32” 0). 
O SAF segue o princípio de um quintal agroflorestal, ou seja, composto por várias espécies 
agrícolas e florestais. Na mesma área de implementação do SAF, parte foi mantida com 
vegetação nativa e outra deixada em pousio. A vegetação nativa predominante são espécies 
florestais do bioma Cerrado. 

O solo é classificado como Latossolo Amarelo distrófico típico, textura média. O 
clima da região é do tipo Aw (Kôppen) com temperatura média anual de 26,5ºC, precipitação 
média anual de 945 mm, com estação chuvosa de outubro a abril, sendo de janeiro a março 
o trimestre mais chuvoso, com ocorrência de veranicos (PRAGANA, 2011). 

A amostragem do solo foi feita em quatro profundidades: 0-10, 10-20, 20-30 e 30-40 
cm. Foi realizado o fracionamento químico da MOS pela técnica da solubilidade diferencial 
para obter as frações humina (HU), ácido húmico (AH), e ácido fúlvico (AF), com base na 
solubilidade em meio ácido e alcalino. A determinação quantitativa de carbono nos extratos 
AF, AH e HU foram realizadas conforme Yeomans e Bremner (1988). 

Os dados foram submetidos à análise de variância, aplicando-se teste F, e quando 
o F foi significativo, as médias foram comparadas entre si pelo teste Skott-Knott a 5% de 
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probabilidade. As análises foram realizadas no programa Sisvar (FERREIRA, 2011). 


31 RESULTADOS E DISCUSSÃO 


Não houve interação entre o tipo de uso e profundidade do solo para nenhuma das 
substâncias húmicas (SHs) da matéria orgânica e para o carbono orgânico total do solo 
(P>0,05). Para a fração humina (HU) foi observado efeito simples (P<0,05) tanto para o 
tipo de uso quanto para a profundidade avaliada, enquanto para o carbono orgânico total a 
variação dos dados foi significativa somente em função da profundidade. 

A fração HU foi a única que apresentou efeito significativo do tipo de uso do solo 
(Figura 1). A igualdade no conteúdo de carbono orgânico total no solo, para os três tipos de 
uso, justifica as similaridades entre os percentuais de SHs observadas. A maior expressão 
da fração HU pode ser verificada na Figura 1 onde observa-se também que as frações 
ácido fúlvico (AF) e ácido húmico (AH) contribuem de forma equivalente. 
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Figura 1. Média das substâncias húmicas da matéria orgânica do solo e carbono orgânico total 
para três tipos de uso do solo no município de Bom Jesus-PI. 


Considerando somente a profundidade do solo não foi observada variação das 
frações AF e AH até a profundidade de 40 cm (Figura 2). Já a camada superficial (0-10 cm) 
do solo apresentou maior percentual da fração HU, seguida das profundidades de 10-20 e 
20-30 cm, enquanto à 30-40 cm a proporção desta fração foi significativamente menor. E, 
uma vez que a maior parte do carbono orgânico do solo é constituída por HU, o carbono 
orgânico total do solo exibiu variação proporcional. No entanto, a profundidade de 0-10 cm 
não diferiu da subsequente (10-20 cm), mas, estas duas possuem conteúdo de carbono 
orgânico significativamente maior que a profundidade de 20 a 40 cm (Figura 2). 

Pode-se dizer que o processo de humificação depende do uso do solo à medida 
que este determina a quantidade de resíduos vegetais mantidos na superfície do solo. 
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Podemos inferir, portanto, que os três tipos de uso avaliados neste trabalho resultaram 
em processos similares de transformação dos resíduos orgânicos, levando a formação de 
percentuais equivalentes de SHs e carbono orgânico total. 

Uma das principais finalidades da avaliação das frações da MOS é com relação a 
modificação na fertilidade do solo. As SHs atuam significativamente na disponibilidade de 
nutrientes para as plantas em consequência da sua capacidade de se ligar à fração mineral 
do solo através de cátions polivalentes (CANELLAS et al., 1999). 

Por fim, conclui-se que o uso do solo com o SAF está sendo uma alternativa 
sustentável, uma vez que está mantendo o conteúdo de carbono orgânico do solo análogo 
ao da área de mata. Uma outra alternativa viável é o pousio, que também apresentou 
percentuais de substâncias húmicas e de carbono orgânico similares tanto à mata quanto 
ao SAF. No entanto, essa condição foi alcançada após um período de quatro anos e, na 
prática, dificilmente deixa-se uma área passível de produção ociosa por tanto tempo. 

A camada superficial de O a 20 cm foi a mais modificada com o processo de 
humificação, verificado pelo maior conteúdo de carbono orgânico e, consequentemente da 
fração HU. Isto se deve ao maior aporte de resíduos orgânicos sob a superfície do solo e 
maior atividade biológica nesta camada. 
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Figura 2. Média das substâncias húmicas da matéria orgânica do solo e carbono orgânico total 
para quatro profundidades do solo no município de Bom Jesus-PI. 


41 CONCLUSÃO 

A distribuição de carbono entre as diferentes frações pouco mudou nos diferentes 
usos do solo e profundidade. As frações ácido fúlvico e ácido húmico foram estáveis e 
semelhantes para os três usos e as todas as profundidades. Por outro lado, a fração 
humina foi mais eficaz em discriminar mudanças no uso do solo, devido à sua labilidade 
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característica. Em comparação com o solo da mata nativa, os usos SAF e pousio mantiveram 
praticamente os mesmos percentuais de substâncias húmicas e carbono orgânico total 
mostrando que ambas as práticas são eficientes em manter a estabilidade da matéria 
orgânica do solo. 
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RESUMO: O objetivo deste trabalho foi calibrar 
o método de dissipação térmica na estimativa 
da transpiração de plantas de pinhão manso. 
Plantas de pinhão manso foram cultivadas em 
quatro lisímetros de pesagem, sendo que estes 
tiveram sua superfície coberta por lona plástica 
para que a única saída de água do sistema 
fosse por transpiração. Sondas de dissipação 
térmica foram inseridas nas plantas a 20 cm 
da superfície do solo e aquecidas com uma 
potência de 0,1 W. Para o aquecimento foram 
construídas fontes de corrente ajustáveis, com 
possibilidade de conexão de até seis SDT, 
com os elementos aquecedores conectados 
em série. As plantas de pinhão-manso foram 
submetidas a 15 dias com fornecimento de água 
até a capacidade de campo do solo, 60 dias de 
supressão da irrigação e retorno da irrigação por 
mais 15 dias. Verificou-se desvio médio entre o 
fluxo de seiva e a transpiração em escala diária 
para a equação geral de Granier de 5,04%, e o 
desvio para a equação modificada foi de 4,88%. 
Diante dos resultados, conclui-se que o método 
de dissipação térmica pode ser utilizado para a 
determinação da transpiração de pinhão-manso. 
PALAVRAS-CHAVE: Evapotranspiração, 
agricultura irrigada, potencial hídrico. 


Capítulo 6 


CALIBRATION OF THE THERMAL DISSIPATION METHOD IN THE 
MEASUREMENT OF SAP FLOW IN PINHAO-MANSO 


ABSTRACT: The objective of this work was to calibrate the thermal dissipation method to 
estimate the transpiration of physic nut plants. Jatropha plants were grown in four weighing 
lysimeters, and their surface was covered with plastic tarpaulin so that the only water outlet 
in the system was through transpiration. Thermal dissipation probes were inserted into the 
plants at 20 cm from the soil surface and heated with a power of 0.1 W. For heating, adjustable 
current sources were built, with the possibility of connecting up to six SDT, with the heating 
elements connected in series. Jatropha plants were submitted to 15 days with water supply 
until the field capacity of the soil, 60 days of suppression of irrigation and return of irrigation 
for another 15 days. There was an average deviation between sap flow and transpiration on 
a daily scale for the general Granier equation of 5.04%, and the deviation for the modified 
equation was 4.88%. In view of the results, it is concluded that the thermal dissipation method 
can be used to determine the physic nut sweat. 

KEYWORDS: Evapotranspiration, irrigated agriculture, water potential. 


INTRODUÇÃO 


Atranspiração é o processo pelo qual a água é absorvida pelo sistema radicular das 
plantas, atravessando-as pelos vasos condutores (xilema) e chega até as folhas onde ocorre 
sua saída para a atmosfera através dos estômatos. De forma resumida, é a difusão da água 
através dos estômatos para a atmosfera (TAIZ et al., 2017). A estimativa da transpiração 
pode ser realizada por meio do método de dissipação térmica, ou simplesmente Granier 
(GRANIER, 1985). envolve a inserção perpendicular de duas sondas na mesma linha axial 
do ramo ou tronco, sendo que uma delas, a superior, é aquecida de forma constante com 
fornecimento de potência elétrica de aproximadamente 0,2 watt. A variação da temperatura 
em cada instante entre as duas sondas (a aquecida e a não aquecida) é ocasionada pelo 
transporte convectivo de calor por meio da seiva. 

A máxima diferença de temperatura entre as duas sondas significa que o fluxo de 
seiva é mínimo ou nulo, enquanto que a mínima diferença significa uma taxa máxima de 
fluxo de seiva através da área de transporte (xilema). 

O método de Granier é amplamente utilizado para determinação da transpiração das 
mais diversas culturas. O método de Granier é amplamente utilizada para determinação 
da transpiração das mais diversas culturas. Em cultivo de manga (LU; CHACKO, 1998) 
encontraram boa concordância entre os valores de transpiração medidas por gravimétrica 
e o fluxo de seiva em escala diária pelo método de dissipação térmica. 

O método de Granier fundamenta-se em uma calibração para se ter o valor do fluxo 
de seiva, ao contrário dos outros métodos térmicos que se baseiam em um tratamento 
teórico, embora no método de pulso de calor possa haver uma calibração adicional ao 
tratamento teórico. Inicialmente o método foi desenvolvido para espécies florestais. O 
próprio Granier; Gross (1987) encontrou relações próximas à do seu primeiro trabalho para 
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outras espécies florestais, o que levou à proposta de que, pela proximidade dos valores 
entre os diferentes trabalhos, a equação de Granier poderia ser adotada como universal, 
pelo menos para plantas lenhosas. 

Sendo assim, o objetivo deste trabalho foi o de calibrar o método de dissipação 
térmica para estimativa da transpiração de plantas de pinhão-manso. 


MATERIAL E MÉTODOS 


Caracterização da área experimental 


O experimento estudo conduzido na área experimental do Departamento de 
Engenharia de Biossistemas da Escola Superior de Agricultura “Luiz de Queiroz” (ESALQ/ 
USP), em ambiente protegido (casa de vegetação), localizada no município de Piracicaba 
- SP (latitude 22º43'31”; longitude 47º38'57”; altitude 547 m). A transpiração das plantas 
foi estimada por medidas de fluxo de seiva, onde optou-se pelo uso da do método de 
dissipação térmica, ou comumente conhecido como Método de Granier Granier (1985). O 
motivo da escolha do método se baseia no baixo custo de obtenção dos sensores e também 
por exigir um menor número de canais de saída dos sistemas automático de aquisição de 
dados, dessa forma, o método viabiliza a medição simultânea de um maior número de 
plantas. Para coleta e armazenamento dos dados, utilizou-se um sistema automático de 
aquisição de dados, com leituras armazenadas a cada um segundo e armazenadas as 
médias a cada 15 min. 


Método de dissipação térmica 

O método de dissipação de térmica envolve a inserção perpendicular de duas sondas 
(Figura 1) na mesma linha axial do ramo ou tronco, sendo que uma delas, a superior, é 
aquecida de a uma potência constante. A variação da temperatura em cada instante entre 
as duas sondas (a aquecida e a não aquecida) é ocasionada pelo transporte convectivo de 
calor através da seiva. A máxima diferença de temperatura entre as duas sondas significa 
que o fluxo de seiva é mínimo ou nulo, enquanto que a mínima diferença significa uma taxa 


máxima de fluxo de seiva através da área de transporte (xilema). 
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Figura 1. Esquema do sensor de Granier inserida perpendicularmente no tronco de uma árvore. 
Observa-se a conexão dos termopares ao sistema de aquisição de dados e da sonda superior 
a uma fonte de energia. 


O fluxo de seiva foi determinado pela equação proposta por Granier (1985), 
apresentada na equação 1: 


F=118,99x10ºxk!23!xSA (1) 


sendo: SA a área do lenho condutor da seiva bruta, dada em m?, normalmente 
considerada a área ocupada pelo xilema e k é o índice de fluxo calculado por meio da 
equação 2, a partir das medidas térmicas registradas por meio do sensor. 
RE (ATM- AT) (2) 
AT 
sendo: ATM (ºC) a diferença máxima de temperatura entre os dois pontos de medida 
e AT (ºC) a diferença de temperatura real. 


Construção e calibração das sondas de dissipação térmica (SDT) 


A sonda dissipadora de calor foi confeccionada com fio de constantan de 50 um de 
diâmetro enrolado em volta de uma agulha de aço inoxidável de 0,9x25 mm, de tal forma que 
após o enrolamento formasse um cilindro de 1,2 mm de diâmetro e um comprimento de 10 
mm. Em cada agulha foi inserido, em sua parte interior média, uma junção de termopar feita 
da união de fios de cobre e com o fio de constantan. Os dois termopares, unidos em série 
e espaçados em 10 cm são responsáveis pelo fornecimento da diferença de temperatura 
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(Tc-Tb), a qual foi registrada por um sistema automático de aquisição de dados. 

As SDT foram inseridas a 20 cm do nível do solo e espaçadas (as agulhas) em 10 
cm. Os furos para inserção da cápsula de alumínio foram feitos com furadeira, utilizando- 
se uma broca de 1,7 mm. Tomou-se o cuidado de utilizar calda bordalesa para evitar a 
infecção do interior da planta com qualquer tipo de fungo. 

Para inserir as sondas notronco, previamente foi introduzida uma cápsula de alumínio 
com diâmetro de 1,67 mm e comprimento de 12 mm. Esta foi inserida duas semanas 
antes da instalação das sondas, afim de que a seção do caule perfurada tivesse tempo 
suficiente para cicatrização e eliminação do estresse causado por esse procedimento. 
Antes da instalação das SDT, foi pincelado nas agulhas uma pasta térmica para auxiliar 
na dissipação e uniformização do calor em volta da mesma. Para garantir que uma boa 
inserção da SDT foi aplicada uma camada de silicone externamente. Após a instalação das 
sondas um papel alumínio foi utilizado para cobrir a maior parte do tronco afim de minimizar 
o efeito da incidência da radiação eletromagnética sobre as medidas. 

O aquecimento da sonda superior ocorreu por meio da resistência criada pelo fio de 
constantan enrolado na parte externa da agulha. Para o aquecimento foram construídas 
fontes de corrente ajustáveis, com possibilidade de conexão de até seis SDT, com os 
elementos aquecedores conectados em série. A resistência elétrica total das sondas foi 
medida por meio de um multímetro e calculada a tensão a ser aplicada de acordo com, a 
lei de Ohm, de modo a se obter uma potência 0,1 W. 

As fontes de corrente foram montadas em uma caixa PB203, onde colocou-se a 
placa de circuito elétrico com os componentes necessários para funcionamento da fonte de 
corrente. Para que fosse possível a montagem de uma fonte de corrente com seis saídas 
em série, foram utilizados seis LM317 TO-220, seis dissipadores DM812 MD, seis trimpots 
MV 2K, seis resistores 1W 22R, um diodo 1N4007, um fusível, um resistor 1/4W 1K2, um 
led 5MM vermelho, dois conectores PTR doze vias, sendo seis machos e seis fêmeas 
catorze e dois conectores PTR 2 vias, sendo um macho e um fêmea. A fonte de corrente foi 
alimentada com voltagem de 12V. 


Calibração das sondas de dissipação térmica e avaliação das plantas 


Plantas de pinhão-manso utilizadas na calibração das SDT foram cultivadas 
em quatro lisímetros de pesagem. Os lisímetros foram confeccionados em vasos com 
capacidade para 250 L, com 0,75 m?, com sistema de pesagem direta, utilizando-se células 
de carga que suportavam diretamente todo o peso do sistema. Estes foram preenchidos 
com camada de pedra brita (em torno de 8 cm), seguida de manta geotêxtil e, por fim, o 
volume restante foi completado com solo. Para evitar perdas por evaporação, os lisímetros 
foram cobertos com lona plástica, de modo que a única saída de água seria por transpiração 
das plantas. 

As plantas localizadas nos quatro lisímetros foram utilizadas para monitorar a 
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transpiração durante 60 dias. Durante os primeiros 15 dias o solo foi mantido com umidade 
próxima da condição de capacidade de vaso (CV) por meio de irrigações frequentes. 
Posteriormente, durante 30 dias a aplicação de água via sistema de irrigação foi suspensa 
e, por fim, os 15 dias finais tiveram retorno da irrigação para a recuperação da hidratação 
das plantas. Todos os dados foram coletados por um sistema de aquisição de dados 
automático. 

Após a calibração das sondas, as mesmas foram instaladas em área experimental 
cultivada com pinhão-manso. Sendo então avaliado o efeito da poda na transpiração das 
plantas de pinhão-manso, com adoção de três diferentes podas de topo, com variações 
na altura de corte e no diâmetro de copa, dessa forma, os tratamentos adotados foram: 
Tratamento 1 (P1) — sem poda; Tratamento 2 (P2) — poda a 1,5 m de altura e 1,5 m de 
diâmetro de copa e; Tratamento 3 (P3) — poda a 1,0 m de altura e 1,0 m de diâmetro de 
copa. 


RESULTADOS E DISCUSSÃO 


Calibração do fluxo de seiva 


Como padrão na calibração do método de dissipação térmica foram utilizados os 
lisímetros de pesagem descritos anteriormente. Para facilitar os cálculos, utilizou-se os 
dados de saída dos lisímetros em escala horária, dessa forma, optou-se por modificar o 
coeficiente linear da equação de Granier, apresentada na equação 1, obtendo-se a seguinte 
equação: 


FS=0,000594xK!2º!xaS (3) 


Na Figura 2 é apresentada a relação entre a transpiração determinada por meio 
dos lisímetros e o fluxo de seiva calculado pela equação geral de Granier e a equação 
modificada (Eq. 3). O desvio médio entre o fluxo de seiva e a transpiração em escala diária 
para a equação geral de Granier foi de 5,04%, e o desvio para a equação modificada foi 
de 4,88%. Os desvios obtidos são aceitáveis, uma vez que a transpiração é um processo 
dinâmico, apresentando grande variabilidade espacial. De maneira geral, considera-se 
que, quando se está de posse de grande número de dados, desvios da ordem de 10% são 
perfeitamente aceitáveis (JUNIOR, 2012). 
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Figura 2. Relação entre as medidas de transpiração realizadas com os lisímetros e as 
determinação do fluxo de seiva diário pela sonda de dissipação térmica. 


Silva et al. (2004) verificaram em plantas de kiwi que a equação geral de Granier 
apresentou uma forte subestimativa dos valores de fluxo de seiva principalmente em 
condições de maior fluxo. Os autores ajustaram uma equação polinomial de grau 2 aos 
valores de k em relação a densidade de fluxo calculada com base em medidas lisimétricas 
e na área da secção condutora do caule, obtendo um coeficiente de determinação de 
0,8252. Rojas (2003) também baseado em medidas lisimétricas verificou que o fluxo de 
seiva pelo método da dissipação térmica subestimou na ordem de 3% a transpiração 
em plantas cítricas. Ferreira; Silvestre (2004) ajustaram uma equação potencial com os 
valores de fluxo de seiva obtidos com os parâmetros originais da equação de Granier e a 
evapotranspiração para vinha em condições de evaporação do solo desprezível, estimada 
pelo método meteorológico da covariância de fluxos turbulentos, obtendo uma boa 
correlação (R2 = 0,84). Vellame (2010) constataram que o método da sonda de dissipação 
térmica, utilizando-se os coeficientes originais da equação de Granier (1985) subestimou 
em média 31% os valores de fluxo de seiva em escala diária. 


Transpiração das plantas 


Na Figura 3 é apresentado o curso do fluxo de seiva estimado pela sonda de 
dissipação térmica e a radiação solar global ao longo ao longo de cinco dias, nos diferentes 
tratamentos de poda adotados. 
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Figura 3. Transpiração estimada pela sonda de dissipação térmica e radiação solar global ao 
longo de cinco dias de avaliação nos tratamentos de poda P1 (A), P3 (B) e PZ(C). 


Nota-se que a transpiração apresenta uma variabilidade entre os tratamentos de 
poda, sendo maior para as plantas com o tratamento de poda 2 (poda de 1,5 m de altura 
por 1,5 m de largura). 

O tratamento de poda P2, por ser o de poda menos drástica, pode ter apresentado 
maior área foliar, este fato explica uma maior transpiração por parte das plantas deste 
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tratamento, pois há uma maior superfície evaporante da planta voltada para atmosfera. 
Aliada à área foliar, a poda aumenta a turbulência da folha com a atmosfera, promovendo 
uma diminuição da camada limítrofe e consequente aumento na perda de água. 

Gush (2008) estudou a transpiração de plantas de pinhão-manso conduzidas sem 
irrigação e em áreas com solos pouco férteis. O autor obteve medidas do fluxo de seiva 
de plantas jovens (4 anos) e maduras (12 anos) para avaliar o consumo hídrico destas. 
As medidas foram realizadas por período de 17 meses usando o método da razão de 
calor (BURGESS et al., 2001). Segundo Gush (2008), esse método foi escolhido devido a 
sua habilidade para medir taxas pequenas de fluxo, o que costuma ser característico das 
plantas que crescem sem irrigação nesses locais avaliados. Considerando um período de 
um ano, o autor encontrou fluxo de seiva de, aproximadamente, 1 940 L pl para plantas 
jovens e 4 460 L pl” para plantas maduras. Esses valores foram extrapolados da escala 
de uma planta para de lavoura considerando população de 740 pI ha”, o que representou 
consumo hídrico de, aproximadamente, 144 e 330 mm ano” para plantas jovens e maduras, 
respectivamente. O autor considerou que o fluxo de seiva variou com o tempo devido em 
grande parte a natureza decidual da espécie, embora também haja influência da duração 
do dia e da precipitação. Valores diários foram bastante influenciados pela disponibilidade 
energética (radiação solar), pelo déficit de pressão de vapor atmosférico e pelos eventos 
de chuva que restituífam a umidade do solo. Para plantas jovens, o maior valor de fluxo de 
seiva diário foi de 32,9 L dia! (2,4 mm dia”), enquanto que a média diária ao longo do ano 
foi de, aproximadamente, 5,3 L dia! (0,4 mm dia'!). Para plantas maduras, esses valores 
foram 79,6 e 12,3 L dia” (5,9 e 0,9 mm dia”), respectivamente. 

Na Figura 4, o valor máximo do fluxo de seiva foi alcançado quando o déficit de 
pressão de vapor atingiu 1,47 kPa. Apesar da máxima transpiração não ter sido verificada 
com o máximo DPV, percebe-se uma característica de proporcionalidade direta entre as 
curvas, ou seja, à medida que o DPV aumenta ocorre um aumento na transpiração das 
plantas, sendo a transpiração diminuída com a diminuição do DPV. Sendo a planta de 
pinhão manso bastante resistente à seca, possivelmente, com boa disponibilidade de água 
no solo, a mesma continua com os estômatos abertos, mesmo nos períodos mais quentes 
do dia. 
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Figura 4. Comportamento diário da transpiração e DPV (da planta de pinhão manso. 


Os valores de fluxo de seiva ao longo do dia em função do DPV podem ser 
observados na Figura 52, por meio das setas. 
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Figura 5. Relação entre transpiração e DPV para a planta de pinhão manso. 


Analisando a Figura 5 percebe-se que a transpiração após a ascensão se normaliza 
e posteriormente decresce, sendo possível verificar uma tendência de histerese na curva. 
Possivelmente com a avaliação sendo feita para todo o dia, a curva se fecharia. Verificou- 
se valor máximo de transpiração de 2,5 L h-1 ocorrendo ao meio dia com DPV de 1,47 kPa 


e valor mínimo de transpiração de 1,3 L-h”, verificado às seis horas da manhã e DPV de 
0,9 kPa. 


CONCLUSÕES 


Nas condições em que o estudo foi realizado, conclui-se que o método de dissipação 
térmica pode ser utilizado para a determinação da transpiração de pinhão-manso. 
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RESUMO: O desenvolvimento dos grãos do 
cafeeiro pode vir a ser prejudicado em função 
das condições climáticas e hídricas. Com isso, 
objetivou-se nesse trabalho, avaliar o tamanho 
de grãos de cafeeiro conilon “Vitória”, em 
diferentes manejos de irrigação. O experimento 
foi conduzido em um esquema de parcela 
subdividida 3 x 13, em um delineamento em 
blocos casualizados com três repetições. Os 
dados experimentais foram submetidos à análise 
de variância e as médias comparadas pelo 
teste de Tukey para comparações entre o fator 
manejo e Scott-Knott para comparações entre o 
fator clone. A irrigação influenciou no tamanho 
dos grãos, sendo os maiores valores de grãos 
retidos em peneira 13 e superiores, encontrados 
nos manejos irrigados e os maiores valores 
fundagem em plantas sem irrigação. 
PALAVRAS-CHAVE: Coffea Canephora, água, 
tamanho, déficit hídrico, qualidade física. 


Capítulo 7 


CLASSIFICATION OF CONILON COFFEE GRAINS IN SIEVE UNDER 
IRRIGATED MANAGEMENT AND DRYLAND 


ABSTRACT: The development of coffee beans may be impaired due to climatic and water 
conditions. Thus, the objective of this work was to evaluate the size of conilon coffee beans 
“Vitória”, in different irrigation managements. The experiment was carried out in a 3 x 13 split 
plot scheme, in a randomized block design with three replications. The experimental data were 
subjected to analysis of variance and the means compared by the Tukey test for comparisons 
between the management factor and Scott-Knott for comparisons between the clone factor. 
Irrigation influenced the size of the grains, with the highest values of grains retained in sieve 
13 and higher, found in irrigated managements and the highest values founding in plants 
without irrigation. 

KEYWORDS: Coffea Canephora, water, size, water deficit, physical quality. 


11 INTRODUÇÃO 


O Brasil é um principal produtor e exportadores de cafés do mundo, sendo o segundo 
maior consumidor da bebida no mundo. A produção de café conilon, serve como base para 
o preparo de cafés solúveis que são consumidos em vários países, tornando este produto 
extremamente importante no mercado nacional e internacional (USDA, 2017; USDA, 2018; 
CONAB, 2019). 

O estado do Espírito Santo possui um grande destaque nesse cenário, pois além de 
ser o segundo maior produtor de café do país com 22,3% da produção, destaca-se também, 
como o responsável por 63,4% de todo o café conilon brasileiro (CONAB, 2019). Entretanto, 
devido à má distribuição de chuvas, condições climáticas não favoráveis e restrições 
hídricas, o desenvolvimento do café conilon pode vir a ser prejudicado, principalmente na 
fase de floração e enchimento dos grãos. 

O déficit hídrico é um dos principais fatores responsáveis pela queda produção, 
sendo também, responsável pela queda da qualidade, pois o déficit hídrico no cafeeiro, 
após a fecundação pode prejudicar o crescimento dos frutos (CAMARGO, 2010). Sendo 
assim, muitos produtores têm adotado a prática da irrigação com o intuito de fornecimento 
de água para melhor desenvolvimento da cultura. Sendo essa prática, quando bem adotada 
e implementada, pode vir a melhorar a produção e a qualidade da produção (FERRÃO 
et al., 2012). Porém, em muitas das áreas cultivadas, não é feito um manejo correto da 
irrigação, o que pode resultar em aplicações insuficientes ou excessivas de água, podendo 
fazer com que as plantas apresentem respostas diferentes do esperado. 

O tamanho do grão é fortemente influenciado pelas condições hídricas da planta, 
visto que em condições adequadas, ocorre a maior expansão do fruto, resultando em 
maior tamanho e melhor tipo (REZENDE et al., 2006). Entretanto, são poucos os estudos 
disponíveis que associem o manejo da irrigação com a classificação física dos grãos de 
café. 
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O presente trabalho teve por objetivo, realizar a classificação física via peneira dos 
grãos de cafeeiro conilon variedade “Vitória, em diferentes manejos de irrigação. 


21 METODOLOGIA 


O experimento foi conduzido em uma lavoura de 0,42 ha cultivada de Coffea 
canephora Pierre, variedade 'Conilon Vitória Incaper 8142”, constituída por treze clones, 
localizada no Instituto Federal do Espírito Santo - Campus de Alegre, Fazenda Caixa 
D'Água, distrito de Rive, localizado na latitude 20º25'53” S e longitude 41º27'25” W, altitude 
média de 137 m e precipitação média anual de 1250 mm. O solo do local é classificado 
como Latossolo Vermelho Amarelo, de textura argilo-arenosa. 

O experimento foi montado em um esquema de parcela subdividida 3 x 13, sendo 
nas parcelas o fator manejo de irrigação em três níveis (irrigado com 100% da ETc, irrigado 
com 50 % da ETc e sem irrigação) e nas subparcelas o fator clone em treze níveis (V1; 
V2; V3; V4; V5; V6; V7; V8; V9; V1O; V11; V12 e V13), em um delineamento em blocos 
casualizados, com três repetições. 

Para realização da irrigação, foi utilizado um sistema de irrigação por aspersão 
convencional, constituído por três linhas laterais, apresentando cada uma, quatro 
aspersores setoriais, espaçados por 18 x 18 m, com bocais 5,6 x 3,2 mm e vazão de 2,66 
mº h'!, sendo o Coeficiente de Uniformidade de Christiansen (CUC) de 81,0%. 
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Figura 1- Esquema do Experimento. 
Fonte: Pereira (2015) 


Para a realização do manejo da irrigação, foi adotado o manejo via clima. Sendo 
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a evapotranspiração de referência (ETo) estimada pelo método de Hargreaves e Samani 
(1985). Para obtenção da ETc (Evapotranspiração da Cultura) foi utilizado um Kc de 1,1. Na 
figura abaixo (Figura 2) é possível observar o gráfico com os valores de evapotranspiração 
de referência, evapotranspiração da cultura, precipitação e lâminas de irrigação aplicadas 
no período de outubro de 2019 a fevereiro de 2020. 
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Figura 2 — Precipitações (P), evapotranspiração de referência (ETo), evapotranspiração 
da cultura (ETc) e lâminas de irrigação (IRN) aplicadas no período de outubro/2019 a 
fevereiro/2020. 


A classificação por peneira foi realizada a partir de 300 g de amostra de café 
beneficiado, segundo as dimensões dos crivos, sendo numeradas de 10, 11 e 12 para 
grãos moca, sendo determinados os percentuais de grãos retidos nas respectivas peneiras, 
seguindo-se a Instrução Normativa Nº 8, de 11 de junho de 20083 (BRASIL, 2003). 

Os dados obtidos no experimento foram submetidos à análise de variância e quando 
significativo, as médias foram comparadas pelo teste de Scott-Knott para comparação 
entre o fator clone e Tukey para o fator manejo, ambos em 5% de probabilidade, utilizando 
o software R (R CORE TEAM, 2020) com auxílio do pacote ExpDes.pt (FERREIRA; 
CAVALCANTI; NOGUEIRA, 2018). 


31 RESULTADOS E DISCUSSÃO 

Pela análise de variância para as variáveis: moca, peneira 13 e superiores e 
fundagem. Verificou-se que houve efeito significativo para a interação manejo x clone, em 
nível de 5% de probabilidade, devendo-se proceder à análise de forma desdobrada. 

Nas Figuras abaixo, encontram-se o percentual de grãos moca, peneiras 13 e 
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superiores e fundagem, em função do manejo de irrigação e clone da variedade *Conilon 
Vitória”. 
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Figura 3 - Percentual de grãos moca em função dos manejos de irrigação para cada clone do 
cafeeiro *Conilon Vitória" 


Médias seguidas de letras iguais não diferem entre si pelo teste de Tukey em nível de 5% de 
probabilidade. 
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Figura 4 - Percentual de grãos moca em função dos treze clones do cafeeiro 'Conilon Vitória” 
para cada manejo (irrigado com 100% da ETc, irrigado 50% da ETc e sem irrigação). 


Médias seguidas de letras iguais não diferem entre si pelo teste Scott-Knott em nível de 5% de 
probabilidade. 
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Analisando-se a Figura 3, é possível observar que comparando os níveis de manejo 
da irrigação, as plantas do nível com 100% da ETc e nível com 50% da ETc, apresentaram 
os maiores percentuais de grãos moca. O clone V12, não apresentou diferenças estatísticas 
quando comparados os manejos, entretanto, quando se analisa o comportamento do clone 
dentro dos níveis de manejo da irrigação, é possível observar que o nível sem irrigação, 
apresentou o maior percentual de grãos moca para este clone. Esses resultados vão de 
acordo com resultados obtidos por Pereira (2015), estudando o efeito de diferentes manejos 
de irrigação sobre o tamanho dos grãos de conilon “Vitória, foi observado que o clone V8, 
apresentou a maior porcentagem de grãos do tipo moca. 

A média obtida da variedade em condições irrigadas para o nível de 100% da 
ETc, 21,20%, em condições irrigadas para o nível de 50% da ETc, apresentou 20,95% 
e em condições para o nível de sem irrigação, a média da variedade foi de 14,75%. Em 
condições irrigadas tanto com o nível de 100% da ETc ou com o nível de 50% da ETc, os 
valores foram próximos ao de 21,40% obtido por Fonseca et al. (2004) como valor médio 
da variedade *Conilon Vitória”. 

O conhecimento desse tipo de grãos dentro é necessário, pois a mistura de grãos tipo 
chato e moca influenciam de forma negativa na qualidade da bebida, pois durante a torra 
do café, os grãos maiores torram lentamente, enquanto os menores torram rapidamente e 
queimam, diminuindo a qualidade da bebida (SILVEIRA et al., 2015). 
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Figura 5 - Percentual de grãos retidos em peneira 13 e superiores para cada manejo de 
irrigação nos clones do cafeeiro Conilon Vitória” 


Médias seguidas de letras iguais não diferem entre si pelo teste de Tukey (A) e Scott-Knott (B) 
em nível de 5% de probabilidade. 


Engenharia Agronômica: Ambientes Agrícolas e seus Campos de Atuação 3 Capítulo 7 





100 





E 9 
q 80 X 

70 a a b b 
8 N q b o CCCc ceuc 
Lu 60 Pe e Cc d 
a E PE d j 
a so el || 

| f 
“is | | | 
30 | | || 
Ê | | ] 
D 2 | | 
É o À | | 
a | | | | 
0 d a 
SEM IRRIGAÇÃO IRRIGADO (50% da E Tc)  IRRIGADO (100% da ETc) 


MANEJOS DE IRRIGAÇÃO 


Vi uV2 uV3 uV4 nVS nV6 mV7 =V8 nV9 nViO 011 nVi2 =V13 


Figura 6 - Percentual de grãos retidos em peneira 13 e superiores em função dos treze clones do cafeeiro 
'Conilon Vitória” para cada manejo (irrigado com 100% da ETc, irrigado 50% da ETc e sem irrigação). 


Médias seguidas de letras iguais não diferem entre si pelo teste de Tukey (A) e Scott-Knott (B) 
em nível de 5% de probabilidade. 


Analisando-se a Figura 5, tem-se o percentual de grãos retidos em peneira 13 e 
superiores. Observa-se que os maiores percentuais foram obtidos em plantas irrigadas para 
o nível de 100% da ETc e plantas irrigadas para o nível de 50% da ETc. Quando comparado 
os níveis do manejo de irrigação, é possível observar que as plantas do nível irrigado de 
100% da ETc, apresentaram maior percentual, não ocorrendo diferença estatística para os 
clones V2 e V8. Quando se estuda o comportamento do clone dentro de cada nível manejo, 
é possível observar que os maiores valores de grãos retidos em peneira 13 e superiores 
foram obtidos no manejo irrigado de 100% da ETc, sendo o V12, o clone que obteve o maior 
percentual do nível de manejo irrigado com 100% da ETc e no nível de manejo irrigado com 
50% da ETc. 

Analisando a Figura 6, é possível observar que no manejo do nível sem irrigação, os 
clones V4 e V12, diferiram estatisticamente dos demais. O clone V7 apresentou o menor 
percentual da variedade em condições sem irrigação. Porém, quando se compara o clone 
V1 e V7 com o nível de 100% da ETc é possível observar que o clone apresentou uma 
melhor resposta. Enquanto o V1 em condições sem irrigadas apresenta 45,20% dos grãos 
retidos em peneira 13 e superiores, para o nível de 100% da ETc, ele apresentou um 
valor de 63,08%. O clone V7, para o nível sem irrigação apresentou um percentual de 
38,55% de grãos retidos em peneira 13 e superiores, enquanto para o nível de 100% da 
ETc apresentou um percentual de 59,49%, mostrando o efeito benéfico da irrigação sobre 
a cultura. Para alguns autores como Dardengo et al (2018), lavouras irrigadas apresentam 
maior percentual de grãos retidos em peneira 13 e superiores, em relação as de sequeiro, 
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isso demonstra maior granação do café irrigado. De acordo Rena e Maestri (2000), 
o tamanho do grão de café é determinado no período compreendido entre a 10º e 17º 
semanas após a florada, quando o fruto expande rapidamente, sendo a água a responsável 
por esse aumento de volume, o que justifica a diferença entre o tamanho dos grãos em 
plantas em condições sem irrigação e irrigadas. 
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Figura 7 - Percentual de grãos retidos em fundagem para cada manejo de irrigação nos clones 
do cafeeiro “Conilon Vitória” 


Médias seguidas de letras iguais não diferem entre si pelo teste de Tukey (A) e Scott-Knott (B) 
em nível de 5% de probabilidade. 
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Figura 8 - Percentual de grãos retidos em fundagem em função dos treze clones do cafeeiro 'Conilon 
Vitória” para cada manejo (irrigado com 100% da ETc, irrigado 50% da ETc e sem irrigação). 


Médias seguidas de letras iguais não diferem entre si pelo teste de Tukey (A) e Scott-Knott (B) 
em nível de 5% de probabilidade. 
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Na Figura 7, observam-se os valores de grãos retidos em fundagem. Comparando- 
se os valores obtidos entre plantas irrigadas para o nível de 100% da ETc, irrigadas para o 
nível de 50% da ETc e para o nível sem irrigação, é possível observar que para o nível sem 
irrigação, apresentou em média o maior valor de fundagem. Quando se compara os clones 
dentro de cada nível do manejo irrigado, observa-se que os clones V1 e V7 apresentaram 
os maiores valores de grãos retidos em fundagem. Para o clone V8 não ocorreu diferença 
estatística entre os níveis de manejos. Quando se analisa os clones em cada nível de 
manejo, é possível observar, que as plantas sem irrigação, apresentaram os maiores valores 
de grãos retidos em fundagem. Os clones V1 e V7, apresentaram os maiores valores retidos 
em fundagem em condições sem irrigadas, porém, quando se compara esses clones com 
os manejos irrigado com nível de 100% da ETc e irrigado com o nível de 50% da ETc, foram 
obtidos valores menores. O clone V1 em condições sem irrigação apresentou 45.65% de 
grãos em fundagem, quando comparado com o nível do manejo irrigado com 50% da ETc, 
apresentou 24,13% e 18.81% quando irrigado com o nível de 100% da ETc. O clone V7 
para o nível sem irrigação apresentou 55,06% de grãos em fundagem, para o nível de 50% 
da ETc, apresenta 40,97% e para o nível de 100% da ETc, 19,07%. O menor valor de grãos 
em fundagem obtidos em plantas irrigadas está relacionado à nutrição das plantas que é 
favorecida pelo fornecimento adequado de água, melhor granação e menor proporção de 
grãos chochos e consequentemente, maior retenção em peneira 13 e superiores. 


41 CONCLUSÕES 


A irrigação influenciou no tamanho dos grãos, plantas irrigadas apresentaram o 
maior valor retido em peneira 13 e superiores. 

Os maiores percentuais de grãos em fundagem podem ser observados em plantas 
não irrigadas. 

Os maiores percentuais de grãos do tipo moca podem ser observados em plantas 
irrigadas com 50% da ETc e sem irrigação. 
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RESUMO: Entre os métodos para estimar 
a Evapotranspiração Potencial (ET;, o de 
Thornthwaite é o menos exigente em relação 
às variáveis de entrada, requerendo, apenas, 
dados de temperatura média do local. Diante 
disso, o presente trabalho teve por objetivo o 
desenvolvimento de um software em linguagem 
Python para viabilizar o planejamento e o manejo 
dairrigação para regiões de difícilacesso a dados 
climáticos, obtidos de estações meteorológicas, 
com apenas a temperatura média mensal como 
dado de entrada, além da informação referente 
a latitude do local. O estudo teve como base os 
dados oriundos de uma estação meteorológica 
do Laboratório de Irrigação e Agrometeorologia 
(LIA), da Universidade Federal de Alagoas 
(UFAL), localizada em Rio Largo, AL. Para 
avaliar o desempenho do software foi realizada a 
comparação entre os dados de ET, obtidos pelos 
métodos de Thornthwaite e Penman-Monteith, 
através dos parâmetros estatísticos: coeficientes 
de determinação (R?), correlação de Pearson (r); 
Erro padrão da Estimativa (EPE) e os índices de 
concordância (d) e confiança (C). A correlação 
entre os valores obtidos apresentou os seguintes 
valores de r, d e c, respectivamente: 0,77; 0,82 
e 0,63. Os resultados permitiram classificar o 
software como “regular” e indicou a viabilidade 
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para seu uso no planejamento e manejo da irrigação. 
PALAVRAS-CHAVE: Água no solo; necessidade hídrica; modelagem; Python. 


DEVELOPMENT OF A PYTHON SOFTWARE TO ESTIMATE THE REFERENCE 
EVAPOTRANSPIRATION THROUGH THE THORNTHWAITE METHOD 


ABSTRACT: Among methods to estimate Potential Evapotranspiration (ETo; Thornthwaite 
methods has minor exigency in relation to input variables, requiring only average temperature 
data. In view of this, the present work aimed to develop software in Python language to enable 
the planning and management of irrigation for regions with difficult access to climatic data, 
obtained from weather stations, with only the average monthly temperature as input, besides 
latitude of the place information. The study was based on data from a meteorological station 
atthe Laboratory of Irrigation and Agrometeorology (LIA), at the Federal University of Alagoas 
(UFAL), located in Rio Largo, AL. To evaluate the performance of the software, a comparison 
was made between the ETO data obtained by Thornthwaite and Penman-Monteith methods, 
using the statistical parameters: coefficients of determination (R2), Pearson correlation (r); 
Standard error of the Estimate (EPE) and the agreement (d) and confidence (C) indices. The 
correlation between the values obtained showed the following r, d and c values, respectively: 
0.77, 0.82 and 0.63. The results allowed to classify the software as “regular” and indicates its 
feasibility for irrigation planning and management using. 

KEYWORDS: Water in the soil; water need; modeling; Python. 


11 INTRODUÇÃO 


Para um planejamento, dimensionamento e manejo de um sistema de irrigação 
eficaz é necessário, a priori, um estudo sobre a água no solo. O conhecimento da demanda 
de água é fundamental uma vez que a necessidade de irrigação é, preliminarmente, 
determinada pelo balanço entre a precipitação e a evapotranspiração da cultura (ET), 
sendo essa, o produto entre a evapotranspiração potencial (ET,) e o coeficiente de cultura 
(kc). (PEREIRA et al., 2007). 

Utilizado por Thornthwaite, o termo “Evapotranspiração” expressa a ocorrência 
simultânea dos processos de perda da água do solo para a atmosfera e da transpiração 
das plantas. Doorenbos e Pruitt (1977) definiram o termo “Evapotranspiração de referência 
(ET,)” como a que ocorre em uma superfície vegetada, com crescimento ativo, cobrindo 
totalmente a superfície do solo e sem deficiência hídrica. 

As estimativas de ET, podem ser calculadas a partir de fórmulas empíricas, 
denominadas de método indireto, baseadas em elementos agrometeorológicos 
disponibilizados através de estações meteorológicas. Dentre as fórmulas empíricas 
utilizadas, a recomendada, conforme a Organização das Nações Unidas para Alimentação 
e Agricultura (FAO), está a do método de Penman-Monteith, utilizada para avaliação 
de outros métodos (ALLEN et al., 1998). Entre os métodos para estimar a ET,, o de 
Thornthwaite é o menos exigente em relação às variáveis de entrada, requerendo, apenas, 
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dados de temperatura média do local. 

Com o avanço tecnológico, linguagens e programas estão sendo utilizados na 
obtenção de resultados, a partir da possibilidade de armazenar e processar dados. Gocic 
& Trajkovic (2010) ao utilizar a linguagem de programação C%, desenvolveram um software 
que apresenta estimativas da evapotranspiração calculadas usando dados meteorológicos 
limitados. 

Dentre as linguagens de programação existentes, Python destaca-se pela facilidade 
de aprendizado e velocidade nos resultados. O presente trabalho teve como objetivo o 
desenvolvimento de um software em Python 3, para disponibilizar informações aos 
produtores e interessados, acerca da ET, em regiões com dados agrometeorológicos 
limitados. Os resultados gerados pelo software são fundamentais ao estabelecimento de 
alternativas de planejamento mais eficazes para utilização de água para a irrigação. 


21 MATERIAL E MÉTODOS 


As variáveis agrometeorológicas utilizadas para calcular as estimativas de ET, e, 
concomitantemente, inseridas no aplicativo, foram obtidas de uma estação meteorológica 
automática pertencente ao Laboratório de Irrigação e Agrometeorologia (LIA), no Campus 
de Engenharias e Ciências Agrárias da Universidade Federal de Alagoas (CECA/UFAL), 
localizada no município de Rio Largo - AL (9'28” S, 35'49” W e 127 metros de altitude), 
situado na região dos Tabuleiros Costeiros de Alagoas. O clima é B,sB',b”, Mesotérmico 
Úmido, conforme classificação de Thornthwaite. A precipitação pluvial anual média é de 
1.800 mm, com o período chuvoso iniciando na primeira quinzena de abril e terminando na 
segunda quinzena de outubro. A temperatura média da região, com uma probabilidade de 
80%, é de 23,14º. Os dados meteorológicos utilizados nesse trabalho correspondem aos 
valores mensais do ano de 2019. 

AET, padrão foi calculada pelo método de Penman-Monteith (FAO56) (Equação 1), 
e denominada de ET,PM seguindo a metodologia proposta por Allen et al. (1998). 


a y 900 U.(e; — eq) 
0,408s (R, — G) + T+ 275 


a s+y(1+034U,) (1) 

Onde ET, — evapotranspiração de referência (mm dia”); R, - saldo de radiação (MJ 
m? dia”); G - fluxo de calor no solo (MJ m? dia"); U, -velocidade do vento a 2 m de altura (m 
s?); es - pressão de saturação do vapor d'água do ar (kPa); ea - pressão do vapor d'água 
do ar (kPa); T - temperatura média diária do ar (ºC); s - inclinação da curva da pressão de 
vapor saturado versus temperatura (kPa ºC"!); y - coeficiente psicrométrico (0,063 kPa ºC"). 

A Rn foi calculada como fração da irradiação solar (Hg) (Rn = 0,572 Hg; com 
coeficiente de determinação (R2) igual a 0,927), feita com medições em grama em 
condições padrão, conforme Allen et al. (1998), realizadas pelo LIA da UFAL. 
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O método de Thornthwaite, empregado no script do programa (ET, MD), foi aplicado de 
acordo com Pereira et al. (2007): 

ETpT = 16(107)º 0<Tn<26,0ºC (2) 

ETpT = —415,85+32,224Tn — 0,43 Tn? Tn > 26,0 (3) 


L= TE (02 TP (4) 


a=6,75107[ — 7,7110º [2 + 1,7912 1021 + 0,49239 (5) 


Onde ETpT — evapotranspiração potencial (mm mês"!); Tn — temperatura média do 
mês (ºC); /— índice que expressão nível de calor na região; a — expoente função do índice. 
Para obter-se a ET, do mês correspondente, o valor de ETpT é corrigido em função 


do fotoperíodo e do número de dias do mês, portanto: 
ETot = ETpT Cor (6) 
Cor = (ND/30)(N/12) (7) 
Onde ND — número de dias do mês em questão; N — fotoperíodo daquele mês. 
Para entrada de dados, foi desenvolvida uma interface interativa para o software, 


visando facilitar a manipulação dos dados e agilizando a interpretação dos resultados, 


conforme Figura 1. 


- 
£ Exvapotranspiração — a >< 


EVAPOTRANSPIRAÇÃO 
Estado Lotituder 





[em 








Fotoperíodo 


Fevereiro : 
Março : 
Abril - 
Maio - 
Junho : 
Julho : 
Agosto : 
Setembro :- 
Outubro : 
Novembro : 


Dezembro : 








Figura 1 - Interface gráfica do software para entrada de dados 
Fonte: Os Autores (2021) 
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Os resultados da estimativa de ET, aparecem, para o usuário, de forma intuitiva 
e de simples visualização e manipulação. Foi realizada a análise dos valores estimados 
pelo método padrão (ET,PM) e os estimados através dos cálculos do software a fim de 
se obter a correlação entre eles, baseando-se nos seguintes indicadores estatísticos: 
coeficiente linear (a); coeficiente angular (b); coeficiente de correlação (r); coeficiente de 
determinação (R2?); índice de concordância (d) de Willmott; erro padrão da estimativa (EPE) 
e pelo coeficiente de confiança ou desempenho (c). O erro padrão da estimativa (EPE) (mm 
mes”), descrito por Allen et al. (1986), foi determinado conforme Equação 8: 


— |[XOi-P)? 
EPE = CE (8) 


O índice de Willmott (d), proposto por Willmott et al. (1985), foi determinado conforme 
Equação 9: 

X(P;-09)2 

Ep rea (9) 

Onde Pi — ETPt (mm mês”); Oi — EtoPM (mm mês”); O — média dos valores de 
EToPM (mm mês"). 

Os valores para “d” podem variar de 0, para nenhuma concordância, a 1, para maior 
concordância. O índice “c”, proposto por Camargo & Sentelhas (1997), é obtido através do 
produto entre o índice de Willmott (d) e o coeficiente de correlação (r). A interpretação do 
resultado deve estar de acordo com os referidos autores, conforme a Tabela 1: 





Índice de desempenho (c) Classificação 
>0,85 Ótimo 
0,76 — 0,85 Muito bom 
0,66 — 0,75 Bom 
0,61 — 0,65 Mediano 
0,51 — 0,60 Regular 
0,41 — 0,50 Mal 
<40 Péssimo 





Tabela 1 - Critério de interpretação do desempenho do método pelo índice “c” 
Fonte: Camargo & Sentelhas (1997). 


31 RESULTADOS E DISCUSSÃO 


A Tabela 2 apresenta os parâmetros estatísticos e os resultados do índice 
classificatório de desempenho da análise feita a partir da comparação dos resultados 
da ET,, calculados pelo modelo padrão (ET,PM), com as variáveis obtidas através da 
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estação meteorológica, e os resultados obtidos através do software, utilizando o modelo de 
Thornthwaite. Os resultados da estimativa da ET, são apresentados pelo software através 
de sua interface gráfica (Figura 2). 


£ Evapotranspiração — a > 


EVAPOTRANSPIRAÇÃO 





Localização 


Cidade Estado 
Rio Largo A -2474861111 





Meses Fotoperíodo 





Janeiro : 12.48 
Fevereiro : az.s 

Março : 12.05 
Abril : E. 
Maio : 21.586 
Junho : EE o 
Julho : EE 
Agosto :- E; 

Setembro : 11.95 
Outubro : 12.22 
Novembro 12.44 
Dezembro 12.54 














Figura 2 - Interface gráfica do software, com as entradas, nos widgets destacados em amarelo 
e as saídas (fotoperíodo médio do mês e ET) 


Fonte: Os Autores (2021) 


Ao analisar os coeficientes linear (a) e angular (b) do método utilizado pelo software, 
observa-se que foram estatisticamente significativos (P>0,01, teste t). O método utilizado 
pelo software (ET,T) apresentou apenas um desempenho “regular”, conforme Camargo & 
Sentelhas (1977b), porém, obteve boa correlação entre os resultados. Conceição & Mandelli 
(2005) avaliando o desempenho dos métodos que utilizam a temperatura do ar como 
variável de entrada e comparando os resultados aos valores obtidos através do método 
padrão-FAO, observaram que os valores de (d) e (c) para o método de Thornthwaite foram, 
respectivamente, 0,79 e 0,58. Esses valores são ligeiramente inferiores ao observados 
neste trabalho, essa divergência ocorre em virtude das diferenças climáticas entre as 
regiões estudadas. Camargo et al. (1999) constataram que, para regiões de clima seco, 


usando a temperatura média convencional (Tmed) como variável de entrada, os valores da 
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ET, foram superestimados. Porém, ao utilizar a temperatura média efetiva (Tef) observou- 


se melhora no desempenho da equação. 








Método R2 EPE a b d c Desempenho*** 
Elk 0,60 14,18 0,8381** 10,68** 0,82 0,63 Regular 
(software) 





Tabela 2 - Índice de Concordância de Willmott (d) e desempenho (c) para valores mensais 
de ET, calculada pelo software e comparados com o método padrão de Penman-Monteith. 
Período de janeiro a dezembro de 2019, Rio Largo, AL (WILLMOTT et al.,1985). 


* Significativo a nível de 1% de probabilidade pelo teste t; *** Classificação do desempenho, conforme 
Indice c, de acordo com Camargo & Sentelhas (1997a). 


Fonte: Os Autores (2021) 


A correlação linear entre as variáveis, ET,T e ET,PM (Figura 3) é classificada, de 
acordo com Ferreira (2018) como correlação positiva alta, com um coeficiente de correlação 
de Pearson igual a 0,78. Gocic & Trajkovic (2010), ao analisar os resultados, encontraram 
uma concordância muito boa, em relação ao método padrão (ET,PM), utilizando apenas a 
temperatura do ar como variável de entrada. O que mostra concordância com os resultados 


deste trabalho. 





50 À 


40 4 


EToPM (mm mês 1 


30 À 


20 4 








00 + 


90 À 











80 90 100 110 120 130 140 


ETpT (mm mês |) 


Figura 3 - Análise de regressão linear para Evapotranspiração Potencial pelo método de 
Thornthwite (ETOT) e Evapotranspiração de Referência, calculada pelo método padrão, 
Penman-Monteith (FAO56) (ETOPM), do ano de 2019 em Rio Largo, Alagoas. 


Fonte: Os Autores (2021) 
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41 CONCLUSÕES 


Para regiões onde as variáveis meteorológicas são limitadas, o uso do software 
através do método de Thornthwaite, mostrou-se viável como uma ferramenta eficaz para 
uso no planejamento e manejo da irrigação. 
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RESUMO: A determinação do real consumo 
hídrico das culturas agrícolas é fundamental para 
um manejo eficiente da irrigação, nesse sentido, a 
determinação da evapotranspiração das culturas 
torna-se peça chave. Dessa forma, objetivou-se 
a construção e calibração de quatro lisímetros de 
pesagem no interior de casa de vegetação. Os 
lisímetros foram constituídos por um vaso de 250 
L, com área de 0,75 m2. O sistema de pesagem 
dos lisímetros foi construído utilizando células de 
carga da marca HBMO, modelo Z6FC3 200kg, 
sendo utilizadas três células de carga, em 
disposição triangular, por lisímetro. A calibração 
foi realizada com 29 massas padrão que foram 
assumidas como referência, sendo sete de 0,47 
kg (0,15 mm), sete de 1,57 kg (0,50 mm), oito 
de 6,28 kg (2,00 mm) e sete de 15,7 kg (5,00 
mm). Os coeficientes angulares encontrados nos 
diferentes lisímetros apresentaram valores de 
305,02; 305,42; 306,35 e 305,85 (kg:(mV-V-1)-1), 
para os lisímetros 1, 2, 3 e 4, respectivamente. 
Verificou-se que os lisímetros apresentaram altos 
valores de coeficiente de determinação após a 
calibração e baixos erros foram observados, 
confirmando que o uso de lisímetros de pesagem 
em casa de vegetação é uma boa opção para 
determinação da evapotranspiração de cultura. 

PALAVRAS-CHAVE: Evapotranspiração, 
agricultura irrigada, conteúdo de água no solo. 
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CONSTRUCTION, CALIBRATION AND PERFORMANCE OF WEIGHING 
LYSIMETERS TO DETERMINE THE EVAPOTRANSPIRATION OF CULTURE 


ABSTRACT: The determination of the real water consumption of crops is essential for 
efficient management of irrigation, in this sense, the determination of the evapotranspiration 
of the crops becomes a key part. Thus, the aim was to build and calibrate four weighing 
lysimeters inside the greenhouse. The lysimeters were constituted by a 250 L vessel, with an 
area of 0.75 m2. The weighing system of the lysimeters was built using load cells of the brand 
HBMO, model Z6FC3 200kg, using three load cells, in a triangular arrangement, per Iysimeter. 
Calibration was performed with 29 standard masses that were taken as reference, seven of 
0.47 kg (0.15 mm), seven of 1.57 kg (0.50 mm), eight of 6.28 kg (2, Mm) and seven of 15.7 kg 
(5.00 mm). The angular coefficients found in the different lysimeters showed values of 305.02; 
305.42; 306.35 and 305.85 (kg:(mV-V-1)-1), for Iysimeters 1, 2, 3 and 4, respectively. It was 
found that the lysimeters showed high values of determination coefficient after calibration and 
low errors were observed, confirming that the use of weighing lysimeters in the greenhouse is 
a good option for determining the culture evapotranspiration. 

KEYWORDS: Evapotranspiration, irrigated agriculture, soil water content. 


11 INTRODUÇÃO 


A determinação da evapotranspiração (ET) é fundamental para o correto manejo 
da irrigação, sendo que a demanda hídrica de uma cultura pode ser obtida com maior 
precisão utilizando-se lisímetros. As principais vantagens do seu uso estão relacionadas à 
calibração, que pode ser validada mecanicamente; ao requerimento de bordadura, que é 
menor quando comparada a outros métodos, como micrometeorológicos, por exemplo; a 
possibilidade de completa automação e a facilidade de inspeções, pois o sistema é quase 
totalmente visível (Allen et al., 2011). 

Na literatura são descritos diferentes tipos de lisímetros, sendo que o tamanho e 
o projeto de montagem de cada um são baseados no requerimento específico do estudo 
e das condições de construção e custos envolvidos (Faria, Campeche e Chibana, 2006). 
Os lisímetros de pesagem, cujo sistema de medição é baseado em células de carga são 
considerados os de maior confiabilidade nas medições (Cruz-Blanco et al., 2014), sendo 
utilizados para calibração dos demais métodos de estimativa de evapotranspiração. 

Talvez o maior fator limitante para o emprego de lisímetros, principalmente os de 
pesagem, são os custos relativamente altos envolvidos em sua construção (Flumignan, 
2012), o que fazem com que esses equipamentos sejam utilizados exclusivamente por 
instituições de ensino e/ou pesquisa. No entanto, a utilização desses equipamentos irão 
fornecer subsídios para que o manejo da irrigação seja realizado de forma o mais precisa 
possível, uma vez que como resultado de saída de um lisímetro bem montado e construído 
são os valores de coeficiente de cultivo (Allen et al., 1998; Allen e Pereira, 2009; Taylor et 
al., 2014). 

É importante destacar que o coeficiente de cultivo de uma determinada cultura diz 
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respeito às características fenológicas e fisiológicas da cultura em questão (Marin, Jones 
e Boote, 2017), sendo, portanto o elo de ligação de conceitos físicos, uma vez que o 
Kc é determinado por meio da evapotranspiração da cultura e da evapotranspiração de 
referência, com a biologia da planta. A grande limitação em se fazer o manejo da irrigação 
por meio da determinação de Kc está no fato de que o mesmo tem que ser um valor que 
represente as condições locais de cultivo, sendo o mesmo totalmente dependente das 
condições meteorológicas decorrentes da estimativa da evapotranspiração de referência. 

Essa afirmativa é verdadeira quando (Doorenbos e Pruitt, 1975)afirmaram que 
existe, obrigatoriamente, a necessidade de calibração dos coeficientes de cultivo para 
as condições locais, além das determinações serem realizadas ao longo de todos os 
estádios de desenvolvimento das culturas. Outros pesquisadores já confirmaram essa 
afirmação (Allen et al., 2011; Payero e Irmak, 2013; Pereira et al., 20083; Silva, Da et al., 
2016; Taylor et al., 2014), e quando se pensa em diminuição de riscos nas safras agrícolas 
principalmente no que diz respeito a culturas que não necessitam de grandes áreas de 
cultivos, a construção e avaliação de lisímetros para serem utilizados no interior de casas 
de vegetação se tornam uma realidade. Entender o comportamento das diferentes espécies 
vegetais, quando submetidas às condições micrometeorológicas no interior de casas de 
vegetação, se mostra fundamental para essa nova fase da agricultura que podemos viver 
futuramente. 

Dessa forma, o objetivo deste trabalho foi realizar a construção, calibração e 
avaliação de quatro lisímetros de pesagem operando no interior de uma casa de vegetação. 


21 MATERIAL E MÉTODOS 


2.1 Caracterização da área experimental e do material solo 


O experimento estudo conduzido na área experimental do Departamento de 
Engenharia de Biossistemas da Escola Superior de Agricultura “Luiz de Queiroz” (ESALQ/ 
USP), em ambiente protegido (casa de vegetação), localizada no município de Piracicaba 
- SP (latitude 22º43'31”; longitude 47º38'57”; altitude 547 m). 

O solo utilizado para preenchimento dos lisímetros foi classificado como sendo um 
Latossolo Vermelho-amarelo distrófico (Embrapa, 2013), obtido na profundidade de 0,30 m 
da camada agricultável, proveniente do campus da ESALQ/USP. As propriedades físicas 
do solo foram determinadas realizando-se análise granulométrica, conforme metodologia 
proposta por (BOUYOUCOS, 1951): densidade das partículas, utilizando o método do 
picnômetro e densidade do solo, determinada pelo método da proveta, os quais são 
recomendados pela (Embrapa, 1997) (Tabela 1). 
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Frações granulométricas 











Camada cc PMP CAD Ds Dp PT 
(cm) (mm) (%) Areia Silte Argila 
cmê cm g cm? % 
0-20 0,226 0,104 24,4 1,48 2,65 44,7 75,1 7,8 17,1 


20-40 0,234 0,109 25,0 1,49 2,65 46,9 74,5 8,6 17,0 





CC — Capacidade de campo; PMP — ponto de murcha permanente; CAD — capacidade de água 
disponível; Ds — Densidade do solo; Dp — Densidade das partículas; PT — Porosidade Total. 


Tabela 1. Caracterização físico-hídrica do solo utilizado para preenchimentos dos lisímetros 


O controle de umidade do solo foi realizado por meio da utilização de tensiômetros 
instalados nas profundidades de 0,10, 0,20 e 0,30 m, com as leituras de tensão sendo 
realizadas por meio de transdutores de pressão, devidamente calibrados para as condições 
de operação. A partir dos valores de tensão obtidos, a umidade foi estimada por meio da 
equação de Van Genuchten (Genuchten, Van, 1980), utilizando os parâmetros de ajustes 
apresentados na . Os parâmetros de ajuste do modelo foram determinados por meio do 
software RETC, desenvolvido pelo US Salinity Laboratory — USDA (Genuchten, Van, Leij e 
Yates, 1991), utilizando os pares de dados obtidos na curva de retenção de água no solo. 








Camada (cm) 0s(m'm?) 9r (mº mº) a (kPa?) n m 
0-20 0,4736 0,0928 0,0817 1,4918 0,329669 
20-40 0,4737 0,0998 0,0657 1,5326 0,347514 





Tabela 2. Valores de umidade de saturação (0s) e residual (Br) e dos parâmetros empíricos (a, 
ne m) do modelo de Van Genuchten, utilizados para estimativa da umidade do solo. 


2.2 Construção dos lisímetros 


Foram construídos quatro lisímetros a partir de vasos com capacidade para 250 
L, área de 0,75 m?, com sistema de pesagem direta, utilizando-se células de carga que 
suportavam diretamente todo o peso do sistema. Estes foram preenchidos com camada de 
pedra brita (em torno de 8 cm), seguida de manta geotêxtil e, por fim, o volume restante foi 
completado com solo descrito no item 2.1. Dessa forma, considerando que todo o sistema 
apresentava peso aproximado de 550 kg (vaso preenchido com solo, planta, água, sistema 
de apoio e drenagem), cada lisímetro foi composto por três células de carga com capacidade 
nominal para 200 kg, as quais quando combinadas resultaram num total de 600 kg. 

O sistema de pesagem dos lisímetros foi construído utilizando células de carga 
da marca HBMO, modelo Z6FC3 200kg. Conforme especificações da célula de carga, a 
mesma suporta uma carga máxima de 300 kg sem que haja ruptura em seu mecanismo 
de pesagem. Construiu-se uma estrutura triangular, fixando-se uma célula de carga em 
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cada vértice para compor o mecanismo de pesagem. A estrutura foi colocada sobre piso de 
concreto tomando-se o devido cuidado para que a mesma permitisse que as três células 
ficassem em nível. Posteriormente, sobre as células foi colocado uma base que serviu de 
suporte para o acomodamento do vaso. 

Para obtenção do sinal de saída das células de carga e a metodologia de ligação 
no sistema de aquisição de dados, resulta na necessidade de se ligar cada célula de carga 
em um canal diferencial e uma porta de excitação para conexão com o sistema automático 
de aquisição de dados, dessa forma optou-se por fazer uma caixa de junção, unindo-se as 
três células de carga, assim foi possível utilizar apenas um canal diferencial por lisímetro. 


2.3 Calibração dos lisímetros 


Com base na análise dos dados de solo, e considerando o solo totalmente seco, 
seriam necessários cerca de 165 L de água para elevar a umidade do solo para além da 
capacidade de campo, ou seja, preencher os espaços porosos. Desta forma, utilizou-se 
175 kg para representar a CAD do solo durante a calibração. Os 175 kg foram divididos 
em 29 massas padrão que foram assumidas como referência, sendo sete de 0,47 kg (0,15 
mm), sete de 1,57 kg (0,50 mm), oito de 6,28 kg (2,00 mm) e sete de 15,7 kg (5,00 mm). 
Os vasos sobre o lisímetro tiveram sua superfície coberta com lona plástica para evitar a 
evaporação da água do solo durante o procedimento de calibração. Sabe-se que ventos 
laterais podem resultar em erros de leituras das células de carga, por exercerem uma 
força vetorial longitudinal ao lisímetro. No entanto, como os mesmos foram construídos e 
instalados no interior de casa de vegetação, desprezou-se os efeitos dos ventos laterais 
na calibração e operação. A disposição das massas durante a calibração foi realizada da 
maneira mais distribuída possível pela área dos lisímetros, procurando-se não as concentrar 
somente em um ponto. 

Asequência de adição e retirada das massas nos lisímetros foi realizada intercalando- 
as, sendo adicionadas as massas de nº1 a 29, visando explorar melhor as características 
das células em cada faixa de utilização das mesmas. Posteriormente estas foram retiradas 
no sentido inverso. A calibração foi realizada em um lisímetro por vez. A calibração foi 
realizada em um lisímetro por vez. Na programação feita para o sistema de aquisição 
de doados foi adicionado uma chave para ligar e desligar a leitura e o armazenamento 
dos dados. Desta forma, sempre antes de se adicionar ou retirar uma massa, a coleta e 
armazenamento eram desligadas, para evitar a interferência das oscilações à adição ou 
retirada das massas. 

Após coleta dos dados de calibração, os pares formados pela média do sinal das 
duas células de carga e a massa acumulada no lisímetro, foram utilizados para ajustar um 
modelo de regressão da calibração. Para a avaliação do desempenho dos modelos dos 
lisímetros, obtidos na calibração, na comparação das estimativas, adotou-se os indicadores 
estatísticos índice de concordância de Willmott “d” (WILLMOTT et al., 1985), índice de 
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concordância e desempenho “c” (CAMARGO; SENTELHAS, 1997), os coeficientes de 
correlação de Person (r) e de determinação (R2) entre os sinais de saída das células de 
carga (Yi) e os somatórios dos valores das massas padrão (Yi), o erro absoluto (EA), o erro 
médio absoluto (ERA) e a raiz do quadrado médio do erro (RQME). 


2.4 Evapotranspiração e coeficiente de cultivo 


No interior de cada lisímetro foi cultivado a cultura do pinhão manso, sendo 
que a evapotranspiração da cultura (ETc) foi determinada por meio da utilização dos 
lisímetros de pesagem. Os valores de evapotranspiração foram determinados por meio 
da contabilização das entradas e saídas de água no sistema lisimétrico. De posse dos 
dados de evapotranspiração, os valores de coeficiente de cultivo (Kc) foram determinados 
de acordo com apresentado em (Allen et a!., 1998). Sendo assim, os valores horários de 
ETo foram calculados no programa REF-ET (Pereira et al., 2015)e acumulados para o total 
diário. Para isso foram utilizados os dados obtidos do monitoramento meteorológico no 
interior da casa de vegetação realizado por meio da instalação, no centro do experimento, 
de sensores de radiação global modelo CM3 Kipp & ZonenQ); temperatura e umidade 
relativa do ar — termohigrômetro modelo HMP45C da Vaisala?), com os dados coletados a 
cada segundo e média a cada quinze minutos (Figura 1B). Para coleta e armazenamento 
dos dados, utilizou-se um sistema automático de aquisição de dados. Os sensores foram 
instalados a 2 m de altura. A estimativa da evapotranspiração de referência foi realizada 
pelo método Penman-Monteith, conforme a padronização disposta em (Snyder e O'Connell, 
2007). 


31 RESULTADOS E DISCUSSÃO 


3.1 Calibração dos lisímetros 


Durante a calibração dos lisímetros, realizado à noite, a temperatura do ar variou 
entre 22,9 e 24,2ºC e a umidade relativa entre 66,2 e 72,3%. Por se encontrar em casa de 
vegetação a velocidade do vento foi considerada nula. 

Os coeficientes angulares encontrados nos diferentes lisímetros apresentaram 
valores de 305,02; 305,42; 306,35 e 305,85 (kg (mV-V'')), para os lisímetros 1,2,3 e 4, 
respectivamente. A comparação de similaridade dos coeficientes angulares entre pares de 
retas não mostrou diferenças estatisticamente significativas pelo teste t (p < 0,05). 

Considerando-se os altos coeficientes de determinação obtidos (Figura 1), o 
modelo de regressão linear foi ajustado para relação entre a massa acumulada e a média 
do sinal de saída das células de carga, evidenciando o comportamento linear existente 
entre estas variáveis. Entretanto, embora tenha sido comprovada a linearidade, isso não 
necessariamente indica que os mesmos sejam exatos ou acurados. Esses conceitos, os 


quais são sinônimos, são diferentes do de linearidade, pois representam a capacidade que 
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o modelo ajustado tem de determinar o valor verdadeiro da massa. Assim, alguns erros 
ocorreram entre a massa determinada por estes e a acumulada sobre os mesmos durante 
a calibração, implicando em diferentes acurácias entre os lisímetros. Estas acurácias são 


representadas pelo erro padrão da estimativa do modelo de regressão ajustado. 
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Figura 1. Equações de ajustes dos dados de calibração dos lisímetros de pesagem por meio da 
adição e retirada das massas padrão 


Na calibração de lisímetros sempre espera-se encontrar equações lineares com 
elevados coeficientes de determinação, uma vez que estas irão resultar em um menor 
efeito da histerese, equações dessa natureza também foram observadas por (Amaral et 
al., 2018; Carvalho et al., 2006; Flumignan, 2012; Silva, Da et al., 2016). A histerese diz 
respeito ao erro que se obtém entre o sinal que é obtido pelo sistema de aquisição de 
dados para uma mesma quantidade de massa, sempre considerando a adição e a retirada 
da mesma massa. 

A presença da histerese, em alguns pontos da faixa de calibração é inevitável e até 
aceitável, porém há a necessidade de se realizar a avaliação dos erros, para obtenção da 
real capacidade do lisímetro em fornecer valores de massa que represente realmente a 
variação de umidade de água no solo. Para este trabalho, verificou-se que os erros foram 
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pouco dispersos em torno do zero, sendo observada variação de -0,21 até 0,23 kg, para 
o lisímetro 1. Valores de erro variando entre -0,42 a 0,57 kg foram observados para o 
lisímetro 2, para o lisímetro 3 os erros foram de -0,46 a 0,42 kg e valores de erro entre 0,39 
a ,48 kg foram verificados no lisímetro 4. As dispersões dos erros em torno do zero podem 


ser observadas na Figura 2. 
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Figura 2. Erros entre a massa acumulada e a massa determinada durante a calibração dos 
lisímetros. 


Os erros encontrados neste trabalho podem ser considerados baixos, o que 
pode ser comprovado pela análise da Tabela 2 onde estão apresentados o coeficiente 
de determinação e correlação linear de Pearson, além dos índices de concordância de 
Willmott e de concordância de desempenho de (Camargo e Sentelhas, 1997), bem como o 
erro absoluto (EA), erro médio absoluto (ERA) e raiz do quadrado médio do erro (RQME). 

Os indicadores de desempenho apresentados na Tabela 2 confirmam que o modelo 
linear apresentou os melhores ajustes, resultando em baixos valores de EA, ERA e RQME. 
De acordo com (Janssen e Heiberger, 1995), quanto mais o valor de R2 e r se aproximarem 
de 1,0 e quanto mais os valores de erro (EA, ERA e RQME) se aproximaram de zero, maior 
é a precisão do modelo ajustado. 








Lisímetros R2 r d c EA (%) ERA(%)  RQME (%) 
1 1,0000 1,0000 0,9970 0,9970 1,6300 0,0016 0,0021 
2 1,0000 1,0000 0,9967 0,9966 2,2600 0,0020 0,0026 
3 1,0000 1,0000 0,9962 0,9961 3,9800 0,0026 0,0032 
4 1,0000 1,0000 0,9967 0,9967 2,7900 0,0021 0,0026 





Tabela 3. Valores de coeficiente de determinação (R?), coeficiente de correlação linear de 
Pearson (r), índice de concordância de Willmott (d), índice de concordância e desempenho 
(c), erro absoluto (EA), erro médio absoluto (ERA) e raiz do quadrado médio do erro (RQME) 

obtidos por meio da análise dos dados de calibração dos quatro lisímetros de pesagem 
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3.2 Evapotranspiração e coeficiente de cultivo 


Na Figura 3 estão apresentados os valores de evapotranspiração de referência e 
evapotranspiração da cultura do pinhão manso, obtidos pelo método de Penman-Monteith 
e por lisímetros, respectivamente, obtidos ao longo de 350 dias de avaliação. 


+Eto P-M (mm) e ETC Lisímetro 


q 


Evapotranspitação (mm dia”) 
+ 





0 50 100 150 200 250 300 350 
Tempo de avaliação (dias) 


Figura 3. Evapotranspiração da cultura e evapotranspiração de referência, ambas em mm dia- 
1, com os lisímetros em funcionamento por um período de 350 dias. 


A ETc e ETo máxima calculada para o período foram de 6,96 e 
7,54 mm-d”, respectivamente. Vale ressaltar que o pinhão-manso é uma planta caducifólia, 
ou seja, ocorre a queda das folhas no inverno, com recuperação das mesmas a partir do mês 
de setembro. Durante essa fase a evapotranspiração da cultura é praticamente governada 
pela evaporação do solo. Observando atentamente Figura 3 percebe-se uma diminuição da 
ETc, a partir do dia 180 após, que se caracteriza pelo início da senescência foliar do pinhão- 
manso. Após as plantas começarem a formar novos ramos e folhas (período vegetativo), 
os valores atingiram um patamar variando entre 4 mm-d'! até aproximadamente 8 mm-d”. 

Os valores de coeficiente de cultivo (kc) médio de 0,75 (Figura 4). Porém, nota-se 
que no dia 157, foi observado um valor de Kc de 2,37, o que pode representar um erro na 
medida desse valor, quer seja por dados inconsistentes obtidos na estação meteorológica 
ou algum erro na obtenção do sinal das células de carga. O que se percebe, é que valores 
de Kc acima de 2,0 é muito raro de se observar, independente da cultura que se está 
avaliando (Carvalho et al., 2016; Ruiz-Pefialver et al., 2015; Silva, Da et al., 2016; Taylor 
et al., 2014). 
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Figura 4. Coeficiente de cultivo (Kc, adimensional) obtido ao longo do período de avaliação dos 
lisímetros. 


Rajaona; Sutterer; Asch (2012)encontraram valores de ETc de 1,9 mm:d” nos 
períodos inicial e final e 5,5 mm:d' no período médio utilizando valores de Kc adaptados 
conforme dados encontrados na literatura. Os valores de Kc encontrados foram de 0,6, 
1,2 e 0,4 para os estádios inicial (vegetativo), médio (reprodutivo) e final (senescência), 
respectivamente. Vale ressaltar que os valores de evapotranspiração são totalmente 
dependes das condições atmosféricas locais, portanto os períodos de avaliação influenciam 
nos resultados, o que torna importante a determinação desses valores in loco. Diante disso, 
o máximo valor de ETo obtido por Rajaona; Sutterer; Asch (2012)foi da ordem de 4,7 mm:d”, 
ao passo que neste trabalho o máximo valor foi de 7,54 mm-d”, indicando que as condições 
de cultivo das plantas de pinhão manso se apresentava diante de uma atmosfera com um 


potencial muito mais turbulento. 


41 CONCLUSÕES 


Nas condições em que este estudo foi realizado, pode-se concluir que, os lisímetros 
de pesagem construídos resultaram em equações de ajustadas que apresentaram alto 
coeficiente de determinação e correlação de Pearson, com baixos valores de erros 
verificados, tendo, portanto, se ajustado adequadamente aos dados amostrais. Os baixos 
erros e o período de avaliação, demonstra que os lisímetros podem ser utilizados, de forma 
confiável, para a obtenção da evapotranspiração de cultura e, consequentemente, de 
coeficiente de cultivo em culturas cultivadas no interior de casas de vegetação. 
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RESUMEN: En los últimos afos se ha 
intensificado el uso agrícola del suelo y con 
ello la práctica de fertilización de los cultivos 
de granos. Con el objetivo de evaluar el efecto 
de los distintos regímenes de fertilización 
de largo plazo sobre: (a) los contenidos de 
materia orgânica de suelos (MOS) y nitrógeno 
potencialmente mineralizable (Nan) en el suelo 
y (b) la relación entre Nan y MOS, se llevó a 
cabo la siguiente experiencia. Se analizaron las 
concentraciones de MOS y Nan en suelo de tres 
estratos de profundidad (0-5 cm, 5-10 cm y 10-20 
cm) de cinco sitios de la Red de Nutrición CREA- 
IPNI de la Región Sur de Santa Fe instalados 
en la campafia 2000-2001. Estos ensayos se 
instalaron en siembra directa (SD) bajo dos 
sistemas de rotación, maíz-trigo/soja (M-T/S) 
(Balducchi y San Alfredo) y maíiz-soja-trigo/soja 
(M-S-T/S) (La Blanca, La Hansa y Lambaré), 
seis tratamientos de fertilización: testigo, PS, NS, 
NP, NPS y completo (NPS+micronutrientes). La 
fertilización con N, P y S provocó aumento en 
los contenidos de MOS (sólo en el estrato de 
0-5 cm de profundidad) pero no afectaron a los 
contenidos de Nan en el suelo. Las rotaciones 
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no afectaron a los contenidos de MOS y Nan. Se 
determinaron relaciones positivas entre el Nan 
y la MOS y sus variaciones fueron atribuibles a 
las diferencias en otras características además 
del contenido de MOS total entre los sitios. Se 
concluye que solo la MOS total no es un buen 
predictor de la capacidad de mineralización de N 
de los suelos. 

PALABRAS CLAVE: Fertilización de cultivos — 
relación Nan: MOS - capacidad de mineralización 
de suelos. 


BALANCE OF ORGANIC MATTER 
AND NITROGEN MINERALIZATION 
CAPACITY OF DIFERENT SOILS WITH 
CONTINUOUS FERTILIZATION 
ABSTRACT: In recent years, the agricultural use 
of the soil has intensified and with it the practice 
of fertilizing grain crops. With the objective 
of evaluating the effect of the different long- 
term fertilization regimes on: (a) the contents 
of soil organic matter (SOM) and potentially 
mineralizable nitrogen (Nan) in the soil and, 
(b) the relationship between Nan and MOS, 
the following experience was carried out. The 
concentrations of MOS and Nan in soil from three 
depth strata (0-5 cm, 5-10 cm and 10-20 cm) of 
five sites of the CREA-IPNI Nutrition Network 
of the South Region of Santa Fe installed in the 
2000-2001 campaign. These trials were installed 
in direct sowing (SD) under two rotation systems, 
corn-wheat / soybean (MT / S) (Balducchi and 
San Alfredo) and corn-soybean-wheat / soybean 
(MST / S) (La Blanca, La Hansa and Lambaré), 
six fertilization treatments: control, PS, NS, NP, 
NPS and complete (NPS + micronutrients). The 
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fertilization with N, P and S caused an increase in the SOM contents (only in the 0-5 cm deep 
stratum) but did not affect the Nan contents in the soil. The rotations did not affect the contents 
of MOS and Nan. Positive relationships between Nan and MOS were determined and their 
variations were attributable to differences in other characteristics in addition to total MOS 
content between sites. It is concluded that only the total SOM is not a good predictor of the 
mineralization capacity of N of the soils. 

KEYWORDS: Crop fertilization - Nan: MOS ratio - soil mineralization capacity. 


INTRODUCCIÓN 


La producción mundial de alimentos debe aumentar en un 70% en 2050 para seguir el 
ritmo de la creciente demanda (Bruinsma, 2009). Para alcanzar este objetivo, la fertilización 
surge como una herramienta para aliviar el efecto de la continua exportación de nutrientes, 
debido a su capacidad para reponer la fertilidad del suelo y mejorar la producción agrícola 
(Wu & Ma, 2015). Esa necesidad de producir más alimentos y el deterioro de la fertilidad 
natural de los suelos ha impulsado a la fertilización en forma continua por parte de los 
productores. En la agricultura Argentina también se ha intensificado el uso de fertilizantes 
para cultivos extensivos. El consumo de fertilizantes tuvo un gran aumento en los últimos 
20 afios, pasando de 300 mil toneladas en 1990 a más de 3 millones de toneladas en el 
2013 (Fertilizar, 2014a). 

La fertilización continua pareciera ser eficaz en el mantenimiento de la materia 
orgánica del suelo (MOS) (Wu et al., 2003) y enriquecimiento con carbono (C) principalmente 
de las fracciones lábiles (Yang et al., 2012). En el suelo se hallan dos fracciones de C con 
distinta dinâmica de descomposición. Una fracción lenta o estable y otra fracción activa 
o lábil (Ciampitti et al., 2011b). A su vez, estos pooles actúan como fuente primaria de 
nitrógeno (N) para las plantas (Haynes, 2005). Las plantas toman nutrientes de la fracción 
lábil, la cual se renueva constantemente a partir de la fracción estable. Una forma de evaluar 
la fracción lábil de N en el suelo es midiendo el amonio liberado durante una incubación 
anaeróbica corta (Nan). Varios estudios mostraron que el Nan estuvo relacionado con la 
capacidad del suelo de proveer N a cultivos de granos (Sharifi et al., 2008; Genovese, 
2009; Reussi Calvo et al., 2014). La posibilidad de contar con experimentos a largo plazo de 
fertilización con secuencias de cultivos en distintos sitios, pueden proveer de información 
valiosa para conocer el impacto del agregado continuo de nutrientes como fertilizantes 
sobre las concentraciones de la MOS y Nan. 

En la campafia 2000-2001, la Región Sur de Santa Fe del movimiento CREA, en 
conjunto con IPNI y Agroservicios Pampeanos (ASP), instalaron una red de ensayos a largo 
plazo. Uno de sus objetivos comunes desde un principio fue determinar las respuestas 
directas y residuales de los cultivos a la aplicación de N, P y azufre (S) y así complementar 
a los métodos de diagnóstico de la fertilidad de suelos. En este trabajo se sintetizan los 
resultados obtenidos de análisis de suelos de tres estratos de profundidad (0-5 cm, 5-10 
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cm, 10-20 cm) después de 14 afos de fertilización continuada. Los objetivos fueron evaluar 
el efecto de los distintos regímenes de fertilización de largo plazo en los distintos suelos y 
rotaciones de la región pampeana sobre: a. los contenidos de MOS y Nan y, b. la relación 
entre el Nan y la MOS. 


MATERIALES Y MÉTODOS 


Se tomaron muestras de suelo de tres estratos de profundidad (0-5 cm, 5-10 em 
y 10-20 cm) de los cinco sitios de la Red de Nutrición del CREA de la Región Sur de 
Santa Fe (Tabla 1), durante la campafia 2013-2014. Las características de los ensayos 
son: dos sitios con rotación maiz-trigo/soja (M-T/S) (Balducci y San Alfredo) y tres 
sitios con rotación maiz-soja-trigo/soja (M-S-T/S) (La Blanca, La Hansa y Lambaré). 
Se evaluaron los seis tratamientos de la red: testigo sin fertilizar, NP, NS, PS, NPS y 
completo (NPS+micronutrientes). El disefo es en bloques completos aleatorizados con 
tres repeticiones en cuatro sitios, y dos en San Alfredo. El tamafio de las parcelas es de 
25-30 m de ancho y 65-70 m de largo. Se utiliza maquinaria del productor y labranza en 
siembra directa (SD) en todos los casos. La dosis de fertilización se estima a partir de los 
rendimientos esperados más un 5-10% de fertilización de enriquecimiento con Py Sa la 
siembra de soja y, a partir de modelos zonales para cultivos de alto rendimiento en el caso 
del N. En la campafia 2013/2014 se aplicaron 30 y 18 kg de Py S respectivamente al cultivo 
de soja de primera en la rotación M-S-T/S y, en la rotación M-T/S se aplicaron 102, 44 y 21 
kg de N, Py S al cultivo de trigo (Boxler et al., 2014). 

Se determinó el N potencialmente mineralizable (Nan) por incubación anaeróbica 
a 40ºC por 7 días, método descripto por Keeney & Bremner (1966) seguido de una 
colorimetría. La MO se determinó por la metodologia propuesta por Walkley & Black (1934). 
Los resultados se analizaron mediante ANOVA, comparación múltiple de LSD y regresión. 
Los ANOVA de concentración de MOS y Nan se analizaron para cada estrato por separado. 
Las diferencias significativas fueron determinadas a un nivel de significancia del 95%. La 
relación entre Nan y MOS se analizó mediante regresión lineal. Se compararon las distintas 
regresiones mediante test de F y, en los casos en que no fueron diferentes, las líneas de 
regresión para esos tratamientos se representaron con una sola función. 
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Sitio Balducchi Alfredo La Blanca La Hansa Lambaré 
Suelo Hapludol Arguidol Hapludol Argiudol Argialbol 
típico tipico típico típico típico 
: Santa a : ; 
Serie Isabel Hughes La Bélgica Bustinza El Trébol 
Labranza siembra directa 
Aros agricultura previo 
a la instalación del + 60 15 12 +20 6 
ensayo 
Rotación M-T/S M-S-T/S 
Antecesor Trigo Maiz 





Textura de los suelos (0-20 cm de profundidad) 





Análisis Unidad 











Arcilla Yo 11,8 18,0 15,5 18,0 20,5 
Limo Yo 58,2 62,0 56,4 78,9 76,5 
Arena Yo 35,1 20,0 28,1 3,1 3,0 





Tabla 1. Información de manejo y de sitio. Red de Nutrición CREA Sur de Santa Fe. Soja. 
Campafia 2013/14. 


RESULTADOS Y DISCUSIÓN 


Concentración de MOS y Nan en el suelo 


Los catorce afios de fertilización continua con N, P y S tuvieron un efecto sobre 
la MOS sólo en el estrato de 0-5 cm de profundidad, mientras que en el resto del perfil 
no hubo diferencia entre los tratamientos de fertilización (Figura 1). La disminución de la 
concentración de MOS por no fertilizar (testigo) fue del 12% en la capa de 0-5 cm respecto 
del tratamiento completo Figura 1). La mayoría de autores han informado que los cambios 
en MOS se producen por efectos de manejo y por los afios de agricultura (Fabrizzi et 
al., 2003; Reussi Calvo, 2013; Genovese et al., 2009). Sin embargo, nuestros resultados 
sugieren que la fertilización también puede ser un factor importante para la regulación 
de los contenidos de MO y por consiguiente para mantener la fertilidad de los suelos en 
el largo plazo. Se han visto disminuciones de los rendimientos con las pérdidas de MOS 
(Ladha et al., 2003), ya que son fuente importante de nutrientes y sustrato de la actividad 
microbiana. En este trabajo, las mayores concentraciones de MOS se vieron en el estrato 
de 0-5 cm en todos los sitios (Figuras 1 y 2) concordante con Wander (2004); Diovisalvi et 
al. (2008) y (Diaz Zorita, 1999). 

San Alfredo y Lambaré presentaron los valores más altos de MOS en todo el perfil 
(promedio 31,5 g kg! y 30,2 g kg"! respectivamente) y Balducchi fue el más pobre en MOS 
(promedio 23,5 g kg')(Figura 2). Estas diferencias entre sitios son previas a la instalación de 
los experimentos en el afio 2000 (Boxler et a!., 2014) y se podrían atribuir a las diferencias 
en afios de agricultura y a las diferencias texturales en los sitios (Tabla 1). La información 
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de los sitios previo a la implementación de los ensayos indican que Balducchi tenía más 
de 60 aros de agricultura continua mientras que Lambaré y San Alfredo tenían 12 y 8 afios 
respectivamente (Tabla 1). Por lo tanto, los contenido de MOS y los afios de agricultura 
en los sitios guardan relación inversa entre sí y, son visibles en el largo plazo (Genovese 
et al. 2009, Reussi Calvo et al., 2013 y Eiza et al., 2005; Divito et al., 2011). Por otro 
lado, Lambaré es el sitio más rico en arcillas y Balducchi el más arenoso. En cuanto a los 
contenidos de (arcilla+limo), Lambaré tiene un 97%, mientras que Balducchi posee un 68% 
de dicha fracción (Tabla 1). Está documentado que la fracción (arcilla+limo) guarda relación 
positiva con los contenidos de COT de los suelos (Galantini et al., 2004), porque en suelos 
más finos el tamafio pequefio de los poros del suelo protegen a la MO de la acción de los 
microorganismos (Van Veen & Kuikman, 1990). Hassink et al. (1993a) indicaron que en 
suelos arenosos, la presencia de grandes poros disminuye la capacidad de retener agua y 
de mantener la humedad para el crecimiento bacteriano. Los cuales explican las mayores 
concentraciones de MOS total en suelos de texturas finas y menos laboreadas en este 
estudio. 

La fertilización con N, P y S durante 14 afios no cambió el contenido de N 
potencialmente mineralizable (Nan) en ninguno de los tres estratos de profundidad (Figura 
3). Fabrizzi et al. (2003) también informaron que el Nan no fue afectado por la fertilización 
con N en dos experimentos en el sudeste de la provincia de Buenos Aires, pero variaba con 
la MOS por efecto de los sistemas de labranza. Diovisalvi et a/. (2008) también informaron 
que la fertilización no afectó los contenidos de Nan en un ensayo de largo plazo. Tang et al. 
(2006) vieron resultados similares en un ensayo de fertilización continua de 15 afios en el 
norte de China, pero aplicaciones orgánicas aumentaban la MOS, y consecuentemente el 
N total y el N como amonio. Manna et al. (2007) también sugieren que se debe incrementar 
la MO si se desea aumentar la disponibilidad de N. En este estudio, el tipo de labranza 
homogénea (SD) en todos los sitios, y la ausencia de aplicación de N al cultivo de soja 
sumado a la alta frecuencia de este cultivo en las rotaciones pudo haber causado que 
las diferencias en los aportes de C al suelo entre los tratamientos no hayan sido lo 
suficientemente grandes como para provocar una diferencia en el Nan. 

El Nan varió con la profundidady conlos sitios (Figuras 3 y 4). Los mayores contenidos 
de este nutriente se vieron en el estrato de O a 5 cm y disminuyeron con la profundidad 
(Figuras 3 y 4). Esto puede atribuirse a la incorporación de residuos con el sistema de SD 
cerca de la superficie, lo cual aumentó la MOS y consigo el Nan en este estrato. Según 
West & Post (2002) la SD generalmente incrementa la MOS en comparación con la LC y es 
más notorio en el estrato superficial. Los valores más bajos en los 3 estratos de profundidad 
correspondieron a Balducchi (promedio 33,85 mg kg!) y Lambaré (promedio 65,13 mg kg!) 
fue el más rico en concentración de Nan (Figura 4). Vale aclarar que los contenidos de COT 
siguen el mismo patrón (Figura 2, Tabla 1). Nuestros resultados sugieren que el Nan está 
relacionado con la MOS, ya que proviene de mineralización del N orgánico. Esto puede 
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atribuirse a los contenidos de MO sumado a las condiciones ambientales favorables para 
la actividad microbiana que favorecieron a la mineralización y disponibilidad de N (Manna 
et al., 2007). 

La secuencia de cultivos de ambas rotaciones tampoco afectó a las concentraciones 
de Nan (Figura 2), coincidiendo con lo reportado por Gregorutti et al. (2014) en varias 
rotaciones que incluyeron a los cultivos de soja, maíz y trigo. Las Figuras 3 y 4 incluyen 
valores de Nan entre 7,32 y 128,44 mg kg”. Estos valores son consistentes con aquellos 
reportados por Reussi Calvo et al. (2014) quienes mostraron valores de Nan entre 12 y 
260 mg kg siendo mayores en el este que en el oeste de la provincia de Buenos Aires. 
Benitende et al. (2007) también reportaron valores de Nan promedio de 183 mg kg! y 104 
mg kg! para suelos Argiudoles ácuicos y vérticos de Entre Ríos. Fabrizzi et al. (2003) 
informaron concentraciones de Nan bajo el sistema SD en promedio de 61,5 mg kg! y 
bajo labranza convencional (LC) de 24,2 mg kg, lo cual indica que tanto la MOS como 
el Nan son afectados por el manejo del suelo. Se han visto que las prácticas de labranza 
afectan principalmente al pool de N en el suelo (Sharifi et al., 2008), pero el agregado de 
N como fertilizante no tiene efectos sobre el Nan (Diovisalvi et a!., 2008). Sin embargo, en 
situaciones de altos niveles de fertilización con N pueden ocurrir disminuciones de Nan 
debido al incremento del reciclaje de la MO sugiriendo pérdidas de N por inmovilización 
(Genovese et al., 2009) lo cual no ocurrió en este estudio. Por lo tanto, la estimación de Nan 
a partir de la MOS podría ser un indicador de la capacidad del suelo de mineralizar N que 
varía ampliamente en los suelos. En las evaluaciones del Nan se deben considerar otras 
características como las historias agrícolas de cada lote, sistema de labranza, fertilización, 
textura y cantidad así como la calidad de residuos de cosecha que se incorporan. 
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Figura 1: Distribución de la concentración de 
MOS a 0-20 cm en los tratamientos de la Red 
de Nutrición CREA Sur de Santa Fe, campafia 
2013/2014. Los círculos encierran tratamientos 

que no son diferentes significativamente 
(p>0.05) y las barras son el error estándar. Los 
cambios de letra indican diferencia significativa 
entre estratos de profundidad. 
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Figura 3: Distribución de la concentración de 
Nan a 0-20 cm en los tratamientos de la Red 
de Nutrición CREA Sur de Santa Fe, campafia 
2013/2014. Los círculos encierran tratamientos 
que no son diferentes significativamente 
(p>0.05) y las barras son el error estándar. Los 
cambios de letra indican diferencia significativa 
entre estratos de profundidad. 
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Figura 2: Distribución de la concentración de 
MOS a 0-20 cm en los sitios de la Red de 
Nutrición CREA Sur de Santa Fe, campafia 
2013/2014. Los círculos encierran sitios que 
no son diferentes significativamente (p>0,05) y 
las barras son el error estándar. Los cambios 
de letra indican diferencia significativa entre 
estratos de profundidad. 
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Figura 4: Distribución de la concentración 
de Nan a 0-20 cm en los sitios de la Red de 
Nutrición CREA Sur de Santa Fe, campafia 
2013/2014. Los círculos encierran sitios que 

no son diferentes significativamente (p>0,05) y 
las barras son el error estándar. Los cambios 
de letra indican diferencia significativa entre 
estratos de profundidad. 


Capítulo 10 


Relación entre Nan y MOS 


La Figura 5 muestra la relación entre el Nan y la MOS (p<0.01). Fabrizzi et al. 
(20083) afirmaron que el C de la MOP está altamente relacionado con el Nan, sus resultados 
sugieren que el Nan proviene de la mineralización de la MOS. Otros autores también han 
reportado la asociación entre MOS y el Nan (Sharifi et al., 2008; Genovese, 2009; Reussi 
Calvo et al., 2014). Contrariamente, Sainz Rosas et al. (2008) estudiaron 26 suelos de tipo 
Argiudol típico y Paleudol petrocálcico (3 sitios) en Balcarce e informaron que el COT no se 
relacionaba con el Nan, por lo cual indicaron que el COT no sería un indicador sensible de 
la capacidad de mineralización de N del suelo. Nuestros resultados también muestran que 
existe diferencias en la relación Nan: MOS entre los sitios. Pese a que los sitios más ricos en 
MOS fueron San Alfredo y Lambaré (Figura 2) no se vieron los mayores contenidos de Nan 
en San Alfredo, sino sólo en Lambaré (Figura 4). Lo cual se refleja en la relación Nan: MOS 
enla Figura 5, donde Lambaré tuvo los mayores valores de la relación Nan:MOS, Balducchi 
los valores menores y los otros 3 sitios presentaron relaciones Nan:MOS intermedias y 
estadísticamente no distintas entre sí. Estos resultados sugieren que el contenido de MOS 
no es un indicador suficiente para predecir la capacidad de proveer N disponible en los 
suelos. Existirían otros factores como la textura, pH, y el contenido de MO lábil que juegan 
papeles importantes en la mineralización de la MOS más allá de los contenidos totales de 
MOS (Baath & Anderson, 2003). Petersen et al. (2013) mostraron una relación positiva y 
muy estrecha entre el Nan y la biomasa microbiana en el suelo. Por otro lado, Diovsalvi et al. 
(2010a) reportaron que existe menor asociación entre el Nan y la MOS en suelos de textura 
fina. Sin embargo, nuestros resultados indican lo contrario, Lambaré pese a ser un sitio de 
textura fina mostró mayor correlación entre el Nan y MOS (Figura 5, r?: 0,46). Esto podría 
ser explicado por las condiciones favorables para la actividad microbiana en Lambaré más 
allá de su riqueza en MO. Además, Kanazawa & Filip (1986) también informaron que la 
mayor parte de los microorganismos se acumulan en la fracción arcillo+limo. 

Los regímenes de fertilización con N, Py S no afectaron a la relación entre el Nan 
y la MOS (Figura 5). Esto podría ser explicado porque no hubo grandes efectos de la 
fertilización sobre los contenidos de MOS. Esto sucedió sólo en el estrato de 0-5 cm, donde 
el testigo fue inferior en MO (Figura 1), dicho efecto no fue significativo sobre la relación 
Nan:MOS. Está documentado que los cambios más notables en contenido de MOS y Nan 
se pueden ver por efectos de manejo y afios de agricultura (Genovese et al., 2009, Reussi 
Calvo etal., 2013 y Eiza et al., 2005; Divito et al., 2011), por el contrario no hay reportes de 
los efectos puros de la fertilización sobre la MO y por consiguiente sobre el Nan. En este 
estudio, el sistema de labranza fue SD en todos los sitios y la frecuencia del cultivo de soja 
en las rotaciones pudieron haber impedido las diferencias significativas del contenido del 
Nan con los tratamientos. Por lo cual, coincidimos con lo que afirmaron Marriott & Wander 
(2006) y Sainz Rosas et al. (2008) que el Nan no es una fracción que varía en paralelo con 
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los contenidos de MOS, sino que existe diferencias dependiendo de las potencialidades de 


mineralización de los suelos y condiciones favorables para la biomasa microbiana. Por lo 


cual no se pueden predecir la capacidad de un suelo de mineralizar N sólo a partir de los 


análisis de MOS total. 
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Figura 5: Relación entre concentración de Nan y la MOS. La línea verde corresponde a la 
función de Nan:MOS en Lambaré, la línea negra a San Alfredo, La Blanca y La Hansa, y la 





La 
Hansa Lambaré 


línea azul corresponde a Balducchi. 


CONCLUSIONES 


Las fuentes de variación evaluadas en estos experimentos (i. e. fertilización con N, 


Py S) mostraron diferencias en los contenidos de MOS sólo en el estrato de 0-5 cm de 


profundidad. Sin embargo, no afectaron a los contenidos de Nan en el suelo. Las rotaciones 


no afectaron a los contenidos de MOS y Nan. 
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Se determinaron relaciones positivas entre el Nan y la MOS y sus variaciones fueron 
atripuibles a las diferencias en otras características además del contenido de MOS total 
entre los sitios. Por lo tanto, solo la MOS total no es un buen predictor de la capacidad de 


mineralización de N de los suelos. 
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RESUMO: O objetivo do presente trabalho foi 
analisar as alterações em atributos químicos 
de um Latossolo Amarelo do bioma Amazônia, 
submetido ao sistema agroflorestal (SAF's) na 
agricultura familiar, sob testemunha de uma área 
de capoeira terciaria com presença de floresta 
remanescente, na região do sudeste paraense, 
município de Marabá. As coletas das amostras 
de solo foram realizadas de forma manual, no 
processo foram abertas pequenas trincheiras no 
formato circular com o auxílio de uma cavadeira 
manual até a profundidade pretendida, e esta 
foi mensurada com a utilização de uma trena 
métrica, posteriormente as amostras de cada 
local eram homogeneizadas em um só recipiente. 
Os resultados foram submetidos a análises 
estatísticas e conclui-se que as concentrações 
dos atributos químicos colaboraram para a 
avaliação da qualidade do solo, subsidiando 
a comprovação de manejos de solo mais 
sustentáveis, desse modo, colaborando para a 
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recuperação de áreas degradadas. 
PALAVRAS-CHAVE: Agricultura familiar, solo, nutrientes. 


HEMICAL ATTRIBUTES OF THE SOIL IN THE SURFACE LAYER UNDER 
DIFFERENT CULTIVATION SYSTEMS IN SOUTHEAST PARAENSE 


ABSTRACT: The objective of the present work was to analyze the changes in chemical 
attributes of a Yellow Latosol of the Amazon biome, submitted to the agroforestry system 
(SAF's) in family farming, under witness of a tertiary capoeira area with the presence of 
remaining forest, in the southeastern region of Pará , municipality of Marabá. The soil samples 
were collected manually, in the process small circular trenches were opened with the aid 
of a manual digger to the desired depth, and this was measured with the use of a metric 
measuring tape, afterwards the samples of each were homogenized in one container. The 
results were subjected to statistical analysis and it was concluded that the concentrations of 
chemical attributes collaborated for the evaluation of soil quality, supporting the evidence of 
more sustainable soil managements, thus collaborating for the recovery of degraded areas. 
KEYWORDS: family farming, soil, nutrients. 


11 INTRODUÇÃO 


O aumento da produção de alimentos nos agrossistemas acarreta a alteração do 
solo resultando na perda de nutrientes e gerando uma deterioração de sua qualidade. 
Assim, novos conceitos e procedimentos agrícolas surgem em busca de sustentabilidade, 
dentre eles o sistema agroflorestal que pode promover um desenvolvimento rural mais 
sustentável, principalmente para a agricultura familiar (PARTELLI ET AL., 2006). 

Com isso, os Sistemas Agroflorestais — SAF's são vistos como opção para 
minimização do desmatamento por quebrarem a predominância do ciclo de agricultura 
migratória e pecuária extensiva praticadas na Amazônia, sendo opção para gerar lucros 
significativos em áreas relativamente pequenas (SERRÃO, 1995; SMITH ET AL., 1998). 
Outro aspecto é o da diversificação na produção, fator contribuinte para obtenção de renda 
contínua, o que pode ser um critério para a adoção desse sistema. 

No estado do Pará, em especial no sul e sudeste, existem poucos incentivos à 
implantação desses sistemas de multe-espécies madeireiras e não-madeireiras. Isso se dá 
por conta das características culturais do estado que se limitam à plantação de espécies 
agrícolas. 

Nesse contexto, a importância do trabalho é a avaliação dos atributos físicos, 
químicos e biológicos do solo em ambientes degradados pela agricultura exorbitante, é de 
extrema importância devido à sua sensibilidade às alterações na qualidade do solo, uma 
vez que pode fornecer subsídios para o estabelecimento de sistemas racionais de manejo 
e contribuir para a manutenção de ecossistemas sustentáveis (CARNEIRO ET AL., 2009). 
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21 MATERIAL E MÉTODOS 


2.1 Descrição da área de estudo 


O trabalho foi realizado no PDS Porto Seguro, no município de Marabá, região 
sudeste do Pará. O local está sob as coordenadas geográficas 5º28'39.55"S e 49º 2'41.86"0, 
138 m de altitude, o solo do local foi classificado como Latossolo. O clima do município de 
Marabá é caracterizado como equatorial, quente e úmido, segundo a classificação Kôppen, 
apresentando temperaturas médias mensais entre 22,9ºC e 32ºC, com média anual de 
26ºC, umidade relativa do ar variável na faixa entre 73% a 93%, e precipitação anual de 
cerca de 1.976 mm. A época de maior precipitação compreende o período entre janeiro e 
março, e o mais seco indo de julho a setembro (RAIOL, 2010) 
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Figura 1. Localização do PDS Porto Seguro e pontos de amostragem. 
Fonte: Autores, 2021. 


2.2 Coleta das amostras de solo 


Para melhor representatividade da área de estudo, foram consideradas duas 
áreas para aplicação do mesmo, com os seguintes sistemas de uso do solo: Sistema 
Agroflorestal (SAF's), com diferentes anos de adoção do sistema, e uma parcela da área de 
capoeira secundária em fase de regeneração natural com histórico de pastagem (Capoeira 
secundária - CS) para servir de padrão de comparação dos teores naturais dos atributos 
químicos do solo do local em estudo (Figura 1). No processo foram abertas pequenas 
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trincheiras no formato circular com o auxílio de uma cavadeira manual até a profundidade 
pretendida, e esta foi mensurada com a utilização de uma trena métrica. Posteriormente as 
amostras de cada local eram homogeneizadas em um só recipiente. 

Nestas amostras buscou-se avaliar os efeitos destes diferentes agroecossistemas 
quanto à qualidade química do solo. Para tal, foram analisados os macronutrientes primários 
e secundários, como pode ser observado na Figura 2 a seguir. 


2.3 Análise química do solo 


Foram determinados os atributos químicos para análise da fertilidade do solo (Figura 
2), conforme metodologia empregada pela Embrapa (2011). 








Atributo Método 
Teor de Matéria Orgânica Determinado pelo método de Walkley e Black (1934) 
Nitrogênio Total (N) Método de Kjeldahl por destilação a vapor 


O pH do solo foi medido por meio de eletrodo combinado 
Determinação do pH (potenciômetro) imerso em suspensão solo: líquido (água, KCI ou 
CaCl2). 





O cálcio e o potássio foram extraídos com a solução de Mehlich 1 
(HCla 0,5 N + H, SO, a 0,025 N) e determinados em colorímetro 
e fotômetro de chama, respectivamente. 


Determinação de Sódio (Na*) e 
Potássio (K*) 





O fósforo foi extraído com a solução de Mehlich e determinado 


Determinação de Fósforo (P) por espectrofotometria (Biospectro — SP - 220). 


Cálcio e Magnésio trocáveis O Ca? e o Mg? foram extraídos em KCl a 1N e determinados por 
(Ca?* e o Mg?) titulação de EDTA 0,025 N. 





Foram determinados por titulometria, usando-se solução de 
acetato de cálcio 1N a pH 7 para sua extração com EDTA- Na 
(acidez potencial). 


Acidez potencial: Hidrogênio 
(H*) + Alumínio (AI **) 





Capacidade de Troca de 
Cátions ou CTC (T) em pH 7 e 
a Saturação por Bases (V). 


T= (Ca? + Mg?* + Nat + K*) + (H* + AS!) 
V= 100 (Ca? + Mg? + Nat + K)/T 





Figura 2 - Atributos químicos e métodos de análise para cada amostra. 
Fonte: Autores, 2021. 


2.4 Delineamento estatístico 


A análise de dados foi feita por meio de estatística descritiva, submetidos a análises 
de variância (ANOVA). A comparação das médias foi realizada por meio do teste de Scott- 
Knott a 5% de probabilidade realizada no software estatístico SISVARO (FERREIRA, 2014) 
para verificar a diferença dos teores e taxas dos atributos estudados entre áreas escolhidas. 
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A partir das análises realizadas, serão discutidos os seguintes macronutrientes 
primários: pH, Matéria Orgânica (M.O), Nitrogênio (N), Fósforo (P) e Potássio (K), que são 
os fatores mais importantes no solo para o desenvolvimento vegetal. 


31 RESULTADOS E DISCUSSÕES 





Tratamento MO N K P Al Ca Ca+Mg pH H+Al Base 





I tiic 14,33b 43,53a 1,01a 0,94a 1,35a 1,7/3a 4,52a 3,59b 39,33a 
I 1,02c 12,89b 45,72a 1,94a 1,06a 1,00c 1,57a 4,73a 4,08b 24,16b 
I 1,09c 42,73a 35,94a 3,93a 0,95a 0,98c 1,21b 483a 5,395a 24,00b 
IH 1,66a 8,63b 56,65a 4,90a 0,89a 1,14b 1,55a 5,17a 397b 34,00a 
IH 1,51a 11,86b 39,31a 4,35a 0,70b 0,91c 1,50a 4,28a 3,98b 29,48b 


IH 1,33b 36,08a 47,86a 0,91a 1,06a 0,89c 1,21b 441a 4,66a 26,17b 





CM. (%) 8,56 42,38 11,68 17,96 9,82 6,36 12,76 7,16 12,61 13,7 





| - Capoeira secundária - CS e Il - SAF's. 


Tabela 01. Atributos químicos da fertilidade do solo em um Latossolo em Marabá, PA, na 
profundidade de 0-0,20 cm 


Fonte: Autores, 2021. 


As amostras do solo apresentaram valores baixos de pH em água, evidenciando 
a acidez elevada em todos os sistemas (Tabela 1). A correção do subsolo ácido pode ser 
feita por meio da calagem profunda. Contudo, essa prática necessita de revolvimento do 
solo, razão por que não é de interesse em áreas já estabelecidas com sistema de cultivo 
que não envolve preparo convencional. Além disso, a calagem profunda exige máquinas 
potentes e equipamentos caros, o que torna a prática onerosa. Em condições de extrema 
acidez (valores de pH muito baixos e altos teores de Al trocável) pode ocorrer limitação na 
decomposição da MO, resultando no seu acúmulo. 

Em relação a quantidade de N no solo apenas dois resultados se diferiram no 
tratamento, e esses resultados revelaram um alto teor. Gonçalves et al. (2001) e D'Andréa 
et al. (2004) afirmam que as maiores reservas de N estão na camada de 0-20 cm do 
solo, concordando com os resultados obtidos. Os outros resultados mais baixos foram 
no sistema SAF'S, isso acontece, pois a perda de N está provavelmente relacionada ao 
maior revolvimento do solo provocado pelos anos de cultivos daquela área. Vários estudos 
apontam correlação negativa entre revolvimento do solo e perda de N (SOLOMON ET AL., 
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2000; SMITH ET AL., 2001; ZIBILSKE ET AL., 2002; MIELNICZUK ET AL., 2003). 

Os teores de fósforo no tratamento 1 e no tratamento 2 não se diferem (Tabela 1), 
foram semelhante à condição verificada na maioria dos solos brasileiros (RAIJ, 1991). O 
baixo teor de fósforo do solo é em função da formação de compostos insolúveis com o 
cálcio. Os solos brasileiros são carentes de P, em consequência do material de origem e 
da forte interação do P com o solo (Raij, 1991), em que menos de 0,1% encontra-se em 
solução (FARDEAU, 1996). 

Os teores de K disponível no solo foram originalmente altos (tabela 1), conforme 
CQFS-RS/SC (2004). Em relação, o fenômeno que está acontecendo nos dois tratamentos, 
são positivos. Assim, os níveis de K relativamente possuem médias aproximadas nos 
dois sistemas. Na solução do solo, o potássio é móvel e, também, sujeito às perdas por 
lixiviação, ocorrendo principalmente em solos ácidos e com baixa CTC, erosão, remoção 
por colheitas (GIRACCA ET/AL., 2016). 

A MO possui a função de complexar os cátions H+ e A+3 livres com compostos 
orgânicos aniônicos dos resíduos e adicionar bases (Ca, Mg e K) que reduzem a acidez do 
solo e aumentam o pH (Franchini et al., 2001; Pavinato, 2008). Os maiores níveis de MO 
foram observados em SAFs isso pode ser explicado porquê a ausência de revolvimento 
do solo nos SAF's e na mata nativa proporciona melhores condições para os organismos 
responsáveis pela fragmentação do material vegetal e ciclagem dos nutrientes (Cunha et 
al., 2012). Vários trabalhos têm evidenciado o potencial dos SAF's, sobretudo daqueles que 
utilizam o cultivo em aleias na ciclagem de nutrientes (FRANZEL ET AL., 2001), que atua 
no aumento dos teores de MO e seus componentes, como C, Pe N. Assim, o SAF induz 
o aumento dos níveis de C do solo em virtude do maior aporte de resíduos que o sistema 
recebe ao longo dos anos (BARRETO & LIMA, 2006) 


41 CONCLUSÕES 


Ambientes de cultivo localizados no Sudeste Paraense apresentam-se em estádio 
avançado de degradação, conforme demonstrado por vários atributos físicos, químicos e 
biológicos destacando-se o uso excessivo do solo, acarretando o baixo teor de matéria 
orgânica e consequentemente a baixa produtividade da área. O sistema agroflorestal tem 
uma média estatística maior que a capoeira secundária, isso ocorre por conta que existe 
mais matéria orgânica no mesmo, fazendo com que a área da capoeira secundária seja 
limitada e apresente dificuldades para produção. Conclui-se que a implantação do sistema 
agroflorestal possibilitou a sucessão ecológica no processo de cultivo concomitante de 
diversas espécies de hortaliças, árvores e frutas. Ao mesmo tempo, contribui com a 
recuperação de áreas degradadas e com a conservação do solo e da água no sistema 
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RESUMO: O objetivo desse estudo foi determinar 
o diâmetro médio ponderado, os atributos 
químicos do solo, o carbono orgânico total e das 
frações granulométricas e químicas da matéria 
orgânica do solo (MOS) de agregados fisiogênicos 
e biogênicos formados em fragmentos florestais 
em distintas condições de relevo (pedoforma 
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convexa e côncava) no município de Pinheiral- 
RJ. Foram selecionadas duas pedoformas, 
a saber: convexo (convexa-divergente) (CX) 
e côncavo (côncava-convergente) (CC), em 
Floresta Estacional Semidecidual Submontana. 
Foram coletadas seis pseudorepetições (torrões) 
na camada de 0,00-0,10 m. As amostras foram 
secas ao ar, submetidas a tamisação, utilizando- 
se um conjunto de peneiras de 9,7 e 8,0 mm de 
malha, selecionando-se somente os agregados 
retidos neste intervalo, sendo separados 
em fisiogênicos e biogênicos. Nestes foram 
determinados o diâmetro médio ponderado 
(DMP) dos agregados, os atributos químicos 
(pH, Ca”, Mg?*, AIS*, H+AI, K e P), os valores do 
complexo sortivo (S, T e V), o carbono orgânico 
total (COT) eo das frações granulométricas (COP, 
particulado; e COAM, associado aos minerais) 
e químicas (C-FAF, ácido fúlvico; C-FAH, ácido 
húmico; C-FHUM, humina; e %SH: porcentagem 
de carbono das substâncias húmicas) da MOS. 
De maneira geral, a maior variabilidade entre os 
atributos foi verificada em CX, principalmente 
nos agregados biogênicos, com destaque para 
os relacionados a MOS. O estudo quantitativo da 
MOS a partir da compartimentalização do carbono 
orgânico nos diferentes tipos de agregados, 
pode ser uma técnica eficaz na mensuração da 
qualidade do solo mesmo em ambientes mais 
estáveis e equilibrados. 

PALAVRAS-CHAVE: Vias de formação de 
agregados, compartimentalização do carbono 
orgânico, qualidade do solo. 


Capítulo 12 


AGGREGATION UNDER DIFFERENT LANDFORMS IN SEMIDECIDUOUS 
FOREST IN SOUTHEASTERN BRAZIL 

ABSTRACT: The objective of this study was to determine the weighted average diameter, the 
soil chemical attributes, the total organic carbon and the organic carbon of the granulometric 
and chemical fractions of soil organic matter (SOM) in physicogenic and biogenic aggregates 
formed in forest fragments in different relief conditions (convex and concave pedoform) in 
the municipality of Pinheiral-RJ. Two pedoforms were selected, namely: convex (convex- 
divergent) (CX) and concave (concave-convergent) (CC), in Submontane Semideciduous 
Seasonal Forest. Six pseudo-repetitions (clods) were collected in the 0.00-0.10 m layer. After 
collection, the samples were dried in the air and subsequently subjected to sieving, using a 
set of sieves of 9.7 and 8.0 mm mesh, selecting only the aggregates retained in this interval. 
These were identified and separated into biogenic and physicogenic. Subsequently, the 
weighted average diameter (DMP), the chemical attributes (pH, Ca?+, Mg?*, Alº+, H+AI, K* and 
P), the values of the sortive complex (S, T and V), the total organic carbon were determined 
(TOC) and organic carbon of particle size (COP, particulate; and COAM, associated with 
minerals) and chemical (C-FAF, fulvic acid; C-FAH, humic acid; C-FHUM, humine; and% 
SH: percentage carbon from humic substances) from soil organic matter (SOM). In general, 
most of the variability between the attributes evaluated was verified in CX, mainly in biogenic 
aggregates, with emphasis on those related to SOM. The quantitative study of SOM from the 
compartmentalization of organic carbon in the different morphological types of aggregates, 
can be an effective technique in measuring soil quality even in more stable and balanced 
environments. 

KEYWORDS: Aggregation formation pathways, compartmentalization of organic carbon, soil 
quality. 


11 INTRODUÇÃO 


Na região do Médio Vale do Paraíba do Sul, Pinheiral — RJ, a Floresta Atlântica é 
representada pela Floresta Estacional Semidecidual Submontana que historicamente foi 
de grande importância para o desenvolvimento econômico do país (Santos et al., 2016). 
Atualmente encontra-se com extensas áreas ocupadas por pastagens degradadas e de 
baixa produtividade. 

A região é caracterizada por apresentar um relevo declivoso com diferentes feições 
topográficas denominadas de pedoformas (côncavas, retilineas e convexas). Associadas a 
outros fatores bióticos e abióticos, essas feições, produzem uma variedade de condições 
ambientais, tais como: solos, transporte de sedimentos, precipitação e gradiente de 
umidade, luminosidade, organização vertical do dossel e geração de aspectos de degraus 
no estrato arbóreo (Santos et al., 2017), levando à formação de diversos microambientes 
na paisagem. 

Nesses microambientes, considerados mais equilibrados e estáveis, em virtude das 
características topográficas específicas e das classes de solos, a dinâmica (deposição, 
acúmulo e decomposição) do material orgânico formador de serapilheira acontece de forma 
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diferenciada, influenciando nos atributos físicos (estrutura), químicas (fertilidade natural) e 
biológicas (biomassa microbiana) do solo. 

Para avaliar tais ambientes, alguns atributos edáficos têm sido utilizados como 
indicadores de qualidade do solo. Nesse sentido, destacam-se a agregação (estabilidade, 
gênese ou vias de formação) e a matéria orgânica do solo (MOS) (acúmulo, proteção ou 
compartimentos). Os mecanismos de formação dos agregados envolvem processos físicos 
e químicos (agregados fisiogênicos), e biológicos (agregados biogênicos) (Loss et al., 
2014). 

A diferenciação entre esses tipos de agregados é realizada de acordo com sua 
gênese ou vias de formação a partir de padrões morfológicos já pré-estabelecidos (Bullock 
et al., 1985; Pulleman et al., 2005; Batista et al., 2013). A MOS é considera o indicador 
chave na avaliação da qualidade do solo, e apresenta uma correlação positiva com os 
mecanismos de formação e estabilização dos agregados do solo. 

A partir do exposto, este estudo teve como hipótese testada que diferentes 
pedoformas influenciam nas vias de formação dos agregados, assim como no seu grau 
de estabilidade, na sua composição química e no conteúdo de carbono orgânico em seus 
diferentes compartimentos. Desta forma, objetivou-se avaliar o diâmetro médio ponderado, 
os atributos químicos, o carbono orgânico total e das frações granulométricas e húmicas 
da MOS em agregados fisiogênicos e biogênicos formados em fragmentos florestais em 
distintas condições de relevo (pedoforma convexa e côncava) no município de Pinheiral-RJ. 


21 MATERIAL & MÉTODOS 


O estudo foi realizado no município de Pinheiral, na região do Médio Vale do Paraíba 
do Sul, Rio de Janeiro, Sudeste do Brasil, na sub-bacia do Ribeirão Cachimbal, que compõe 
a bacia hidrográfica do rio Paraíba do Sul, localizado entre as coordenadas 22º29'03” S a 
22º35'27" S e 43º54'49" W a 44º04'05” W. 

A área de estudo está inserida em um fragmento florestal com aproximadamente 
204 ha, em área de relevo que varia de ondulado a forte ondulado e altitudes variando 
de 360 a 720 m. Segundo Kôppen o clima da região pode ser classificado como Cwa — 
clima temperado de inverno seco e verão chuvoso e Am — clima tropical chuvoso com 
inverno seco. As médias anuais de precipitação e temperatura são de 1,117 mme 20,9º C, 
respectivamente e o solo foi classificado como Cambissolo Háplico Tb distrófico (Santos 
et al., 2018). 

Foram selecionadas duas pedoformas adjacentes com formato do tipo convexo 
(convexa-divergente) (CX) e côncavo (côncava-convergente) (CC), com orientação da 
vertente predominante a sudeste e cobertura vegetal tipo Floresta Estacional Semidecidual 
Submontana. Cada pedoforma foi subdividida em três terços (superior, médio e 
inferior), sendo em cada um destes coletadas duas amostras compostas indeformadas 
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(pseudorepetições), formadas a partir de três amostras indeformadas simples (torrões) na 
camada de 0,00-0,10 m de profundidade mediante abertura de trincheiras, totalizando seis 
pseudorepetições por pedoforma. 

Após a coleta, as amostras foram secas ao ar e, submetidas a tamisação, sendo 
utilizado um conjunto de peneiras de 9,7 e 8,0 mm de malha, selecionando-se somente os 
agregados retidos neste intervalo. Os agregados foram examinados sob lupa e separados 
em frações morfológicas, identificando-se duas classes (fisiogênicos e biogênicos) 
mediante método adaptado de Pulleman et al. (2005) a partir dos padrões morfológicos 
estabelecidos por Bullock et al. (1985) e validados por outros autores (Batista et al., 2013; 
Loss et al., 2014, 2017). 

A diferenciação foi feita através de padrões morfológicos, obedecendo aos 
seguintes critérios: agregados biogênicos - aqueles em que é possível a visualização 
de formas arredondadas, providas do trato intestinal dos indivíduos da macrofauna do 
solo, principalmente Oligochaeta (minhocas) e/ou aqueles em que é possível visualizar a 
presença e atividade de raízes; e agregados fisiogênicos - aqueles que possuem formas 
angulares resultantes da interação entre carbono, argila, cátions e ciclos de umedecimento 
e secagem do solo. 

Após a identificação e separação, os diferentes tipos de agregados foram submetidos 
à análise de estabilidade via úmida, utilizando-se tamisação vertical, com conjunto de 
peneiras com malha decrescente, a saber: 2,00, 1,00, 0,50, 0,250 e 0,105 mm, por 15 min, 
no aparelho de Yooder. O material retido em cada peneira foi transferido para placas de 
Petri, e seco em estufa, a 105 ºC, até obtenção de massa constante. A partir dos dados de 
massa de agregados, foi calculado seu diâmetro médio ponderado (DMP) dos agregados 
de acordo com Teixeira et al. (2017). 

Para a realização das demais análises, uma outra amostra de agregados foi 
destorroada e passada por peneira de 2,00 mm de malha, obtendo-se assim a fração 
terra fina seca ao ar (TFSA). A partir da TFSA dos diferentes tipos de agregados foram 
realizadas as análises dos atributos químicos, carbono orgânico total (COT) e das frações 
granulométricas e húmicas da MOS. A quantificação dos teores de Ca?t, Mg?*, K+, AI, 
H+AI, P, pH em água, sendo calculados a soma de bases (S), a saturação por bases (V%) e 
a capacidade de troca catiônica (T) segundo Teixeira et al. (2017) e o COT foi determinado 
conforme Yeomans & Bremner (1988). 

Para o fracionamento físico granulométrico foi utilizado o método proposto por 
Cambardella & Elliot (1992), separando a MOS em duas frações, fração particulada 
(relacionada à fração areia do solo) e fração associada aos minerais (relativa às frações 
argila e silte). A determinação do carbono orgânico da fração particulada (COP) foi 
realizada segundo Yeomans & Bremner (1988) e o carbono orgânico da fração associada 
aos minerais (COAM) da MOS foi quantificado pela diferença entre o COT e o COP. 

As substâncias húmicas (SHs) foram obtidas pela extração e separação por meio 
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da solubilidade diferencial da matéria orgânica em meio básico ou ácido (ácidos fúlvicos e 
húmicos) e o resíduo (humina) (Benites et al., 2003). A determinação do carbono orgânico 
em cada uma das frações húmicas, ácidos fúlvicos (C-FAF), ácidos húmicos (C-FAH) e 
humina (C-FHUM), também foi realizada de acordo com Yeomans & Bremner (1988). 

Os resultados foram analisados quanto à normalidade e homocedasticidade dos 
dados pelos testes de Shapiro Wilk e Barttlet, respectivamente. Posteriormente, os resultados 
foram analisados como delineamento inteiramente casualizado, sendo submetidos à 
análise de variância com aplicação do teste F e os valores médios comparados entre si 
pelo teste de Tukey a 5% de probabilidade utilizando o programa R 3.3.1. 

Também foi realizada a análise de componentes principais (ACP), utilizada para 
reduzir as dimensões dos dados e, consequentemente, facilitar a análise por meio do 
gráfico do círculo de correlações, fornecendo com isso uma melhor interpretação do padrão 
dos agregados em função dos atributos avaliados por meio do programa PAST, e análise de 
correlação de Pearson pelo programa Excel. 


31 RESULTADOS & DISCUSSÃO 

Na Figura 1 são apresentados os resultados de diâmetro médio ponderados (DMP) 
dos agregados, sendo observadas diferenças somente entre as pedoformas. Em CX para 
os agregados fisiogênicos foram quantificados os maiores valores de DMP em comparação 
a CC. 
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Figura 1. Diâmetro médio ponderado de agregados fisiogênicos e biogênicos sob pedoforma 
convexa (CX) e côncavo (CC) em Floresta Estacional Semidecidual Submontana no Sudeste 
do Brasil, Pinheiral — RJ, na camada 0,00-0,10 m. 

Médias seguidas pela mesma letra maiúscula não diferem as pedoformas para o mesmo tipo de 
agregado, mesma letra minúscula não diferem os tipos de agregados para a mesma pedoforma 


avaliada (Teste de Tukey, p<0,05). CV1: Coeficiente de variação entre os tipos de agregados; e CV2: 
Coeficiente de variação entre as pedoformas. 
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Em estudo envolvendo agregados biogênicos e fisiogênicos de diferentes classes 
de solos do Estado da Paraíba sob cana-de-açúcar e pastagem de capim braquiária, Silva 
Neto et al. (2010) verificaram que os biogênicos apresentaram maior estabilidade física em 
comparação aos agregados fisiogênicos, como um reflexo da atividade biológica na gênese 
dos agregados biogênicos. Em áreas de floresta secundária em estágio avançado de 
sucessão (FSEA) e pasto misto manejado (PMM) em Pinheiral - RJ, Silva Neto et al. (2016) 
quantificaram maiores valores de DMP em agregados biogênicos quando comparados aos 
intermediários e fisiogênicos. 

Também em Pinheiral - RJ, Fernandes et al. (2017) verificaram maiores valores 
de DMP para agregados biogênicos em área de FSEA em relação aos demais tipos 
de agregados estudados. Mergen Junior et al. (2019 a) observaram que os agregados 
biogênicos em área de sistema plantio direto com aplicação de cama sobreposta de suínos 
na camada de 0,05-0,10 m apresentaram maiores valores de DMP em comparação aos 
fisiogênicos. 

Para Silva Neto et al. (2010), as maiores proporções de argila, cátions e matéria 
orgânica geralmente encontradas em agregados biogênicos (coprólitos) podem favorecer 
interações mais intensas, comparativamente aos agregados fisiogênicos, aumentando sua 
estabilidade. 

As gramíneas quando presentes, também contribuem nos processos de formação 
e estabilização dos agregados, devido, principalmente, à elevada densidade de suas 
raízes, que promove a aproximação de partículas pela constante absorção de água do 
solo, bem como, às periódicas renovações do sistema radicular e à uniforme distribuição 
dos exsudados orgânicos ao solo (Mergen Junior et al., 2019 a). Estes fatores, estimulam a 
atividade microbiana edáfica, cujos subprodutos produzidos também atuam na estabilização 
dos agregados (Brandão & Silva, 2012; Mergen Junior et al., 2019 a). 

Os maiores valores de DMP observados nos agregados fisiogênicos em CX em 
comparação à CC (Figura 1), podem estar relacionados às condição ambientais da 
pedoforma (maior incidência de radiação solar e ventos) em época chuvosa (maior oferta 
de umidade) associada a diversidade e estruturação do sistema radicular das espécies 
vegetais do ambiente (Santos et al., 2017). Visto que a absorção de água continuamente 
pelas raízes promove uma secagem local, contribuindo para aumentar a coesão entre as 
partículas do solo (Brandão & Silva, 2012), favorecendo a gênese de agregados fisiogênicos. 

Em estudo realizado em áreas sob diferentes sistemas de manejo agroecológico 
em Seropédica - RJ, Rossi et al. (2016) observaram que os resultados de DMP foram 
influenciados apenas nos sistemas manejo e não dos diferentes tipos de agregados 
avaliados. Também nesse sistema de manejo, Schuliz et al. (2019) não observaram 
diferenças nos valores de DMP entre as vias de formação avaliadas. 

Os resultados da caracterização dos atributos químicos nos diferentes tipos de 
agregados são apresentados na Tabela 1. As maiores variações foram verificadas com 
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maior expressão entre as pedoformas do que entre as vias de formação, principalmente 


para Ca?*, soma de bases trocáveis (S), total de cátions adsorvidos (T) e saturação por 


bases (V). 


Na comparação entre as pedoformas, em CX foram quantificados os maiores teores 


de Ca? e valores S, T e V em agregados fisiogênicos. Quanto aos tipos de agregados, 


verificou-se que nos biogênicos foram observados os maiores valores S e Vem CC, e nos 


fisiogênicos maiores valores T em CX (Tabela 1). 






































pH (H,0) Ca? (cmol, dm*) Mg?* (cmol, dm) 
Pedoformas 
Fis Bio Fis Bio Fis Bio 
CX 4,81 NS ns 4,85 NS ns 4,16 Ars 4,52 NS ns 3,20 NSns 3,78 NS ns 
CC 4,61 NS ns 4,60 NS ns 1,96 Brs 2,85 NS ns 3,26 NS ns 3,90 NS ns 
CVI% 6,5 32,4 
CV2% 14,6 35,7 
AIS* (cmol, dm?) H+AI (cmol, dm?) 
Pedoformas 
Fis Bio Fis Bio 
CX 0,86 NS ns 0,81 NS ns 8,34 NS ns 7,24 NS ns 
CC 1,05 NS ns 0,95 NS ns 9,12NS ns 8,21 NS ns 
CVI% 68,5 19,8 
CV2% 84,1 38,5 
K* (mg dm?) P (mg dm?) 
Pedoformas 
Fis Bio Fis Bio 
CX 103 NS ns 118NS ns 2,72 NS ns 2,25 NS ns 
CC 97 NS ns 106 NS ns 2,24 NS ns 2,62 NS ns 
CVI% 21,6 30,2 
CV2% 53,6 34,0 
S (cmol, dm?) T (cmol, dm?) V(%) 
Pedoformas 
Fis Bio Fis Bio Fis Bio 
CX 8,05 A"s 8,62 NSns 18,09Aa 15,86NSb 54 Ars 59 NSns 
cc 4,74Bb 7,12NS a 15,422 B"rs 14,24NSns 35Bb 50NSa 
CVI% 24,1 20,9 
CV2% 20,3 20,8 





Médias seguidas pela mesma letra maiúscula na coluna não diferem as pedoformas para o 
mesmo tipo de agregado, mesma letra minúscula na linha não diferem os tipos de agregados 
para a mesma pedoforma avaliado (teste de Tukey, a 5% de probabilidade) e "s Não 
significativo a 5% de probabilidade. S: Soma de bases trocáveis; T: Total de cátions adsorvidos, 
corresponde a CTC à pH 7,0; V: Saturação por bases; CV1: Coeficiente de variação entre os 

tipos de agregados; e CV2: Coeficiente de variação entre as pedoformas. 


Tabela 1. Caracterização dos atributos químicos de agregados fisiogênicos (Fis) e biogênicos 
(Bio) sob pedoforma convexa (CX) e côncavo (CC) em Floresta Estacional Semidecidual 
Submontana no Sudeste do Brasil, Pinheiral — RJ, na camada 0,00-0,10 m. 
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Em relação aos atributos químicos (Tabela 1), os maiores teores de Ca?* e valores 
S, Te V nos agregados fisiogênicos em CX estão relacionados as condições ambientais da 
pedoforma. Avaliando a influência da curvatura da superfície na variabilidade espacial dos 
atributos químicos do solo na região de Mar de Morros, Pinheiral - RJ, Santos et al. (2016) 
quantificaram maiores valores de pH, Ca?, Mg?*, P e K* na pedoforma côncava e teores 
mais elevados de AlS* na convexa. Esses resultados diferem dos verificados nesse estudo. 

Quanto aos tipos de agregados, para os resultados dos valores S e V entre as 
vias de formação (Figura 1) verificou-se padrão semelhante ao observado nos trabalhos 
de Silva Neto et al. (2010) e Melo et al. (2019). Os agregados biogênicos, são formados 
por processos biológicos, o que contribui para diferenças nos atributos químicos quando 
comparados aos agregados formados por outros mecanismos (Silva Neto et al., 2010), 
consequentemente influenciando nos valores do complexo sortivo desses agregados. 
Para Mergen Junior et al. (2019 b), os agregados biogênicos apresentam-se como sítios 
potenciais de acúmulo de nutrientes, fruto da maior atividade biológica e presença de 
sistema radicular. 

Na Tabela 2 são apresentados os teores de carbono orgânico total (COT), carbono 
orgânico particulado (COP) e carbono orgânico associado aos minerais (COAM), na qual 
não foram verificadas diferenças na comparação entre CX e CC. Em relação as vias 
de formação, foram observadas diferenças somente para COT, sendo nos agregados 
biogênicos em CX quantificados os teores mais elevados de carbono em comparação a 
via fisiogênica. 

Mesmo com a ausência de significância estatística encontrada nos resultados das 
frações COP e COAM nos agregados, pode-se observar tendência de aumento nos teores 
de carbono nessas frações em agregados biogênicos em relação aos fisiogênicos nas 
duas pedoformas avaliadas (Tabela 2). Para COP, o aumento nos teores de carbono nos 
agregados biogênicos foi de 25 e 24% em CX e CC, respectivamente. Já para COAM, o 
aumento foi de 15 e 9% em CX e CC, em sequência. 

Para o conteúdo de COT (Tabela 2), os resultados corroboram aos encontrados por 
Silva Neto et al. (2010), Batista et al. (2013), Loss et al. (2014; 2017), Ventura et al. (2018), 
Melo et al. (2019) e Schultz et al. (2019). Silva Neto et al. (2016) e Fernandes et al. (2017) 
quantificaram maiores teores de COT nos agregados biogênicos quando comparados aos 
fisiogênicos em áreas de floresta secundária em estádios médio (FSEM) e inicial (FSEI) e 
PMM. 

Em área de floresta nativa, no município de Guaíra - PR, Ferreira et al. (2020) 
verificaram maiores teores de COT em agregados biogênicos em relação aos agregados 
intermediários e fisiogênicos. De modo geral, nas áreas sob vegetação nativa, em virtude da 
maior estabilidade do ambiente, pressupõem-se uma maior atividade dos microorganismos, 
plantas e animais edáficos responsáveis pela gênese de agregados biogênicos. Tais 
resultados sugerem que esses agregados, além de indicarem maior atividade biológica, 
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também auxiliam no acúmulo de carbono no solo. 




















COT (g kg!) 

Pedoformas — ——— 

Fisiogênicos Biogênicos 
Cx 36,65 NSb 44,32 NSa 
Cc 38,56 NS ns 43,94 NS ns 
CVi% 12,7 
CV2% 23,4 

COP (g kg *!) COAM (g kg !) 
Pedoformas 

Fisiogênicos Biogênicos Fisiogênicos Biogênicos 

CX 6,33 NS ns 8,51 NS ns 30,23 NS ns 35,77 NS ns 
Cc 5,45 NS ns 7,20 NS ns 33,11 NS ns 36,73 NS ns 
CVI% 37,6 19,9 
CV2% 30,2 30,1 





Médias seguidas pela mesma letra maiúscula na coluna não diferem as pedoformas para o 
mesmo tipo de agregado, mesma letra minúscula na linha não diferem os tipos de agregados 
para a mesma pedoforma avaliado (teste de Tukey, a 5% de probabilidade) e "s Não 
significativo a 5% de probabilidade. COT: Carbono orgânico total; COP: Carbono orgânico 
particulado; COAM: Carbono orgânico associado aos minerais; CV1: Coeficiente de variação 
entre os tipos de agregados; e CV2: Coeficiente de variação entre as pedoformas. 


Tabela 2. Carbono orgânico total e das frações granulométricas da MOS de agregados 
fisiogênicos e biogênicos sob pedoforma convexa (CX) e côncavo (CC) em Floresta Estacional 
Semidecidual Submontana no Sudeste do Brasil, Pinheiral — RJ, na camada de 0,00-0,10 m. 


O aumento de carbono em estruturas biogênicas indica seletiva alimentação do 
material orgânico por organismos do solo, elevando as concentrações de carbono nesses 
agregados (Pulleman et al., 2005), fato este observado com os maiores valores de COT 
quantificados nos agregados biogênicos em CX (Tabela 2). Do exposto, é possível inferir 
que os agregados de origem biológica contribuem mais efetivamente na proteção da MOS, 
diminuindo as taxas de decomposição e elevando o potencial de sequestro de carbono pelo 
solo (Silva Neto et al., 2010). 

Dentre as frações granulométricas, o COP é considerado um dos compartimentos 
mais lábeis, altamente sensível a práticas de revolvimento do solo e facilmente mineralizável 
caso o manejo ou uso não esteja dentro dos preceitos da agricultura conservacionista 
(Rossi et al., 2016), podendo assim ser apontado como um potencial indicador da qualidade 
da MOS em diferentes ambientes edáficos. Todavia em ambientes mais equilibrados e 
estáveis, variações significativas no conteúdo de COP podem ser menos frequentes, como 
observado nesse estudo. 

Já a fração COAM nem sempre funciona como um indicador sensível e eficiente na 
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avaliação da qualidade edáfica, uma vez que alterações no estoque desse compartimento 
da matéria orgânica demoram muitos anos para serem detectados (Carmo et al., 2012), em 
virtude do elevado grau de estabilidade dessa fração física da MOS. 

Os resultados de carbono orgânico das substâncias húmicas da MOS variaram em 
função das pedoformas e das vias de formação em todas as frações estudadas (Tabela 3). 




















C-FAF (g kg ') C-FAH (g kg !) 

Pedoformas E T T O 

Fisiogênicos Biogênicos Fisiogênicos Biogênicos 
Cx 3,20 NSb 513Aa 2,61 NSb 5,05Aa 
Cc 3,44 N$ns 3,29 B's 3,20 NSns 2,53 B's 
CVI% 18,5 26,4 
CV2% 22,7 38,6 

C-FHUM (g kg "') SH (%) 

Pedoformas T T a E 

Fisiogênicos Biogênicos Fisiogênicos Biogênicos 
Cx 12,80Bb 19,93Aa 52NS b 70Aa 
Cc 17,79 Ars 14,62 Brs 64NSa 48Bb 
CVI% 17,5 17,6 
CV2% 18,8 22,9 





Médias seguidas pela mesma letra maiúscula na coluna não diferem as pedoformas para o 
mesmo tipo de agregado, mesma letra minúscula na linha não diferem os tipos de agregados 
para a mesma pedoforma avaliado (teste de Tukey, a 5% de probabilidade) e "* Não 
significativo a 5% de probabilidade. C-FAF: Carbono da fração ácido fúlvico; C-FAH: Carbono 
da fração ácido húmico; C-FHUM: Carbono da fração humina; SH: Porcentagem de carbono 
nas substâncias húmicas; CV1: Coeficiente de variação entre os tipos de agregados; e CV2: 
Coeficiente de variação entre as pedoformas. 


Tabela 3. Carbono orgânico das frações húmicas da MOS de agregados fisiogênicos e 
biogênicos sob pedoforma convexa (CX) e côncavo (CC) em Floresta Estacional Semidecidual 
Submontana no Sudeste do Brasil, Pinheiral — RJ, na camada de 0,00-0,10 m. 


No que tange às pedoformas, em CX nos agregados biogênicos foram quantificados 
os maiores teores de carbono orgânico nas frações ácido fúlvico (C-FAF), ácido húmico 
(C-FAH), humina (C-FHUM) e na porcentagem de carbono nas substâncias húmicas (SH) 
em relação a CC. Entretanto nos agregados fisiogênicos em CC foram verificados os teores 
mais elevados de C-FHUM em comparação a CX. 

Para as vias de formação, foram observados que nos agregados biogênicos 
foram quantificados os maiores teores de C-FAF, C-FAH, C-FHUM e SH em CX quando 
comparados aos fisiogênicos, exceto em CC, na qual nos fisiogênicos foram quantificados 
os maiores valores de SH (Tabela 3). 

De maneira geral, foi verificado que nos agregados biogênicos da pedoforma 
convexa os maiores teores de carbono orgânico das frações húmicas da MOS (C-FAF, 
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C-FAH e C-FHUM) e no total de carbono humificado (SH), tanto entre as pedoformas 
quanto entre as vias de formação, corroborando aos valores de COT observados. Outros 
autores observaram resultados semelhantes para a agregação biogênica (Loss et al., 2014; 
Mergen Junior et al., 2014; Fernandes et al., 2017; Ventura et al., 2018; Melo et al., 2019 e 
Ferreira et al., 2020). Apontando a mesma como uma potencial indicadora na avaliação da 
qualidade do solo, justamente por sua maior sensibilidade as diferentes condições edafo- 
ambientais (Silva Neto et al., 2010; Loss et al., 2014; e Fernandes et al., 2017). 

Em relação as frações C-FAF e C-FAH, Loss et al. (2014), no estado do Paraná, 
em áreas de sistema de preparo convencional (SPC) e plantio direto (SPD), pastagem e 
floresta secundária, verificaram que os agregados biogênicos propiciam um ambiente mais 
favorável para a formação de substâncias húmicas. Os autores atribuem este padrão à 
atuação dos fatores como fauna do solo (macro e microfauna) e sistema radicular. 

Para Ferreira et al. (2020), as frações ácido fúlvico e ácido húmico, são os 
compartimentos mais importantes das frações húmicas da MOS em termos de reatividade 
e ocorrência em ecossistemas. Sendo assim, os elevados teores de C-FAF e C-FAH em 
agregados biogênicos nesse estudo, demonstram a sua importância na qualidade química 
do solo. 

Para o C-FHUM, a maior concentração de carbono nessa fração em agregados 
biogênicos pode ser atribuída a ação das raízes, e principalmente pela fauna edáfica. Nas 
diferentes vias de formação dos agregados, Ferreira et al. (2020) quantificaram maiores 
teores de C-HUM em áreas de floresta nativa e SPD com 23 anos de implantação em 
comparação aos SPDs com 7 e 14 anos de instalação. Para Fernandes et al. (2017), a 
presença da macrofauna edáfica, principalmente as minhocas, têm ação direta sobre a 
ciclagem biogeoquímica no solo (Jouquet et al., 2006), influenciando positivamente a 
ecologia do processo de humificação. 

Os maiores teores de C-FAF, C-FAH, C-FHUM e SH em agregados formados pelos 
agentes biológicos do solo, podem estar relacionados a capacidade que esses agregados 
têm de proteção e estabilização das substâncias húmicas, favorecendo assim o acúmulo 
dessas frações. Silva Neto et al. (2016), sugerem que essa via de formação é mais eficiente 
em termos de estabilidade estrutural do solo e sequestro de carbono. 

A análise de componentes principais (ACP) foi realizada a partir das interpretações 
feitas na distribuição dos valores de DMP, atributos químicos, COT e carbono das frações 
granulométricas e húmicas da MOS nas vias de formação dos agregados (Figura 2). 
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Eixo 2: 31,19% + 





Eixo 1: 50,01%. 


Figura 2. Gráfico bipolar da análise de componentes principais (ACP) da estabilidade, atributos 
químicos e conteúdo de carbono orgânico total e das frações granulométricas e químicas da 
MOS de agregados fisiogênicos (F) e biogênicos (B), sob pedoforma convexa (CX) e côncavo 
(CC) em Floresta Estacional Semidecidual Submontana no Sudeste do Brasil, Pinheiral — RJ, 
na camada de 0,00-0,10 m. 


DMP: Diâmetro médio ponderado; Ca: Cálcio trocável; Mg: Magnésio trocável; Al: Alumínio 
trocável; K: Potássio trocável; P: Fósforo trocável; H+Al: Acidez potencial; S: Soma de bases 
trocáveis; T: Total de cátions adsorvidos, corresponde a CTC à pH 7,0; V: Saturação de bases; 
COT: Carbono orgânico total; COP: Carbono orgânico particulado; COAM: Carbono orgânico 
associado aos minerais; C-FAF: Carbono da fração ácido fúlvico; C-FAH: Carbono da fração 
ácido húmico; C-FHUM: Carbono da fração humina; e %SH: Porcentagem de carbono nas 
substâncias húmicas. 


Para essa análise, os atributos avaliados são representados por setas, que indicam 
a direção do gradiente máximo, em que o comprimento da seta é proporcional à correlação 
do atributo com os eixos e à sua importância na explicação da variância projetada em cada 
eixo. Verificou-se que as distribuições dos atributos apresentaram variância acumulada 
para os eixos 1 e 2 de 81,20%. 

Através da ACP verifica-se a formação de três grupos, sendo os agregados 
biogênicos (B) e fisiogênicos (F) na pedoforma convexa (CX Be CX F, respectivamente) 
separados dos agregados da pedoforma côncava (CC), ao longo do eixo 1, que melhor 
explica a dissimilaridade entre as diferentes condições topográficas e vias de formação 
avaliadas (50,01%). Enquanto o eixo 2, com menos relevância, separou CX B de CX F 
(31,19%). 

Amaioria das variáveis correlacionaram-se positivamente (>0,8) ao eixo 1 e estiveram 
mais associadas aos agregados da pedoforma CX. Ou seja, o Mg?*, K*, COT, COP, COAM, 
C-FAF, C-FAH, C-FHUM e %SH estão relacionados a via biogênica; e DMP, pH, Ca?*, S, 
TeVavia fisiogênica. Verificou-se também que os agregados em CC apresentaram-se 
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próximos dos autovetores Alº* e H+AlI, e distante dos agregados da pedoforma CX. 

Com o emprego da ACP, os diferentes tipos de agregados e pedoformas puderam 
ser separados de acordo com os atributos avaliados. Houve clara separação entre 
as pedoformas CX e CC, e tendo a CC para ambas as vias de formação, correlação 
positiva somente com Alt e H+AI. O agrupamento dos autovetores (Mg?*, Kt, COT, COP, 
COAM, C-FAF, C-FAH, C-FHUM e %SH) próximos aos agregados biogênicos em CX 
demostra que essa condição de relevo, favorece a decomposição do material orgânico. 
Consequentemente aumenta os teores de carbono orgânico, contribuindo na formação de 
agregados biogênicos. 

Segundo Melo et al. (2019), quanto maior a entrada de carbono ao sistema, mais 
intensa será a formação de agregados biogênicos. O padrão observado na APC associado 
aos testes estatísticos aplicados evidenciam o favorecimento da agregação biogênica nas 
diferentes frações da matéria orgânica do solo. 


41 CONCLUSÕES 


As diferentes condições de relevo influenciaram nas vias de formação dos agregados, 
com destaque para a composição química e partição dos teores de carbono orgânico total 
e das substâncias húmicas da MOS encontrados nos mesmos. 

Na pedoforma convexa foram observados os maiores valores de diâmetro médio 
ponderado nos agregados fisiogênicos. 

Nos agregados biogênicos foram verificados melhor qualidade dos atributos 
químicos, representada por maiores valores do complexo sortivo, carbono orgânico total e 
das frações húmicas da MOS em comparação aos agregados fisiogênicos. 

O estudo quantitativo da MOS a partir da compartimentalização do carbono orgânico 
nos diferentes tipos de agregados, pode ser uma técnica eficaz na mensuração da qualidade 
do solo. 
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RESUMO: A curva de retenção de água do 
substrato é importante instrumento na eficiência 
das irrigações. Entretanto, no Brasil, ainda não 
existe um método padrão, confiável, para a 
determinação da curva de retenção de água 
em substratos, destinados para o cultivo 
de ornamentais. Assim, neste trabalho, o 
objetivo foi obter a curva de retenção de água 
em substratos, para o cultivo de ornamental, 
utilizando tensiômetro. O trabalho foi realizado 
no laboratório de hidráulica da Universidade 
Federal do Pampa-Campus de Alegrete/RS. 
Foram testadas diferentes misturas, entre 
solo e substrato da marca Carolina SoilB, 
sendo essas: T1 (100% substrato); T2 (100% 
Solo); T3 (90% solo+10% substrato); T4 (80% 
solo+20% substrato); T5 (70% solo+30% 
substrato); T6 (60% solo+40% substrato); T7 
(50% solo+50% substrato), T8 (40% solo+50% 
substrato); T9 (30% solo+70% substrato) e T1O 
(20% solo+80% substrato). O procedimento 


Engenharia Agronômica: Ambientes Agrícolas e seus Campos de Atuação 3 


adotado foi a determinação conjunta da tensão 
de água nos substratos, pela equação proposta 
por Brito et al. (2009), e a pesagem dos vasos 
para, por diferença, determinar a quantidade de 
água correspondente. Os dados de tensão e 
quantidade de água do substrato foram usados 
para a confecção das curvas de retenção de 
água pelo substrato. Os resultados obtidos 
indicam a possibilidade do uso de tensiômetros 
para determinar as curvas de retenção de água 
nos substratos testados. 

PALAVRAS-CHAVE: Potencial mátrico, água no 
solo, tensiômetro, substratos, irrigação. 


ABSTRACT: The substrate water retention 
curve is an important instrument in the irrigation 
efficiency. However, in Brazil, there is still no 
standard (and reliable) method for determining 
the water retention curve in substrates, intended 
for the cultivation of ornamentals. Thus, in this 
work, the objective was to obtain the water 
retention curve in substrates, for ornamental 
cultivation, using a tensiometer. The work was 
carried out in the hydraulics laboratory of the 
Federal University of Pampa-Campus of Alegrete 
/ RS. Different mixtures were tested, between soil 
and substrate of the Carolina Soil? brand, these 
being: T1 (100% substrate); T2 (100% Solo); 
T3 (90% soil + 10% substrate); T4 (80% soil + 
20% substrate); T5 (70% soil + 30% substrate); 
T6 (60% soil + 40% substrate); T7 (50% soil + 
50% substrate), T8 (40% soil + 50% substrate); 
T9 (30% soil + 70% substrate) and T1O (20% soil 
+ 80% substrate). The adopted procedure was 
the joint determination of the water tension in the 
substrates, by the equation proposed by Brito et 
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al., (2009), and the weighing of the vessels to (by difference) determine the corresponding 
amount of water. The data of tension and quantity of water of the substrate were used for 
the preparation of the water retention curves by the substrate. The results obtained indicate 
that it is possible to use tensiometers to determine the water retention curves in the tested 
substrates. 

KEYWORDS: Potential, soil water, tensiometer, substrates, irrigation. 


11 INTRODUÇÃO 


O solo apresenta gênese, perfil e densidades característicos, com processos de 
formação envolvendo milênios, estando fortemente relacionado com a paisagem e as 
condições ambientais circundantes (Kampf, 2000). É fato a ser considerado, também, que 
as plantas quando cultivadas no solo in situ têm a sua disposição, volumes de solo de 
dimensões praticamente ilimitadas para explorar. 

No cultivo em recipientes, independentes do substrato utilizado, observa-se uma 
limitação de espaço para a expansão das raízes. A limitação deste exige que o substrato 
seja capaz de manter um volume de água facilmente disponível às plantas sem, no entanto, 
comprometer a concentração de oxigênio no meio (Bunt, 1961). 

Conhecer o comportamento da água no substrato é essencial para a tomada de 
decisões referentes ao manejo racional dos cultivos irrigados. Nesse aspecto, uma 
importante variável a ser mensurada é a energia com que a água está retida pela fração 
sólida do substrato, conhecida na literatura que trata do tema por potencial mátrico. Sendo 
chamada de curva de retenção de água no substrato, o gráfico que relaciona o potencial 
mátrico de água no substrato com à umidade a base de volume do mesmo. 

Conhecer somente os valores totais de água contidos em um substrato, ou seja, sua 
capacidade de retenção de água, não é suficiente, pois esta pode apresentar-se retida com 
uma força superior a da sucção exercida pelas plantas, e não estar acessível. 

A água deve estar disponível às plantas em baixas tensões, para que não ocorra 
estresse hídrico ou desvio de energia que seria utilizada para o crescimento e produção 
de biomassa, para a absorção de água. Isso significa que a sucção no substrato não deve 
exceder uma tensão correspondente a 100 cm de coluna de água. De Boodt & Verdonck, 
(1972) estabeleceram como referência os pontos de sucção de 0, 10, 50, e 100 cm de 
coluna de água, na construção da curva de retenção de água. No ponto de sucção zero, o 
substrato encontra-se totalmente saturado, sendo que o volume de água presente indica 
a porosidade total. A sucção de 10 cm foi escolhida por corresponder à altura média da 
camada de substrato em que as raízes das plantas ornamentais se desenvolvem em casa 
de vegetação. O ponto 100 cm foi considerado o limite superior da disponibilidade de água, 
de forma a se evitar limitações sobre o crescimento das plantas (BELLÉ & KÂMPF, 1994). 

O volume liberado entre 10 e 100 cm de sucção é considerado como água disponível, 
sendo este dividido em água facilmente disponível (10 e 50 cm ) e água tamponante ou 
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de reserva (50 e 100 cm). Um substrato ideal possui água disponível e água tamponante 
em torno de 20 a 30%, e 4 a 10%, respectivamente (DE BOODT & VERDONCK, 1972). A 
escolha adequada das matérias-primas e o correto manejo dos substratos é fundamental 
em sistemas de cultivo sem solo, pois diferentes substratos apresentam 17 diferentes 
proporções entre a água fracamente e a fortemente retida pelo meio, o que é de extrema 
importância em um programa de irrigação (MARTINEZ, 1992). Portanto, o fornecimento de 
água às culturas deve ocorrer com intensidade e frequência apropriados para cada tipo de 
substrato, de forma a maximizar o crescimento e desenvolvimento das plantas. 

Na determinação da curva de retenção de água em substratos agrícolas, entre 
outros, têm sido utilizados métodos que envolvem equipamentos de pressão e placas 
porosas. Também se utilizam funis de “Buchener”, equipados com placas porosas ou outros 
materiais, por exemplo, gesso. Uma das limitações desses métodos, tema de discussões, 
é o problema do contato da amostra com a placa porosa e com outros tipos de materiais. 

Vieira e Castro (1987) afirmam que em relação a outros métodos o uso de cápsulas 
porosas para determinação de curvas de retenção de água em componentes é bastante 
vantajoso, tanto em termos de equipamentos necessários, tempo de obtenção de pontos, 
número de pontos possíveis e, principalmente, quanto à proximidade de curvas determinadas 
por outros métodos. Esses mesmos autores determinaram curvas de retenção de água 
em componentes bastante contrastantes em granulometria, usando cápsulas porosas de 
10 mm de diâmetro. Os resultados que obtiveram foram bons e comparáveis a métodos 
tradicionais, como o da câmara de Richards. 

Segundo Vieira e Castro (1987), com a utilização do método da cápsula porosa foram 
necessários cerca de 40 dias para cada curva de retenção enquanto que a determinação 
em câmara de Richards se estendeu por 150 dias. Estes, relataram que quando usaram 
cápsula porosa de 10 mm na determinação da curva de retenção de água em Latossolo 
Roxo Distrófico, Podzólico Vermelho Amarelo de textura média argilosa e de textura 
arenosa média também não tiveram problemas com o contato entre as amostras de solo e 
os meios de tensão. 

No entanto, em relação à curva de retenção de água, percebe-se que existe 
variabilidade nos resultados encontrados por diferentes autores. Possivelmente porque 
vários são os fatores que influenciam na determinação dessa variável física em substratos 
agrícolas, como por exemplo o equipamento utilizado nas determinações. 

No Brasil, a determinação da curva de retenção de água em substratos agrícolas 
ainda não é padronizada e os métodos mais utilizados demandam muito tempo, são 
tediosos e necessitam de equipamentos caros. 

Neste contexto, esse trabalho teve por objetivo mensurar a tensão de água em 
diferentes substratos, por meio do uso de tensiometria, para determinar a curva de retenção 


de água nos mesmos. 
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21 MATERIAL E MÉTODOS 


O trabalho foi desenvolvido no laboratório de hidráulica da Universidade Federal do 
Pampa/Campus de Alegrete-RS, com solo coletado na localidade do Capivari, no referido 
município, e com substrato comercial da marca Carolina Soil. 

O solo utilizado possui classificação de Argissolo (Streck et al., 2008). A amostragem 
do mesmo foi feita por meio da abertura de uma mini trincheira, na qual retirou-se uma 
amostra deformada de solo na camada de 0,00-0,20 m. de profundidade. As amostras de 
solo e substrato foram passadas em peneira de 2 mm e seca ao ar, caracterizando terra 
fina seca ao ar. 

Foram testadas diferentes misturas de solo com substrato comercial. As misturas 
testadas foram: S1 (100% substrato); S2 (100% solo); S3 (90% solo + 10% substrato); S4 
(80% solo + 20% substrato); S5 (70% solo+30% substrato); S6 (60% + 40% de substrato); 
S7 (50% solo + 50% substrato); S8 (40% solo + 60% substrato); S9 (30% solo+70% 
substrato) e S10 (20% solo+80% substrato). Estas misturas foram alocadas em vasos de 
plásticos na cor preta com diâmetro de 0,3 m, altura de 0,2 m e capacidade de 2,6 L. 

Cada recipiente foi preenchido com a mistura seca ao ar, até um limite, de modo que 
uma borda do vaso ficou acima do nível do solo. A mistura foi colocada suavemente solto no 
vaso em camadas de 0,05m, onde uma leve compressão foi aplicada. Em seguida os vasos 
foram postos para saturar em água no período de 24 horas, logo após foram instalados os 
tensiômetros. 

Os tensiômetros foram preparados para instalação, segundo a metodologia da 
Embrapa (1999). Estes foram instalados no centro do vaso a 5 cm de distância do fundo do 
recipiente plástico, ficando com 12 cm de profundidade dentro do solo. 

Os vasos, recipiente contendo a mistura e o tensiômetro, foram colocados sobre 
uma bancada no laboratório e pesados diariamente, em horário fixo no turno da manhã, 
até que a variação de massa dos mesmos fosse insignificante, utilizando-se uma balança 
eletrônica, com precisão de 0,019. Concomitante a leitura da massa era realizada a leitura 
do tensiômetro. O experimento teve duração de 70 dias. 

Após encerradas as pesagens e as leituras de potencial de água no solo (kPa). 
O solo foi retirado dos recipientes, levado para secagem em estufa a 105 “C por 24 he 
pesado novamente para determinação da massa de solo seco. 

A massa de água presente no solo em cada pesagem foi determinada pela diferença 
entre as leituras diárias da balança e a massa do vaso (recipiente contendo a mistura seca 
e tensiômetro). A umidade volumétrica do solo foi obtida pela relação entre o volume de 
água em cada pesagem (considerando a densidade da água de 1000 kg.m?) e o volume de 
mistura em cada recipiente. 

A determinação do potencial mátrico da água no solo para os tensiômetros, com 
vacuômetro de Bourdon, foi realizada, pela equação 1, proposta por Brito et al., (2009): 
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bm = (—(L x 0,0136) =h +, + 2) x 9,8066 (1) 


Em que: 

&,= potencial mátrico da água no solo (kPa); L =leitura no visor analógico (mmHg); 
h = altura da coluna de água acima do ponto de inserção do manômetro (mca); h, = altura 
de inserção, em relação à superfície do solo, do manômetro de Bourdon no tubo de PVC 
(mca); Z = distância entre o centro da cápsula e a superfície do solo (mca). 

Após a obtenção da umidade do solo com os respectivos valores de tensão, ajustou- 
se os dados por meio do modelo de Van Genuchten (1980), a partir da equação 2. 

0 = 6,+ TT (2) 

Em que: 

O = a umidade volumétrica (mº*.m*); 9r = umidade residual (mº.m**); Os = umidade de 
saturação (mº.m*); 4, = potencial mátrico (kPa); a = escalonador do p ; me n = parâmeiros 
do modelo. 

O ajuste da curva de retenção de água no solo foi realizado por meio do software 
Soil Water Retention Curve — SWRC (Dourado Neto et al., 1990), optando-se pelo modelo 
de Van Genuchten (1980). 

A partir dos valores de umidade volumétrica, obtidos pelo modelo de Van Genuchten, 
obtiveram-se as variáveis água facilmente disponível (AFD), água tamponante (AT) e água 
residual (AR), onde o valor de AR será igual ao obtido pela umidade volumétrica a 10 
kPa, AFD equivale a subtração de umidade volumétrica na tensão 1kPa pela umidade 
volumétrica na tensão 5kPa e o valor de AT obtém-se a partir da redução do valor da 
umidade volumétrica na tensão 5kPa pela umidade volumétrica em 10kPa, conforme De 
Boodt e Verdonck (1972) e Wilson (1983). 


31 RESULTADOS E DISCUSSÃO 


Os parâmetros de ajuste dos pares de dados de umidade e tensão ao modelo de 
Van Genuchten (1980), para as diferentes misturas, estão apresentados na Tabela 1. 
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Mistura [si 9** a n m 2 
s 


r 





S1** 0,243 0,000 0,009 0,693 0,009 0,947 
S2 0,199 0,025 0,089 29,670 0,032 0,954 
S3 0,223 0,000 0,230 2,726 0,184 0,963 
s4 0,230 0,000 0,158 11,039 0,052 0,981 
sS5 1,019 0,023 0,238 4,603 0,261 0,966 
s6 1,019 0,003 1,025 2,674 0,262 0,989 
S7 1,056 0,023 0,507 3,587 0,260 0,971 
s8 0,256 0,000 0,141 0,862 0,705 0,945 
S9 0,265 0,042 0,161 10,920 0,081 0,943 
S10 0,253 0,000 0,001 0,485 6,993 0,824 





*8, — umidade residual; **9, — umidade de saturação; ***S1 - 100% substrato); S2 - 100% solo; 
S3 - 90% solo + 10% substrato; S4 - 80% solo + 20% substrato; S5 - 70% solo+30% substrato; 
S6 - 60% solo + 40% de substrato; S7 - 50% solo + 50% substrato; S8 - 40% solo + 60% 
substrato; S9 - 30% solo+70% substrato e S10 - 20% solo+80% substrato. 


Tabela 1. Parâmetros de ajuste ao modelo de Van Genuchten (1980) para curvas de retenção 
de água no solo, elaboradas para diferentes tipos de misturas de solo e substrato. 


Observa-se, na Tabela 1, que o ajuste obtido (r?) para a equação de Van Genuchten 
(1980) foi próximo a 1, em todas as misturas estudadas. Segundo GARCIA (1989), valores 
de coeficiente de determinação superiores a 0,70 indicam bons ajustes, ou seja, com estes 
resultados, pode-se usar o modelo testado para estimar a umidade volumétrica destas 
misturas, para potenciais matriciais entre O e 80 kPa. 

Não observa-se relação entre umidade residual com as diferentes porcentagens 
de solo e substrato, estes resultados corroboram com os encontrados por Gonçalves 
(2012). Esta autora estudou curvas de retenção e condutividades hidráulicas, publicadas 
na literatura, para diferentes tipos de solos e quantificou a variação dos parâmetros de Van 
Genuchten (1980) com a distribuição granulométrica desses solos. Ainda, Schaap et al. 
(2001) assumiu que 6, é zero para uma grande quantidade das classes texturais estudadas, 
para solos sedimentares da América do Norte e Europa. 

Os maiores valores de umidade de saturação foram observados nas misturas que 
contém em sua composição valores entre 50 a 70% de solo com porcentagens entre 30 a 
50% de substrato comercial. Isso, segundo Silva (2005), do ponto de vista capilar, significa 
a presença de poros menores. 

As curvas de retenção de água, para as diferentes misturas, obtidas com 
tensiômetros, no laboratório, e determinadas por secamento, são apresentadas na Figura 
1.Nota-se, que os tensiômetros foram eficientes para a obtenção da curva de retenção de 
água, para as diferentes misturas. Conseguindo-se obter a umidade volumétrica para o 
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potencial matricial de até aproximadamente 80 kPa, exceto para a mistura composta por 
20% solo+80% substrato. Observando-se, assim, a eficiência do uso de tensiômetros, na 
estimativa da curva de retenção de água em substratos, uma vez que a faixa de tensão 
de trabalho, dos tensiômetros, para maioria das culturas fica entre O até 75 kPa (LIBARDI, 
1999). 

Nos maiores potenciais matriciais ocorreram as maiores variações de umidade 
volumétrica. Resultado semelhante foi observado por Farias, et al. (2004), estes para 
substrato constituído de 30% de solo+30% de pó de xaxim+40% de casca de pínus, 
registraram uma redução de aproximadamente 3% da umidade volumétrica nos potenciais 
inferiores a 10 kPa, com posterior tendência a estabilização, da umidade, até o potencial 
matricial de 100 kPa. 

Segundo Diaz (2010) a sucção no substrato não deve exceder uma tensão 
correspondente a 100 cm de coluna de água (9,8 kPa). De Boodt & Verdonck, (1972) 
estabeleceram como referência os pontos de sucção de 0, 10, 50, e 100 cm de coluna de 
água, na construção da curva de retenção de água, para substratos. 

Observa-se, ainda, na Figura 1 que a umidade volumétrica quando o potencial 
matricial estava próximo de zero, tensão inicial, variou de 0,223 mº.mº na mistura S3 a 
0,256 m*.m* na composição S8. Na mistura com 100% de solo, na tensão de 78,7 kPa, foi 
registrado o menor valor de umidade, 0,053m'.m*. 

Fermino (2003) determinado a curva de retenção, através do uso de mesa de tensão, 
funil de tensão e cilindros de pressão, para diferentes substratos comerciais, encontrou 
umidade volumétrica média de 0,6 m*.mº, na tensão inicial de 10 hPa, e de 0,34 mº.mº, 
para a tensão final de 100 hPa. 
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Figura 1. Curva de retenção de água para os substratos preparados com 100% substrato (A), 
100% solo (B), 90% solo + 10% substrato (C); 80% solo + 20% substrato (D); 70% solo+30% 
substrato (E); 60% solo + 40% de substrato (F); 50% solo + 50% substrato (G); 40% solo + 
60% substrato (H); 30% solo+70% substrato (I) e 20% solo+80% substrato (J). 


Segundo Andriesse (1988) a retenção de água das turfas pode chegar até 20 vezes a 
sua massa. Também as cascas, como as de Pinus (Handreck e Black, 1999) e a vermiculita 
(Burés, 1997), têm alta capacidade de retenção de água. 

O volume de água retida pelos substratos em diferentes tensões é apresentado na 
Tabela 2, observa-se que apenas os substratos S5, S6 e S7, têm valores de água facilmente 
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disponível (AFD), entre 26,16 e 58,01%, superiores ao de referência (25% de AFD) segundo 
De Boodt et al. (1974). Assim, os substratos com valores de AFD inferiores ao de referência 
não apresentam condições satisfatórias para o desenvolvimento das plantas, uma vez que 
segundo Klein et al. (2000), o importante para o bom desenvolvimento das plantas é a água 
facilmente disponível. 

Manfron et al. (2005), avaliando a disponibilidade de água, que o substrato composto 
por 60% de húmus e 40% de casca de arroz natural, disponibiliza para as plantas de alface, 
observou valores de água facilmente disponível, inferior a faixa ideal, em média de 10,46% 








Substrato AR (%) AFD (%) AT (%) 
Si 14,09 4,43 3,07 
s2 19,88 0,00 0,02 
s3 14,43 3,33 4,47 
s4 17,72 0,08 5,19 
s5 37,27 26,16 38,43 
s6 20,24 51,72 12,35 
S7 25,07 58,01 20,30 
s8 14,05 5,39 3,29 
s9 18,80 0,16 7,53 

Sto 12,86 4,63 2,64 





Tabela 2. Água residual (AR), água facilmente disponível (AFD) e água tamponante (AT) as 
plantas para os substratos avaliados nas diferentes tensões. 


Azevedo et al. (2014) caracterizando as propriedades físico-hídricas de diferentes 
substratos, observaram que os substratos de casca de pinus e misto artesanal, apresentaram 
valores de AFD abaixo do considerado satisfatório por Fermino et al. (2014) e Cadahia 
(1998), com valores de 12,21 e 15,50%, respectivamente. 

Costa et al. (2017) em estudo, com distintos substratos, registraram, para a água 
facilmente disponível, um modelo quadrático, com o ponto de mínima de 0,6% com a 
utilização de 37,23% de casca de arroz carbonizada. Estes, observaram apenas para o 
substrato formado pela mistura de 75% de casca de arroz carbonizada + 25% composto 
orgânico, valores compatíveis com os de referência, proposto por Fermino (2003) e Cadahia 
(1998). 

De acordo com Costa et al. (2017) a AFD é liberada sob baixas tensões, indicando 
que é rapidamente perdida, não se constituindo um reservatório de água para as plantas. 
Desta forma, quando no substrato há grande percentual de AFD, faz-se necessário o 
prolongamento nas regas. Para Zorzeto et al. (2014) valores críticos, ou seja, abaixo do 
considerado ideal, podem inibir o crescimento vegetal na falta de irrigações frequentes. 
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Ainda, Schafer (2004) observou que substratos com valores abaixo do recomendado, 
tanto para AFD como para a AT, apresentam limitações quanto ao suprimento adequado de 
água de fácil absorção para as plantas, indicando que, as regas devem ter durações mais 
curtas com intervalos de tempo menores. 

Avaliando os valores de água residual (AR), que é a água que não será disponibilizada 
às plantas, o maior e menor valor foram observados nos substratos S5 (37,27 %) e S10 
(12,86%), respectivamente (Tabela 2), ou seja, esses percentuais de água não estão 
disponíveis para absorção das raízes. Esses resultados corroboram com os observados por 
Feraz et al. (2005), estes, trabalhando com diferentes substratos comerciais, encontraram 
valores de AR variando de 14 % para o substrato Germina (10%) e 26 % para o substrato 
de turfa. 

Azevedo et al. (2014) encontraram valores para AR variando de 36,37 a 86,68 %, 
para fibra de coco e vermiculita, respectivamente, em estudo das características físico- 
hídricas de diferentes substratos encontrados para comercialização, 

Costa et al. (2017) observou, para a variável AR, comportamento linear decrescente 
e as composições S3 (50% CAC + 50% CO) e S4 (75% CAC+ 25% CO) atingiram valores 
de 31,09% e 19,95%, respectivamente. 

Nota-se, na Tabela 2, para água tamponante (AT) que os substratos S3, S4, S5, 
S6, S7 e 89, apresentaram valores acima do valor ideal. De Boot & Verdonck (1972) e 
Haynes & Gohg (1978) estimaram que entre 4 a 10% seria um valor de referência para 
água tamponante. De acordo com os autores substratos com valores iguais ou próximos ao 
recomendado proporcionam suprimento hídrico adequado para as plantas. Já, os demais 
substratos apresentam valores inferiores a referência, portanto estes não proporcionariam 
um possível suprimento hídrico as plantas nele estabelecidas caso haja stress hídrico. Para 
os mesmos autores o volume liberado acima desta tensão é considerado água dificilmente 
disponível, remanescente ou residual, portanto, desconsiderada. 

Estes resultados são semelhantes aos observados por Costa et. al. (2014), os 
mesmos observaram valores de AT inferior a 5% para os substratos compostos por: 100% 
composto orgânico; 25% casca de arroz carbonizada + 75% composto orgânico; 50% casca 
de arroz carbonizada + 50% composto orgânico; 75% casca de arroz carbonizada + 25% 
composto orgânico e 100% casca de arroz carbonizada. Manfron et al. (2005) observou 
para as amostras de substratos volume de água tamponante em média de 1,30%. 

As diferenças de resultados encontrados no trabalho e na literatura, certamente 
ocorrem devido a grande variabilidade existente entre os substratos, nas suas composições 
e estruturas, bem como, no arranjamento de suas partículas no decorrer do tempo, 


resultando em efeitos diversos para as plantas. 
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41 CONCLUSÃO 


A determinação da curva de retenção de água de substratos através do método 
do tensiôimetro, pelos resultados deste estudo, apresenta-se como uma alternativa 
potencialmente viável aos métodos tradicionais, da mesa de tensão e da câmara de 
pressão de Richards. 

O substrato composto por 70% solo e 30% substrato apresenta a melhor qualidade 
físico-hidrica, sendo uma interessante alternativa para produção de flores ornamentais com 
qualidade e baixo custo. 

As diferentes composições dos substratos revelam o contraste do comportamento 
físico-hídrico dos mesmos, mostrando a importância de sua caracterização, podendo assim 
definir a frequência e quantidade das irrigações. 
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